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FOREWORD 


As we enter the sixth year of the program to improve the structural integ- 
rity and durability of reusable space propulsion systems, it is time to reflect 
on where we were, where we are, and where we hope to be. 

The program was initiated in late 1983 when the NASA Office of Space Flight 
(OSF) recognized a deficiency in the space shuttle main engine (SSME) technology 
base and decided to augment the Ear th-to-Orbi t Propulsion Technology Program 
sponsored by the Office of Aeronautics and Space Technology (OAST). That defi- 
ciency was the inability to predict the life of structural components in high- 
performance reusable rocket engines. It was manifested by the SSME component re 
placement and repair requirements. Development of the analytical models for 
aerothermodynamic loads, structural response, and fatigue/fracture effects was 
therefore undertaken to produce more accurate life prediction design tools. 

Today, several hundred government, industry, and university technologists 
are involved in establishing the reusable rocket engine design data base. Pro- 
ducts of their involvement can be seen in the designs of the advanced launch sys 
tern (ALS), heavy lift launch vehicle ( HLLV ) , and chemical transfer or orbit 
transfer vehicle. Developments of a three-dimensional viscous code for multi- 
stage turbines, of probabilistic structural analysis codes, and of the under- 
standing of the interaction between high and low cycle fatigue are examples of 
the strides made to date. Advanced instrumentation techniques have been proceed 
ing in parallel to enable validation/calibration of those codes in technology 
test-bed applications. Fiber optic systems, heat flux gauges, and thin-film sen 
sors developed under this program are now finding their way into propulsion test 
facilities * 

Life prediction algorithms will be utilized in health management and diag- 
nostic expert systems to reduce ground testing maintenance operations and turn- 
around times. As confidence grows in accuracy, flight systems will be developed 
Eventually, these expert systems will be allowed to "close the loop" and provide 
real-time engine control to prevent operational failures and guarantee mission 
success . 

This third biannual conference provided a forum for representatives from 
government, industry, and academia to learn about and discuss the latest find- 
ings and progress. Extended abstracts and the figures used for each of the pre- 
sentations at the two-day conference are included in these proceedings. 


Sol Gorland 
Conference Chairman 
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PREFACE 


A major part of the NASA technology program is focused on advanced propul- 
sion systems designed to meet the future Earth-to-orbit (ETO) space transporta- 
tion needs of the agency. The objective of this program is to provide the tech- 
nology base that will allow the design, development, and operation of advanced 
propulsion systems with pre-defined, predictable performance, service life, and 
operational characteristics needed to significantly reduce the life-cycle-costs 
of the ETO vehicles in which they are to be used, whether these vehicles are ex- 
pendable, partially or fully reusable, manned or unmanned. 

A major objective of the ETO propulsion program is to signi f icantly enhance 
knowledge and understanding of fundamental rocket engine chemical and physical 
processes. From this understanding will come the ability to analytically simu- 
late internal engine environments; to define steady and unsteady thermal and me- 
chanical loads that these environments impose on critical engine components, as 
well as the resulting component material behavior and structural response; and 
to define failure mechanisms and to develop accurate component life predictions. 
With these tools in hand, means for making internal environments less hostile, 
as well as for developing design methodologies that will result in more load re- 
sistant component designs, will be explored, leading to greatly increased design 
margins and longer component life and/or increased component performance. An- 
other major objective of the program is to develop technologies for reducing 
ground operations costs and for ensuring safer, more reliable flight operations 
through the development and demons tr ation of on-board integrated condition and 
safety monitoring systems. 

The Structural Integrity and Durability Program is a key part of the first 
objective and is currently focused on the hot gas sections of typical rocket en- 
gines, particularly the hot gas flow through the turbine stages. Specific areas 
of effort, the status of which will be reviewed during this conference include 
hot gas flow characterization (CFD), aerothermo loads definition, material behav- 
ior and structural analyses (including probabilistic theories), fracture mechan- 
ics, high cycle/low cycle fatigue life prediction, and research ins trumentation. 

The Space Shuttle Main Engine (SSME) development program has provided a 
unique opportunity to identify technology deficiencies and to initiate the pur- 
suit of technology solutions to overcome these deficiencies. These advanced 
technologies will be used primarily for the design and development of advanced 
ETO engines of the future, although the SSME has already benefited to a signifi- 
cant degree from this technology focus. The SSME will also be used as a testbed 
to verify and validate many of the resulting analytical techniques and design 
methodologies, as will the large scale combustor and turbomachinery subsystems 
testbeds that are now an integral part of the ETO propulsion technology program. 


Earl VanLandingham 
NASA Headquarters 
Washington, D.C. 
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OVERVIEW OF THE FATIGUE/FRACTURE/LIFE WORKING GROUP PROGRAM 

AT THE LEWIS RESEARCH CENTER 

Michael A. McGaw 
NASA Lewis Research Center 
Cleveland, Ohio 


The objective of the Fatigue/Fracture/Life Working Group at Lewis is to 
develop and verify constitutive and life prediction models for materials typi- 
cally used in hot gas path components of reusable space propulsion systems 
over the range of relevant operating environments . The efforts at Lewis have 
centered on the development of crack initiation life prediction methods, while 
the efforts of our counterpart group located at the Marshall Space Flight Cen- 
ter have centered on the development of cyclic crack propagation life predic- 
tion methods. 

The complexion of active tasks is shown in figure 1. The program is a 
blend of in-house, contract, and grant research. Significant progress has 
been made during the past year; the papers presented in this session represent 
contributions from most, but not all, of the tasks currently active in this 
program, 

A significant new task started this year will incorporate the various 
material constitutive and life prediction models (including materials proper- 
ties data) developed in this program into a comprehensive creep-fatigue damage 
analysis and life assessment computer code (fig. 2). The program will func- 
tion as a postprocessor to general structural analysis programs (such as fi- 
nite element or boundary element codes) using the output of such analyses 
(stress, strain, and temperature fields as functions of time) as the input to 
the damage analysis and life assessment code. The code will be designed to ex- 
ecute on engineering/scientific workstations and will feature a windowing, 
mouse-driven user interface. Current plans call for the code to be finished 
and made available for use in mid 1991. 
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CURRENT STRUCTURAL DURABILITY 
PROGRAMS AND GOALS 

MACROSCOPIC EFFECTS OF HIGH PRESSURE HYDROGEN ON 
CONSTITUTIVE AND LIFE MODELS 

• ESTABLISH VALIDITY OF STRUCTURAL DURABILITY MODELS 
DEVELOPED UNDER THE PROGRAM 

HCF/LCF INTERACTION UNDER MISSION-RELATED LOADING 

• INTERSPERSED HCF AND LCF LOADINGS 

• TMF INTERACTION WITH HCF LOADINGS 

• CUMULATIVE CREEP-FATIGUE INTERACTIONS , 

APPLICATION OF CONSTITUTIVE AND LIFE MODELS IN STRUCTURAL 
ANALYSIS 

• APPLY MODELS DEVELOPED TO SSME COMPONENTS 
THERMAL RATCHETING DAMAGE ANALYSES 

• INTERACTION OF RATCHET STRAINS WITH TMF AND HCF 

• APPLICATION TO MAIN COMBUSTION CHAMBER MATERIALS 

CD-89-39970 


CURRENT STRUCTURAL DURABILITY 
PROGRAMS AND GOALS 


ADVANCED CONCEPT MATERIALS MODELS DEVELOPMENT 

• METAL MATRIX COMPOSITES MODELING: MATERIAL 
CONSTITUTIVE AND LIFE MODELS 

CYCLIC CRACK GROWTH UNDER THERMOMECHANICAL LOADING 
CONDITIONS 

• APPLICATION OF PATH INDEPENDENT INTEGRALS TO THE 
DESCRIPTION OF THE CRACK DRIVING FUNCTION 

SINGLE CRYSTAL CONSTITUTIVE AND LIFE MODEL DEVELOPMENT 

• ANISOTROPIC CONSTITUTIVE AND LIFE MODELS 

• ANALYSIS METHODOLOGY FOR ADDRESSING BLADE ROOT 
ATTACHMENTS 

CREEP-FATIGUE DAMAGE ANALYSIS PROCESSOR 

• INTEGRATE CONSTITUTIVE AND LIFE PREDICTION MODELS, AS 
WELL AS RELEVANT MATERIALS PROPERTIES DATA INTO A 
DAMAGE ANALYSIS AND LIFE ASSESSMENT COMPUTER CODE 


CD-89- 3997 1 
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CREEP-FATIGUE DAMAGE ANALYSIS PROCESSOR 


• COMPREHENSIVE COMPONENT DAMAGE ANALYSIS AND LIFE 
ASSESSMENT SYSTEM 

• MOUSE-DRIVEN WINDOWING USER INTERFACE IN A PORTABLE 
SOFTWARE PACKAGE 


MATERIALS PROPERTIES DATABASE 

1 

> , , : 

1 I" 

1 
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USER INPUT 


CO-89-39972 

Figure 2 
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STRUCTURAL 
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FAILURE 

MODELS 
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LIFE 

ASSESSMENT 
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application of cyclic damage accumulation life prediction model 

TO HIGH TEMPERATURE COMPONENTS 


Richard S. Nelson 


/ 935*7 


United Technologies Corporation 
Pratt & Whitney 
E. Hartford, Connecticut 


INTRODUCTION 

A high temperature, low cycle fatigue life prediction method has been 
developed by Pratt & Whitney under the sponsorship of the Lewis Research 
Center (Moreno et al, 1984). This method, Cyclic Damage Accumulation, has 
been developed for use in predicting the crack initiation lifetime of gas 
turbipe engine materials, but it can be applied to other materials as 
well. The method is designed to account for the effects on creep-fatigue 
life of complex loadings such as thermomechanical fatigue, hold periods, 
waveshapes, mean stresses, mul t iaxial i ty, cumulative damage, coatings, and 
environmental attack (Nelson et al, 1986). Several features of the model 
were developed to make it practical for application to actual component 
analysis, such as the ability to handle non-isothermal loadings (including 
TMF), arbitrary cycle paths, and multiple damage modes. 

The CDA life prediction model was derived from extensive specimen tests 
conducted as part of the contract activity on the cast nickel-base 
superalloy B1900+Hf. These included both monotonic tests (tensile and 
creep) and strain-controlled fatigue experiments (uniaxial, biaxial, TMF, 
mixed creep-fatigue, and controlled mean stress). Additional specimen 
tests were conducted on wrought INCO 718 to verify the applicability of 
the final CDA model to other high-temperature alloys. The model will be 
available to potential users in the near future in the form of a 
FORTRAN-77 computer program. 

FEATURES OF CDA LIFE PREDICTION MODEL 

The CDA model is based on the fundamental assumption that an 
instantaneous damage rate can be calculated and integrated on a 
cycle-by-cycle basis until a damage of unity is achieved. The current 
constants in the program reflect the definition of damage as the 
initiation of detectable cracks (0.030 in. {0.76 mm.}) at the location 
being considered. Other definitions of damage may also be modeled by 
recalibration of the various constants. Different rates are calculated 
simultaneously for whatever modes of damage may apply to a particular 
material. Currently defined modes in the CDA program include transgranular 
and intergranular fatigue, oxidation, and coating cracking. 

One of the most significant features of the CDA model is its use of 
ratios of stress/strain parameters rather than absolute levels. This 
technique greatly simplifies practical application of the model to actual 


Work performed under NASA Contract NAS3-23288; Dr. Gary R. Halford, LeRC, 
serves as NASA Technical Monitor. 
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components, since these ratios are generally easier to calculate than the 
absolute levels* The model also makes use of total strain rather than 
inelastic strain; this avoids the problem of having to calculate the very 
small inelastic strains associated with the typically high design lives 
for components made of modern superalloys. The t ransgranular damage model 
incorporates a loading history term derived from primary creep data. This 
captures the effects of overloads and other previous load excursions which 
can set up compliant dislocation structures and thereby affect subsequent 
fatigue behavior. Mean stress effects are handled through an exponential 
term based on the ratio of maximum cycle stress to a reference maximum 
stress. Subcycles are determined by the computer program using the 
traditional rainflow cycle counting method, based on reversals of stress 
range ratio rather than stress or strain alone. 

Temperature effects are very important to any high temperature model, 
and these are included in the CDA transgranular model in several ways. 
First, time dependence (hold time, strain rate, etc.) is introduced 
through the use of a multiplier term based on the integral of an Arrhenius 
function. This will automatically handle any arbitrary periodic waveform, 
eliminating the need to estimate average strain rate, frequency, or other 
such parameters from a component loading history. Second, an additional 
strain term is postulated to result from microscopic thermal expansion 
mismatch within the material structure. This is added to the macroscopic 
total strain, and the resulting strain range is used by CDA to compute the 
life. Finally, temperature effects are included implicitly through the 
temperature dependence of certain reference stress/strain quantities. 

The second major damage mode currently active in the CDA model is the 
intergranular mode. This is based on an integral of an Arrhenius function 
which is calibrated to data from creep testing. This was found adequate 
for many components which fail simply from cyclic loadings in the creep 
regime. Other damage modes include cyclic oxidation (which modifies the 
transgranular damage term) and coating cracking (based on the coating life 
model from a companion NASA contract {Swanson et al, 1987}). 

APPLICATION TO HIGH TEMPERATURE COMPONENTS 

The CDA computer program is designed to accept one or more "defined 
cycles" for which a complete stress/strain/temperature/time history is 
known. These may be obtained from several sources, including simplified 
methods or full inelastic analysis programs. The program accepts this data 
(interpolating where required) and directly integrates the rates of the 
various damage modes until one of them reaches unity. 

REFERENCES 

Moreno, V., Nissley, D. M. , and Lin, L. S., 1984, "Creep Fatigue Life 
Prediction for Engine Hot Section Materials (Isotropic) Second Annual 
Report," NASA CR-174844. 

Nelson, R. S., Schoendorf, J. F., and Lin, L. S., 1986, "Creep Fatigue 
Life Prediction for Engine Hot Section Materials (Isotropic) Interim 
Report," NASA CR-179550. 

Swanson, G. A., Linask, I., Nissley, D. M., Norris, P. P., Meyer, T. 

G., and Walker, K. P., 1987, "Life Prediction and Constitutive Models for 
Engine Hot Section Anisotropic Materials Program," NASA CR-179594. 
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OVERVIEW 


• Summary of Contract Activities 

• Features of CDA Life Prediction Model 

• Application to High Temperature Components 

• Future Research 


THE CDA METHOD IS BASED ON EXTENSIVE HIGH 
TEMPERATURE SPECIMEN TESTING 
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CDA MODEL BEGINS WITH INTEGRATION OF DAMAGE 
RATES FOR MULTIPLE MODES 


For use as a practical design method, “Damage = 1” is defined as initiation of 
a 0.030 in. (0.76 mm.) surface length crack 


• The basic equation for CDA is simply an integration of damage rate from 
0 to 1: 


' -/ 


Ni dD dN 
1 dN 


0 ) 


• This same calculation is performed simultaneously for multiple modes of 
damage, as required for a given material. Currently defined modes 
include: 


1. Transgranular fatigue 

2. Intergranular fatigue 

3. Oxidation / environmental damage 

4. Interactions of the above 


TRANSGRANULAR DAMAGE MODEL UTILIZES A 
STRESS AND STRAIN RATIO CONCEPT 


Absolute inelastic strains can be small and dihicuH to predict reliably for many 
nickel-base superalloys 

• The damage rate is calculated using the ratio of current parameters to 
reference values: 


— = 

( 1 ) 

/ AtT 

1 / A o Fi: 

- / epref \ 

dN 

\ Nref f 

' Aet, REF ' 

\AorefZ 

' £P ) 


• A loading history parameter derived from primary creep data is also 
included: 

MM = f £ p 

VERREF/ 'OERREF / MAX, HISTORY (3) 
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TWO ADDITIONAL TERMS WERE NECESSARY TO 
PREDICT BASELINE FAST RATE DATA 


• The first predicts mean stress and R-ratio effects using a ratio based on 
maximum tensile stress: 


dD 

dN basic 


dD 

dN 


/ °MA X \ 
6 '°MAX,REF/ 


(4) 


• A threshold strain range was necessary to predict the high life, low strain 
data for both B1900+Hf and INCO 718: 


Aej = ^tactual -A£ tthreshold 


(5) 


STRAIN RATE EFFECTS CAN BE PREDICTED FOR ANY 
ARBITRARY WAVEFORM 


Permits direct application of CDA to practical component cycles 


• The basic transgranular damage rate (product of the 
modified by a time & temperature dependent term: 

fTD = 1 + / Aie B1o e (m-) dt 


previous terms) is 


(6) 


• This eliminates the need to estimate hold time, average frequency, or 
other similar parameters for complex loading cycles typically experienced 
by modern high temperature components 
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CLASSICAL CUMULATIVF DAMAGE EFFECTS ARE ALSO 
PREDICTED BY THE CDA MODEL 


Real components see a spectrum of loads during their service life 

• Non-linear damage accumulation is included in a modifying function for 
the basic transgranular damage rate: 

G(D) - (l-f L )(/M I)D' J + f L (7) 

where n -/ — \° (8) 

'dD/dN BA sic' 

• Subcycles are computed using the rainflow cycle counting method based 
on reversals of stress range ratio 

• The primary creep ductility term automatically reduces the damage done 
by subcycles when their maximum stress is lower than that of the major 

cycle 


THE RATIO CONCEPT PERMITS CDA TO HANDLE 
NON-ISOTHERMAL LOADING 


Many practical applications involve varying temperature 

• The reference stresses and strains are functions of temperature for each 
material in the CDA library 

This is sufficient for variable temperature components below the range 
where TMF is a concern 


• For TMF, an additional strain due lo thermal expansion coefficient 
mismatch was necessary to correlate the data- 

‘ TMF - />' A “(T) f dl (9) 
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THE SECOND DAMAGE MODE PREDICTS 
INTERGRANULAR FAILURE MECHANISMS 


Often components fail by cyclic loads which result in failure by creep mechanisms, 
not fatigue 

• The damage rate for this mode is based on an Arrhenius function, but 
with constants which are different from those used for the transgranular 
mode: 

dJ2 - /A 2 e B2o e-(^j:)dt 

dN (10) 

• This rate is integrated cycle-by-cycle along with the transgranular damage 
rate (and any other active modes) until some mode reaches a value of 
unity 


ENVIRONMENTAL DAMAGE CAN BE INCLUDED FOR 

CYCLIC AND STATIC OXIDATION 

1 "" 1 

Certain components can switch between these modes, depending on how they 

are operated 

• For cyclic oxidation, a multiplier factor based on a time-temperature 
dependent integral is added to the transgranular damage model 

• For static oxidation, a separate third damage mode can be calculated 
which depends only on time and temperature, not cyclic variables such 
as stress or strain 
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INTERACTIONS BETWEEN DAMAGE MODES 
CAN BE EASILY INCORPORATED 


• Where interactions are observed among the various damage modes 
(creep/fatigue/environment), these are handled by combining damage 
levels from those modes and integrating the result as a separate mode 

• Almost any type of interaction may be modeled using this technique: 

1. Linear (additive) 

2. Non-linear 

3. Power law 


COATING LIFE PREDICTIONS UTILIZE THE MODEL 
FROM A COMPANION CONTRACT 


Marry modern superalloys must be coated for oxidation protection in actual sen/ice; 

this must be accounted for by a practical system 

• First, a viscoplastic constitutive model is used to predict the local 
stress/strain history in the coating itself 

• This coating stress/strain history is next used to predict coating cracking 
life (through the coating thickness) 

• Finally, this life is converted to an equivalent damage rate which can be 
integrated along with the other active damage modes for the substrate 
material 
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MULTIAXIAL STRESS/STRAIN STATES ARE CONVERTED 
TO EQUIVALENT CDA VALUES 


Many actual components are analyzed using 3-D methods which produce full 
tensors; these must be interpreted by a practical life system 

• The applied stress/strain tensors are transformed into 100 different axis 
orientations in the plane of the component surface to search for the crack 
plane (with the “worst” orientation) 

• Two different algorithms are currently implemented for choosing the plane 
which experiences the worst damage: 

1. Maximum normal strain range 

2. Socie parameter (combination of shear and normal) 

• New materials may use these or choose some other appropriate 
parameter 


THE EFFECTS OF COMPLEX STRAIN AND 
TEMPERATURE HISTORIES MUST BE DETERMINED FOR 
ACCURATE LIFE PREDICTIONS 


Hot Spot 

created by impingement 
of gas path air 



Outer vane support 
TMF crack location 



strain-temp-time histories 
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CDA USES THE GIVEN STRESS/STRAIN HISTORY 
DIRECTLY TO PREDICT UFE 


• The CDA program accepts as input one or many “defined cycles” for 
which a full stress/strain history is known. These may come from many 
sources: 

1. Actual rig data taken during specimen testing 

2. Simplified methods using elastic input 

a. Neuber's rule 

b. Integrated energy density 

3. Full inelastic analysis methods 

a. Finite element 

b. Boundary integral equation 

• The program performs semi-log interpolation of damage rates to 
determine intermediate values between the defined cycles 


THE CDA PROGRAM CAN BE EASILY EXPANDED FOR 
NEW MATERIALS OR METHODS 


“The only thing constant is change ..." 

• The current CDA equations are sufficient for most current nickel-base 
alloys; a new set of constants can be generated for materials other than 
B1900+Hf or INCO 718 

• The program is modular so that new equations, damage modes, or even 
completely different life prediction methods can be added if desired 

• The program is written in standard FORTRAN-77 and has run on 
machines from PC’s to mainframes; source code will be available for 
distribution soon 
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FUTURE RESEARCH 


Practical methods must undergo continuous improvement and continue to adapt 
to new materials and analytical methods 

• The extensive database for B1900+Hf, including both constitutive and 
creep-fatigue data, provides an excellent opportunity for evaluating new 
prediction methods 

• TMF cycle path effects (such as CW vs. CCW elliptical cycles) are not 
fully explained by the current CDA formulation; possible new methods 
include: 

t. Equilibrium stress/strain methods 
2. Micromechanical models 

• The CDA model could be applied to the constituent materials in 
advanced systems such as metal matrix composites 
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NOTCHED FATIGUE OF SINGLE CRYSTAL PWA 1480 

AT TURBINE ATTACHMENT TEMPERATURES 5 3 ‘ 3 7 


T, G. MEYER, UNITED TECHNOLOGIES RESEARCH CENTER 
D. M. NISSLEY AND G. A. SWANSON, PRATT & WHITNEY 
EAST HARTFORD, CONNECTICUT 


This paper presents results obtained under Option 1 of NASA contract NAS3-23939, 
“Life Prediction and Constitutive Models for Engine Hot Section Anisotropic Material”. 
The Option 1 portion of this contract is focused on the lower temperature, uncoated and 
notched features of gas turbine blades. Constitutive and fatigue life prediction models 
applicable to these regions are being developed. This paper will present fatigue results 
obtained thus far. The program will continue throughout the current year. 

Fatigue tests are being conducted on PWA 1480 single crystal material using smooth 
strain controlled specimens and three different notched specimens. The smooth fatigue 
specimens have a circular cross section and a one inch gage length. The notched speci- 
men configurations are shown in Figure 1. The threaded ends of the thick specimen are 
not shown in Figure 1. The relative orientation of the geometric and crystallographic axes 

are uniquely defined by a specimen’s primary orientation, P , (corresponding to the load- 
ing direction and also to the crystal growth direction of the bars from which the specimen 
is rpachined) and its secondary orientation, S^, (corresponding to the direction normal to 
the notch). Six different notched specimen orientations will be tested in the program. 
Thus far, notched fatigue tests have been conducted on three of these orientations: <001 > 
<I00>, <001 > <210>, and <11 1> <01 !>. 

Isothermal fatigue tests have been conducted at 1200°F, 1400°F and 1600°F. The 
bulk of the tests have been conducted at 1200°F and are presented here. The strain 
controlled tests were conducted at 0.4% per second strain rate and the notched tests were 
cycled at 1.0 cycle per second. Previous tests have shown no material rate dependence at 
1^00°F. A clear orientation dependence is observed in the smooth strain controlled fa- 
tigue results as shown in Figure 3. Figure 4 indicates that much of this apparent orienta- 
tion dependence is reduced when the life data is plotted versus stress range. The fatigue 
lives of the thin, mild notched specimens agree fairly well with this smooth data as shown 
ip Figure 5 when elastic stress range is used as a correlating parameter. Finite element 
analyses were used to calculate notch stresses. The finite element meshes used in the 
analyses are shown in Figure 1. 

In the notched specimens, the locations of maximum principal stress, maximum 
principal strain and maximum octahedral slip system shear stress do not always coincide. 
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f atigue cracks initiate at micropores at the maximum principal stress location and the 
crack plane is normal to the principal stress. Final fracture occurs along a <11 1> crystal 
plane. Figure 6 shows the initiation site and final crack planes for a thin, mild notched 
specimen. In this specimen, the maximum principal strain occurs at the minimum cross 
section while the maximum principal stress is calculated to occur 26 degrees from the 
minimum section point. In the smooth specimens, cracks initiate at pores near the surface 
on planes perpendicular to the applied stress. 

Even though the data shown in Figure 5 appears to be fairly well correlated, a closer 
examination shows that some orientation dependence is still manifested. This can be seen 
clearly in the statistical distribution of the data about the mean trend line as shown in 
Figure 7; all <11 1 > specimens have lives greater than the mean trend line. This indicates 
that the fatigue models being developed cannot depend solely upon stress range. 

Two of the notched specimens being tested can be considered to be thick in that the 
ratio of actual thickness to notch radius is greater than or equal to 1.0. Figure 8 shows 
that these specimens exhibit a somewhat different life trend from the thin specimens. 
These results indicate that correlating parameters should include the effect of the trans- 
verse stress. 

When the material is hot isostatically pressed (HIP) to remove casting micropores, a 
clear improvement in fatigue life is observed as shown in Figure 9. The initiation sites are 
again at the maximum principal stress location with the crack plane normal to the princi- 
pal stress. 

Fatigue testing will continue to further explore the trends observed thus far. Constitu- 
tive and life prediction models are being developed. ’ <• 


• Thin, mild notched specimen 



Notch radius - 0.125 in. 
Thickness = 0.050 in. 


• Thin, sharp notched speciment 



Notch radius = 0.050 in. 
Thickness = 0.050 in. 


• Thick, mild notched specimen s\ 



Thickness = 0.200 in. 


Figure 1. Notched fatigue specimen geometries and finite element meshes used to 
determine notch stresses. 


P s 


<101 > < 101 ^ 



P s 
<iii> <oii> 


Figure 2. Notched specimen orientations that will be tested. 
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1200F PWA 1480 STRAIN CONTROLLED FATIGUE 

2.0 
1.8 
1.6 
1.4 

Strain 

range, % 12 

1.0 
0.8 
0.6 
0.4 

10 10 2 10 3 10 4 10 5 10 6 
Cycles to failure 

Figure 3. A clear orientation dependence is seen in the smooth strain controlled fatigue 
lives when plotted versus strain range. 


1200F PWA 1480 STRAIN CONTROLLED FATIGUE 
400 


300 

Stress 

range, 

ksi 

200 


100 

10 10 2 10 3 10 4 10 5 10 6 
Cycles to failure 


Figure 4. Stress range is a better life correlating parameter. 
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black and white photograph 


1200F PWA 1480 FATIGUE 
Uniaxial and mild notched specimens 



Figure 5. Mild notched and smooth strain controlled lives follow the same trend line. 



| ORIGINAL PAGE IS 

<ioo> 0F P0 °* quality 


Figure 6. Fatigue cracks initiate at the maximum principal stress location with crack 
planes perpendicular to the principal stress. Cracks progress along <11 1> 
planes in the final stages of crack growth. 
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STATISTICAL DISTRIBUTION ABOUT MEAN TREND LINE 


Cumulative 
% failed 



Figure 7. The statistical distribution of failure lives relative to the mean trend line 
shows that all < 1 1 1> specimens have lives greater than the mean trend. 


1200F PWA 1480 FATIGUE 
Effect of thickness 



Figure 8. Thin, sharp notched and thick, mild notch specimens display a life trend 
different than the thin mild notch and smooth specimens. 
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1200F PWA 1480 FATIGUE 


Effect of HIP on thin and smooth specimens 



Figure 9. Fatigue life is impi ved by hot isostatic pressing. 
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TENSILE AND FATIGUE BEHAVIOR OF TUNGSTEN/ 
COPPER COMPOSITES 


-V-2/ 

/$$&/ 


M. J. Verrilli, T.P. Gabb, and Y.S. Kim 
NASA Lewis Research Center 


Work on W/Cu unidirectional composites was initiated to study the behav- 
ior of this ductile-ductile composite system under thermomechanical fatigue 
and to examine the applicability of fatigue-life prediction methods for ther- 
momechanical fatigue of this metal matrix composite. 

The first step has been to characterize the tensile behavior of four ply, 
10 vol.% W/Cu plates at room and elevated temperatures. Since it is envisioned 
that unidirectional composites in service will likely see off-axis stresses, 
the tensile behavior of four fiber orientations are being studied: 0°, 5°, 15°, 
and 90° (where fiber orientation is the angle between the applied stress and 
the fiber direction). The 0° and 90° tests are the upper and lower bounds of 
tensile strength and modulus. The 5° and 15° tests are intermediate cases. 
Strain-controlled experiments at 25 and 260 °C have been performed to date. At 
25 °C the tensile strength, yield strength, and modulus of the 0°, 5°, and 15° 
specimens were about the same. The 90° specimen had lower strengths and modu- 
lus. At 260 °C the strengths and the modulus were largest for the 0° specimen 
and decreased as the orientation angle increased. These properties were con- 
sistent with the rule of mixtures at both temperatures. The yielding behavior 
of the 0°, 5°, and 15° specimens was discontinuous at 25 °C but was smooth at 
260 °C. 


The fracture surfaces of the 0°, 5°, and 15° specimens had the same pre- 
dominant features. Fiber pullout predominated, but, because of a very strong 
fiber-matrix bond, the W fiber surfaces have tightly adherred Cu matrix scales, 
covering about 20 to 25 percent of the fiber surface. Ductile microvoid co- 
alescence predominates in the copper matrix. Longitudinal sections of the 
failed specimens reveal that the fibers each have multiple breaks, up to 25 to 
50 mm from the fracture surface. Upon closer examination, each individual fi- 
ber break is found to be two closely spaced fiber fractures. 

The 90° specimens fail by transverse splitting of the fibers. The matrix 
between the fibers fails by ductile microvoid coalescence. The fibers remain 
very well bonded to the matrix, with no significant debonding even on the frac- 
ture surface. The fibers are all split up to 35 mm away from the fracture sur- 
face. Some of these splits start inside the fibers, away from the machined 
surface, although most of the observed fibers had at least one wide split, 
which extended across the width of the specimen. 

Fatigue experiments were conducted in load control on 0° specimens at 
260 °C. The maximum cyclic stress was varied but the minimum cyclic stress 
was kept constant. Cycle frequency was about 3 cpm, which gave strain rates 
similar to those used in the tensile tests. All tests were performed in vacuum. 
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The cyclic mean strain increased throughout the fatigue tests, suggesting 
that there may be a creep fatigue process controlling the failure. The elas- 
tic modulus remained constant until just before failure. The strain at fail- 
ure increased with increasing maximum cyclic stress. 

Future experimental work will include fatigue at 550 °C and thermomechani- 
cal fatigue between 260 and 550 °C. In addition, the availability of a 12-ply 
composite plate should enable fully reversed, strain-controlled fatigue to be 
performed. The applicability of currently available life prediction techniques 
to the data generated will be examined. 


OUTLINE 

1 . MATERIAL AND TEST PROCEDURES 

2. TENSILE BEHAVIOR 

3. FATIGUE BEHAVIOR 

4. SUMMARY 

5. FUTURE WORK 
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TUNGSTEN-FIBER-REINFORCED COPPER 



500 pm 


9 VOLUME PERCENT UNIDIRECTIONAL W FIBERS, 4 PLIES 

MATRIX: OFHC COPPER, ASTM GRAIN SIZE #8-9 
FIBERS: G.E. 218 W WIRE (0.008 in. DIAM.) 


CO-89-39143 


TEST CONDITIONS 



TENSILE 

FATIGUE 

TEMPERATURE, °C 

25, 260 

260 

FIBER ORIENTATIONS 

0°, 5°, 15°, 90° 

00 

TEST CONTROL 

STRAIN 

LOAD 

TEST RATE 

2.0x1 0* 4 /SEC 

3 CPM 

ENVIRONMENT 

VACUUM 

VACUUM 
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TENSILE DATA FOR 9 VOLUME PERCENT W/Cu 

T = 25 °C 



TENSILE DATA FOR 9 VOLUME PERCENT W/Cu 

T = 260 °C 
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TENSILE FAILURE BEHAVIOR FOR W/Cu AT 25 ®C 

0°, 5°, AND 15® SPECIMENS 


FIBER FRACTURES 



5.6 X 


PAIRS OF FIBER BREAKS 
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TENSILE FAILURE BEHAVIOR FOR W/Cu AT 25 °C 

90® SPECIMEN 


FIBER SPLITTING 



50 X 


INTERGRANULAR FIBER SPLIT 



■ 

20 /xm 
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FATIGUE DATA FOR 9 VOLUME PERCENT W/Cu 

T = 260 ®C 



MAXIMUM CYCLIC STRESS VERSUS LIFE FOR 
FATIGUE OF 9 VOLUME PERCENT W/Cu 

T = 260 °C 
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SUMMARY 


TENSILE BEHAVIOR 

1. AT 25 °C, 00,50,15° SPECIMENS HAD THE SAME TENSILE 
PROPERTIES. THESE VALUES WERE CONSISTENT WITH RULE 
OF MIXTURES. 

2. THESE SPECIMENS HAD DISCONTINUOUS YIELDING CAUSED 
BY PROGRESSIVE FIBER FAILURE. 

3. AT 25 AND 260 °C THE TRANSVERSE (90°) SPECIMEN HAD 
THE LOWEST PROPERTIES. FAILURE AT 25 ®C WAS CAUSED 
BY FIBER SPLITTING. 

FATIGUE BEHAVIOR 

1 . CYCLIC MEAN STRAIN INCREASED CONTINUOUSLY TO FAILURE. 

2. ELASTIC MODULUS REMAINED CONSTANT UNTIL FAILURE. 

3. STRAIN AT FAILURE INCREASED WITH INCREASING MAXIMUM 
CYCLIC STRESS. 


CO-89-39171 


FUTURE WORK 

1. ISOTHERMAL FATIGUE AT 550 »C 

2. THERMOMECHANICAL FATIGUE 

3. FULLY REVERSED FATIGUE 

4. LIFE PREDICTION 
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CUMULATIVE CREEP FATIGUE DAMAGE IN 316 STAINLESS STEEL 

Michael A. McGaw 
NASA Lewis Research Center 


The cumulative creep-fatigue damage behavior of 316 stainless steel at 
1500 °F has been experimentally established for the two-level loading cases of 
fatigue followed by fatigue, creep fatigue followed by fatigue, and fatigue 
followed by creep fatigue. The two-level loadings were conducted such that 
the lpwer life (high strain) cycling was applied first for a controlled number 
of cycles and the higher life (low strain) cycling was conducted as the second 
level, to failure. The target life levels in this study were 100 cycles to 
failure for both the fatigue and creep-fatigue low life loading, 5000 cycles 
to failure for the higher life fatigue loading and 10 000 cycles to failure for 
the higher life creep-fatigue loading. Results of the fatigue followed by fa- 
tigue loading experiments follow the behavior predicted by the double linear 
damage rule and damage curve approach of Manson and Halford (ref. 1), while 
the remainder of the two-level experiments produced results that were inade- 
quately predicted by this approach. These experiments all involved creep- 
fatigue loading and likely introduced a different damage mechanism from a 
cumulative damage standpoint. The failed specimens are being examined both 
f ractographically and metallographically to ascertain the nature of the damag- 
ing mechanisms that produced failure. Models of creep-fatigue damage accumula- 
tion are being evaluated and knowledge of the various damaging mechanisms is 
necessary to ensure that predictive capability is instilled in the final fail- 
ure model. 


EXPERIMENTAL RESULTS AND DISCUSSION 

A strainrange partitioning (SRP) characterization of this material was 
performed for the PP and CP damage modes at 1500 °F. While a complete SRP 
characterization for this material at this temperature has been established by 
Kalluri (ref. 2), the emphasis in that work was on the effect of exposure time 
on the low-cycle fatigue life; the life data so developed are over a range 
that was considered too narrow for the purposes of the present cumulative 
creep-f atigue damage study. Since, in fatigue, the nonlinear cumulative ef- 
fects depend strongly on the relative life levels involved in the loading, the 
present study is focused on obtaining failure data at two widely separated 
life levels, l?oth under pure fatigue and creep-fatigue loading. To this end, 
baseline SRP experiments were performed to establish the inelastic strain 
range levels for 100 and 5000 cycles to failure in PP and 100 and 10000 cycles 
to failure in CP. Results of these baseline characterization experiments are 
shown in figure 1. The resulting life relationships expressed in terms of in- 
elastic strain range versus cyclic life are substantially different in the low- 
er life regime and appear to approach each other in the longer life regime. 
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These life relationships were used to determine the life levels corresponding 
to the loading levels applied in the two-level loading experiments. 

The results of the two-level loading experiments involving PP followed by 
PP loadings are shown in figure 2. Substantial cyclic hardening was observed 
at the low life level, which served to influence the cyclic behavior observed 
at the (subsequent) high life level. The material response at the second load 
ing level was different from that observed for experiments involving that load 
ing level alone, in the respect that it exhibited much higher stresses for the 
applied strains. As the second level cycling progressed, the material cycli- 
cally softened and tended to approach the stabilized cyclic response observed 
in the baseline fatigue characterization experiments. The failure data in fig 
ure 2 are well represented by the damage curve approach of Manson and Halford 
(ref. 1), and are unconservative ly represented by the linear damage rule of 
Miner ( ref . 3 ) . 

The results of the two-level loading experiments involving PP followed by 
CP loadings are shown in figure 3. In this case as well, tl\e influence of the 
substantial amount hardening induced by the first level, lower life cycling 
manifested itself in terms of a somewhat different cyclic response at the sec- 
ond, higher life level, CP loading. In the second level the material exhibit- 
ed much higher cycle times, which, as the cycling went on, tended to approach 
the cycle times observed in the baseline CP behavior. The data in figure 3 
show a trend that is reasonably represented by the linear damage rule. A pre- 
diction based on the damage curve approach is also shown, which, while provid- 
ing a very conservative estimate for the results, does not adeguately repre- 
sent the observed behavior. The presence of creep damage in the second load- 
ing level appears to influence the cumulative damage behavior, implying that 
another damage mechanism is operative. Fractographic and metallographic 
analysis of the fractured specimens is currently underway to, ascertain the 
nature of the damaging mechanisms operative in this type of cumulative fatigue 
loading. 

The results of the two-level loading experiments involving CP followed by 
PP are shown in figure 4. In this series of tests the effect of the CP low 
life level cycling on the second, high life level PP cycling was to introduce 
a significant tensile mean stress as well as to produce a "harder”, cyclic 
stress strain response than was observed in the baseline PP experiments at the 
same applied total strain range. This mean stress tended to relax somewhat as 
the second load level cycling progressed, but a significant mean stress was 
nevertheless maintained throughout the test. The life data shown in figure 4 
were calculated assuming R = -1, and the results of these experiments indi- 
cate that the influence of the earlier CP cycling on the subseguent PP cycling 
is, at worst, virtually insignificant, and, at best, an enhancer of fatigue re- 
sistance at the second level. These results also indicate the possibility of 
another damage mechanism. Fractographic and metallographic analyses are in 
progress . 

The results obtained thus far indicate that for this material in cumula- 
tive creep fatigue loading, the interaction between PP and CP loadings is not 
as deleterious as the interaction that occurs between PP and PP loading. This 
result suggests that another cumulative damage mechanism is operative when CP 
loading is followed by PP loading or vice versa. While this is not surprising 
(the SRP approach is predicated on the possibility of different damaging modes 
of straining), the direction of the influence in the cumulative case, is an 
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unexpected experimental result. Future work will include characterization of 
the PP and PC and the PP and CC cumulative damage behaviors so that the devel- 
opment of a physically based model that describes the general cumulative da- 
mage behavior for this material under conditions of fatigue and creep-fatigue 
loading is achieved. 
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SRP PP AND CP CHARACTERIZATIONS 

316 SS AT 1500 »F 
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Figure 1 
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PP PLUS CP INTERACTION 
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Figure 3 

CP PLUS PP INTERACTION 

316 SS AT 1500 °F 



Figure 4 
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FATIGUE BEHAVIOR OF A SINGLE-CRYSTAL SUPERALLOY 

Sreeramesh Kalluri 
Sverdrup Technology, Inc. 

NASA Lewis Research Center Group 




and 


Michael A. McGaw 
NASA Lewis Research Center 


Single-crystal superalloys have been used extensively during the past dec- 
ade in the manufacture of turbine blades of aircraft engines. At the present 
time, a single crystal superalloy, PWA 1480 is under consideration as a re- 
placement material for the turbine blades of the high pressure fuel turbopump 
(HPFTP) of the space shuttle main engine (SSME). It is necessary to character- 
ize the fatigue behavior of this material for the safe and reliable operation 
of the SSME. In order to achieve this goal, NASA Lewis Research Center pur- 
chased during the early 80 • s cast single crystal PWA 1480. The material was 
cast by TRW, Inc., both in the form of cylindrical bars (11/16 In. diam. and 6 
in, long) and slabs (6.5 by 2.5 by 0.5 in.) along the (00l) orientation. Even 
though each bar and slab was cast individually, all nominally had the same 
chemical composition (table I). The cast single crystal PWA 1480 material was 
observed to contain microporosity, especially between the dendrites (refs. 1 
and 2 ) . 

Three separate experimental programs were conducted to characterize the 
fatigue behavior of this alloy at two laboratories (refs. 1 to 3). Fatigue 
tests were conducted at room temperature (in air) and at 1000 °F (in vacuum) 
on smooth specimens machined from both the cast bars and slabs. In this study 
the data from all of these programs are consolidated to provide a broader char- 
acterization pf the fatigue behavior of the single crystal PWA 1480. The zero- 
mean-stress fatigue life relationships are expressed in terms of stress range 
versus cyclic life lines on log-log plots (figs. 1 and 2). At room temperature 
the slope of the life relationship of this single-crystal superalloy is larger 
than that observed typically in polycrystalline alloys. For' a polycrystalline 
alloy, the slope of this line is usually between -0.07 and -0.14, whereas for 
PWA 1480 single crystal the slope is -0.17. A large amount of scatter was ob- 
served in the room temperature fatigue data generated at NASA Lewis Research 
Center. This scatter is attributable to the porosity observed both within the 
material and on the surface of the fatigue specimens (ref. 1). However, it is 
interesting to note that the data generated from both the cast bars and slabs 
are within the same scatter band (fig. 1). Also, the fact that two different 
laboratories and three separate test systems that were involved in generating 
the data does not appear to have affected the scatter in the data to a signifi- 
cant extent. Only two fatigue data points were available for the zero mean 
stress condition in vacuum at 1000 °F. The life relationship at this tempera- 
ture was determined by using the slope of the room temperature zero-mean-stress 
fatigue data (fig. 2). 


39 


PRECEDING PAGE BLANK NOT FILMED 


Characterization of the fatigue behavior of (OOl) oriented PWA 1480 sin- 
gle crystal under conditions of tensile mean stress was an objective of the 
work done at Argonne National Laboratory (ANL) (ref. 3). Fatigue data with 
tensile mean stresses were also generated by Rockwell International on this 
material (ref. 4). The material tested at ANL had the same chemical composi- 
tion as the NASA Lewis material, while the Rockwell material had a chemical 
composition that was not significantly different from the material used at 
NASA Lewis (table I). In the current study an attempt was made to character- 
ize the effect of tensile mean stress on the fatigue life of PWA 1480 single 
crystal by using the unified approach proposed by Heidmann (ref. 5). In this 
approach the fatigue life is modified by a mean stress parameter so that a sin- 
gle life relationship can be used to represent both zero and tensile mean 
stress data. A better correlation of the fatigue data was obtained by modi- 
fying the functional form proposed originally in reference 5. When tensile 
mean stresses are properly accounted for, the data with and without mean 
stress collapse to a single line (figs. 3 and 4). The implications of this 
type of treatment of mean stress effects, in terms of the conventional ap- 
proach, are shown in figure 5 (room temperature in air) and Figure 6 (1000 °F 
in vacuum). These plots indicate that tensile mean stresses are likely to be 
more detrimental in low-cycle fatigue than in high-cycle fatigue for PWA 1480 
single crystal. However, this type of behavior is exactly the opposite of the 
behavior exhibited by most polycrystalline alloys (ref. 6). The effect of ten- 
sile mean stress on the fatigue life appears to diminish at 1000 °F compared 
with room temperature. This result may be due to the differences in the defor- 
mation mechanisms at these two temperatures and environmental conditions. 
Additional experimental fatigue data on PWA 1480 single-crystal superalloy 
with different levels of mean stress and cyclic lives are necessary to confirm 
the trends of mean stress effects observed in this analysis. 
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CHEMICAL COMPOSITION OF CAST 
PWA 1480 SINGLE CRYSTAL 

[WEIGHT PERCENT] 


ELEMENT 

NASA LEWIS 
BARS AND 
SLABS 
(REF. 1) 

ROCKETDYNE 
MATERIAL 
(REF. 4) 

CHROMIUM 

10.40 

10.16 

COBALT 

5.30 

5.35 

TUNGSTEN 

4.10 

4.13 

TANTALUM 

11.90 

11.95 

ALUMINUM 

4.80 

4.91 

TITANIUM 

1.30 

1.35 

NICKEL 

BALANCE 

BALANCE 


Table I 
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PWA 1480 SINGLE CRYSTAL ZERO-MEAN-STRESS 
(R = -1) FATIGUE DATA 

ROOM TEMPERATURE IN AIR 
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PWA 1480 SINGLE CRYSTAL ZERO-MEAN-STRESS 
(R = -1) FATIGUE DATA 

1000 OF IN VACUUM 



LIFE. N fo 


Figure 2 n 


PWA 1480 SINGLE CRYSTAL 
TENSILE-MEAN-STRESS FATIGUE DATA 

1000 OF, IN VACUUM; ARGONNE 


ZERO-MEAN-STRESS LIFE LINE 
LIFE WITH MEAN STRESS, N, 

LIFE CORRECTED FOR MEAN STRESS, N fo 


STRESS 
RANGE, 10 2 
ksi 



LIFE, Nj AND N fo 

1/b WHERE 

S m TENSILE MEAN STRESS 
S f TRUE FRACTURE STRENGTH 
b FATIGUE STRENGTH EXPONENT 
G (N ( ) = 0.007 (log N,) 2 - 0.583 (log N ( ) + 2.0 c — p-mk 

Figure 3 
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PWA 1480 SINGLE CRYSTAL 
TENSILE-MEAN-STRESS FATIGUE DATA 

ROOM TEMPERATURE, IN AIR; ARGONNE AND ROCKETDYNE 


10 3 
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RANGE, 10* 
ksl 
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O LIFE WITH MEAN STRESS, N, 
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G (N ( ) = 0.142 (log N,) 2 - 0.928 (log N,) + 


WHERE 

S m TENSILE MEAN STRESS 
S f TRUE FRACTURE STRENGTH 
b FATIGUE STRENGTH EXPONENT 

co-tf-nm 


Figure 4 


EFFECT OF TENSILE MEAN STRESS ON 
FATIGUE LIFE OF PWA 1480 SINGLE CRYSTAL 

ROOM TEMPERATURE, IN AIR 


LIFE, 

Nf 



S m TENSILE MEAN STRESS 
S f TRUE FRACTURE STRENGTH 
S g ALTERNATING STRESS AMPLITUDE 

S r FULLY REVERSED STRESS AMPLITUDE WHICH PRODUCES THE SAME 
FATIGUE LIFE AS A COMBINED MEAN AND ALTERNATING STRESS 


Figure 5 
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EFFECT OF TENSILE MEAN STRESS ON 
FATIGUE LIFE OF PWA 1480 SINGLE CRYSTAL 

1000 °F, IN VACUUM 



S m TENSILE MEAN STRESS 

S f TRUE FRACTURE STRENGTH 

S a ALTERNATING STRESS AMPLITUDE 

S r FULLY REVERSED STRESS AMPLITUDE WHICH PRODUCES THE SAME 
FATIGUE LIFE AS A COMBINED MEAN AND ALTERNATING STRESS 

CD-B7-3?2tV 


Figure 6 
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The propagation of a crack in hot path components such as combustor liners can 
occur under complicated thermo-mechanical loading histories which may include 
substantial inelastic deformation. Current methods of predicting crack growth 
behavior utilizes linear elastic fracture mechanics (LEFM) which is not accurate for 
these circumstances. Therefore, it is desirable to determine suitable non-linear 
fracture mechanics parameters. The most likely parameters appear to be 
path- independent (P-I) integrals, several of which can be applied to high 
temperature inelastic crack growth problems. A combined analytical and experimental 
evaluation of these parameters is being performed at elevated temperatures under 
isothermal and thermo-mechanical loading, both with and without thermal gradients 
and hold times . 

Several path- independent (P-I) integrals have been reviewed to determine their 
limitations under thermo -mechanical loadings. The basis for selecting candidate P-I 
integrals was that they were pa th - independent for these complex loading conditions. 
Based on these results, the following integrals were selected for more extensive 
evaluations 

1. Blackburn's J*- Integral and its rate form J* 

A A 

2. Kishimoto's J-Integral and its rate form J 

3. Atluri's AT*- and T*-Integral 

Alloy 718 c rack growth experiments were conducted to assess the ability of the 
selected P-I Integrals to describe the elevated temperature crack growth behavior. 
These tests were performed on single edge notch (SEN) specimens under displacement 
control with multiple extensome ters to monitor the specimen and crack mouth opening 
displacement (CMOD) . The displacements in these tests were sufficiently high to 
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induce bulk cyclic inelastic deformation of the specimen. Under these conditions, 
the LEFM parameter K does not correlate the crack growth data. The experimentally 
measured displacements gradient at the end pf specimen gage length were used as the 
boundary conditions in elastic -plast ic finite element method (FEM) analyses. These 
analyses were performed with a node release approach using CYANIDE, a GEAE FEM code, 
which included a gap element which is capable of efficiently simulating c rack 
closure. Excellent correlation was obtained between the experimentally measured and 
predicted variation of stress and CMOD with crack length and the stress-CMOD loops 
for Alloy 718 tests conducted at 538°C. This confirmed the accuracy of the FEM 
crack growth simulation approach. The experimentally measured crack growth rate 
data correlated well the three selected P-I integrals. The selection of the most 
accurate parameter will be based on analysis of temperature gradient and thermal 
mechanical fatigue crack growth experiments performed using Alloy 718. 

It Is currently planned to model c rack growth under time dependent deformation 
using the superposition of the crack growth from cyclic or time- independent 
deformation and that from static or time-dependent deformation. The latter crack 
growth may be described using rate integrals (J*, and T*) , The software for 
computation of these P-I Integrals from FEM analyses has been developed. The P-I 
integrals in an SEN specimen have been calculated for both uniform load and uniform 
strain under isothermal conditions and for uniform load with a linear temperature 
gradient, The FEM analyses showed the relaxation and redistribution of the normal 
stress which occurs as a result of the time -dependent deformation. The three rate 
P-I integrals (J*, and AT*) and their rate integrals (J*, *j\ and T*) were 
obtained with reasonable accuracy and showed path - independence for the conditions 
evaluated . 

These investigations have produced significant progress in developing P-I 
integrals as non-linear fracture mechanics parameters. The analytical and 
experimental results to date suggests that this methodology has the potential of 
accurately describing elevated temperature crack growth behavior under the combined 
influence of thermal cycling and bulk e last ic - ine las tit deformation states. 


This work was performed under contract NAS3- 23940 with the NASA- Lewis Research 
Center. T.W. Orange is the Program Manager. 
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TIME INDEPENDENT P-I INTEGRALS 


P-LJNIEGRAL 

NON- PROPORTIONAL 
_CYCLIC_LQADIN(L 

TEMPERATURE 

MATERIAL 

INHOMOGENEITY 

ELASTIC-PLASTIC 

STRAINS 

RICE 

NO 

NO 

NO 

YES 

AINSWORTH 

NO 

YES 

NO 

YES 

WILSON & YU 

NO 

YES 

NO 

NO 

GURTIN 

NO 

YES 

NO 

NO 

BLACKBURN 

YES 

YES 

YES 

YES 

KISHIMOTO 

YES 

YES 

YES 

YES 

ATLURI 

YES 

YES 

YES 

YES 



TIME DEPENDENT P-I 

INTEGRALS 


P-I INTEGRAL 

NON-PROPORTIONAL 
CYCLIC LOADING 

CREEP 

DEFORMATION 

TEMPERATURE 

GRADIENT 

MATERIAL 

INHOMOGENEITY 

BLACKBURN 

YES 

YES 

YES 

YES 

KISHIMOTO 

YES 

YES 

YES 

YES 

ATLURI 

YES 

YES 

YES 

YES 

HUTCHINSON 
GOLDBERG 
LANDES 
& BEGLEY 

NO 

STEADY 

STATE 

ONLY 

NO 

NO 
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DC 

Potentail 

Drop 



SCHEMATIC DIAGRAM OF TESTING METHOD 



HYSTERESIS LOOPS FROM 
538C ALLOY 718 TEST WITH 
A STRAIN RANGE OF 1.15% 
AND A MEAN STRESS OF ZERO 
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718 TESTS 
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FINITE ELEMENT MESH 

AND LINEAR DISPLACEMENT BOUNDARY CONDITIONS 
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CMOD RANGE fmm) 
oio ais 0.20 



COMPARISON OF EXPERIMENTALLY MEASURED 
AND PREDICTED VARIATION OF STRESS AND CMOD 
WITH CRACK LENGTH IN 538C ALLOY 718 CRACK GROWTH TESTS 



CMOD (mm) CMOD (mm) 

COMPARISON OF EXPERIMENTALLY MEASURED AND PREDICTED 
STRESS - CMOD HYSTERESIS LOOPS FOR TWO CRACK LENGTHS 
IN 538C ALLOY 718 CRACK GROWTH TESTS 
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ANALYSIS OF DAMAGE IN MMC COMPONENTS USING 
AN INTERNAL STATE VARIABLE MODEL 

V.K. Arya 

University of Toledo 



The need to model the rate dependence and observed creep and plasticity 
interactions exhibited by materials, particularly at elevated temperatures, 
has greatly stimulated the development of numerous viscoplastic models. These 
models, in general,, provide a better description of high-temperature time- 
dependent inelastic behavior of materials. 

Owing to their lightweight and enhanced strength, the metal-matrix com- 
posite (MMC) materials are attracting considerable attention for high-tempera- 
ture applications. As a result of concerted and leading efforts in this direc- 
tion at NASA Lewis, a metal-matrix composite model was developed by Robinson 
and Duffy (1989). (See figs. 1 and 2.) The concept of damage evolution has 
recently been included by Robinson (personal communication) in the above 
model. The evolution of damage is assumed to be governed by a Kachanov-type 
equation* (See fig. 3.) 

The highly nonlinear and mathematically ’’stiff" nature of the constitutive 
equations of viscoplastic models renders analytical solutions of problems, in 
general, Intractable. It, therefore, becomes mandatory to develop suitable fi- 
nite element or other numerical solution technologies to be able to use these 
models in component design. With this objective in mind, the above viscoplas- 
tic damage model was implemented in the finite element code, MARC. Both uniax- 
ial (creep) and multiaxial (an internally pressurized thick-walled cylinder) 
problems were analyzed using this implementat ion . Some preliminary results 
are presented here. These results consider monotonic (constant) loadings. The 
study of damage accumulation under variable (cyclic) loadings is being under- 
taken. The results of this study will be presented later. 


Figure 4 shows the experimental data (Cooper, 1966) utilized to determine 
the values of anisotropic parameters that appear in the model. 


The creep curves including damage for four fiber orientations are depicted 
in figure 5, As expected, the minimum creep occurs when load is applied in a 
direction parallel to the fibers. 

The tangential strains at the inner radius of a thick-walled MMC-cylinder 
for four fiber orientations are plotted in figure 6. The damage is included. 
The cylinder exhibits the maximum creep resistance when the fibers are oriented 
in the circumferential direction, perpendicular to the axis of cylinder. 

Figure 7 shows the time-to-f ai lure for tbe thick-walled cylinder for the 
same fiber orientation angles. As expected, the life of the cylinder can sub- 
stantially be increased by orientating the fibers in the circumferential direc- 
tion, perpendicular to the axis of cylinder. 
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The results, although qualitative, indicate that significant benefits in 
c reep - res i stance and service life can be achieved by using MMC materials as 
structural materials for high-temper ature design. 

The finite element technology developed herein is proposed to be applied 
to aerospace structural components subjected to (cyclic) thermomechanical load- 
ings. The results of these analyses will be reported subseguently. 

REFERENCES 

Robinson, D.N.; and Duffy, S.F. (1989): A Continuum Deformation Theory for 

Metal -Matrix Composites at High Temperature. Journal of Engineering 
Mechanics, ASCE (under publication). 

Cooper, G.A. (1966): Orientation Effects on Fiber-Reinforced Metals. Journal 

of Mechanics and Physics of Solids, vol. 14, pp. 103-111. 


METAL MATRIX COMPOSITE MODEL 
INCLUDING DAMAGE 


FLOW LAW 
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Figure 1 
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METAL MATRIX COMPOSITE MODEL 
INCLUDING DAMAGE 


EVOLUTION LAW 
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POINT OF THE MATERIAL 
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Figure 2 


METAL MATRIX COMPOSITE MODEL 
INCLUDING DAMAGE 


DAMAGE LAW 


\j/ = — Ao^ /ipq 


o = x (N,S) 


IN WHICH DAMAGE VARIABLE, 


D = 1—^ 


WHERE 

N = MAXIMUM TENSILE STRESS TO THE FIBER DIRECTION 
S = MAXIMUM LONGITUDINAL SHEAR STRESS ALONG FIBER DIRECTION 


FORM OF FUNCTION x USED: 


a = X(W}+ (1 -X) S (0 < A < 1 ) 

Figure 3 


CO-89-3930 I 
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YIELD STRESS OF TUNGSTEN/COPPER FOR 
VARIOUS FIBER ORIENTATION ANGLES 

a = 7; EXPERIMENTAL DATA FROM COOPER (1966) 
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Figure 5 
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TANGENTIAL STRAIN AT INNER RADIUS OF 
METAL-MATRIX COMPOSITE CYLINDER 
(INCLUDING DAMAGE) 


INELASTIC 

STRAIN 



TIME, t, hr 

CO-8?- 39304 


Figure 6 


ROBINSON’S METAL-MATRIX COMPOSITE 
VISCOPLASTIC DAMAGE MODEL 

THICK-WALLED CYLINDER 
(800 °F, p = 3.65 ksi) 



CP-89-39303 


Figure 7 


57 



N91-24316 


TI IE EFFECT OF HYDROGEN ON THE DEFORMATION BEHAVIOR OF A SINGLE 
CRYSTAL NICKEL-BASE SUPERALLOY 


W.S. Walsion and A.W. Thompson 
Carnegie Mellon University 
Pittsburgh, PA 15213 




I.M. Bernstein 

Illinois Institute of Technology 
Chicago, IL 60616 


The major goal of this study is to determine how hydrogen affects the deformation 
and fracture behavior of PWA 1 480, Tensile tests, hydrogen-induced crack growth and 
fracture toughness tests are being performed in order to meet this goal. The role of 
hydrogen trapping sites is also being incorporated by varying the levels of porosity and 
eutectic y. This presentation will concentrate mainly on the effect of hydrogen on the 
tensile deformation behavior. 

All tensile samples tested were within 5° of [0011. The samples were given the 
normal heat treatment resulting in a cuboidal y 1 size of approximately 0.5 Jim. Samples 
were gas-phase charged for 15 days at 350°C and 20 kpsi resulting in a uniform hydrogen 
concentration of 300 ppm by weight (17500 appm). The concentration in the samples was 
analyzed by vacuum hot extraction at 900°C. Tensile tests were interrupted at different 
plastic strain levels to observe the development of the dislocation structure. TEM foils 
were cut perpendicular to the tensile axis to allow the deformation of both phases to be 
simultaneously observed as well as parallel to [ 1 1 1 ] to show the superdislocations on their 
slip planes. TEM observations were performed on a JEOL 120 CX operating at 120 
KeV. 

Similar to other nickel base superalloys, hydrogen was detrimental to the room 
temperature tensile properties of PWA 1 480. There was little effect on strength, however 
the material was severely embrittled. Even without hydrogen, the elongation-to-failure was 
only approximately 3%. The tensile fracture surface was made up primarily of ductile 
voids with regions of cleavage fracture These cleavage facets are the eutectic y in the 
microstructure. It was shown by quantitative fractography that hydrogen embrittles the 
eutectic y and causes the crack path to seek out and fracture through the eutectic y. There 
was two to three times the amount of cleavage on the fracture surface of the hydrogen- 
charged samples than on the surface of the uncharged samples. 

The effect of hydrogen can also be seen in the dislocation structure. There is a 
marked tendency for dislocation trapping in the y matrix with and without hydrogen at all 


This work is supported by NASA Lewis Research Center 
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plastic strain levels. Without hydrogen there is a high dislocation density in the y matrix 
leading to strain exhaustion in this region and failure through the matrix. The dislocation 
structure at failure with hydrogen is slightly different. TEM foils cut parallel to [1 1 1] 
showed dislocations wrapping around Y precipitates. This is in contrast to other 
observations on high Y Vf superalloys which report long straight screw dislocations 
cutting through the Y precipitates. 1 [001] foils show that there is a lower dislocation 
density in the y matrix which can be linked to the effects of hydrogen on the fracture 
behavior. By seeking out the eutectic Y arid preferentially failing through this phase, it is 
possible that the maximum accommodating strain in the y matrix is not reached. This 
suggests that removing the eutectic Y through heat treatment would result in less hydrogen 
embrittlement. Preliminary results indicate that this is true and also that eliminating the 
eutectic Y results in a significantly higher ductility in all specimens. 

The primary activity in the Y precipitates is in the form of superlattice intrinsic 
stacking faults (SISFs). These faults have been reported in other ordered alloys and 
superalloys. 2 ’ 3 It is not believed that these SISFs play a large role in the deformation of 
this alloy. In hydrogen-charged samples, "extended Z" configurations have been observed 
in the Y precipitates similar to those in superalloys at high temperatures. 4 This suggests 
that solute hydrogen promotes the cross-slip of superdislocation segments at room 
temperature. This phenomena was observed in previous work at Carnegie Mellon on a 
similar superalloy, 1 and recently room temperature cross-slip has been reported in N^Al. 5 

References 
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2. D.P. Pope and S.S. F.zz, bit. Met. Rev .,29, 136, (1984) 
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60 



Table I. Effect of Hydrogen on <00 1> Room Temperature Tensile Properties. 



a ys (MPa) 

°UTS < MPa > 

e t (%) 

R.A. (%) 

uncharged 

1001 

1120 

2.9 

3.1 

charged 

1014 

1060 

0.38 

2.7 


MPa 




% elong. 


Figure 1. Effect of heat treatment solutionizing time on tensile properties. At 20 hours all of 
the eutectic Y is removed, and at longer times incipient melting becomes a factor. 
This shows that removing the eutectic Y can significantly increase ductility. 


61 
















Figure 2. (a) Overall [001] tensile fracture (b) small ductile voids (c) cleavage facet 

(d) plateau etch showing cleavage on fracture surface and eutectic y in 
microstructure. 
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lack and white photograph 



a 



b 


Figure 3. [001] micrographs of uncharged sample (a) Sp = 1.2% (b) £f = 3.0%. 
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a 



b 

Figure 4. [001] micrographs of charged sample (a) £p = 0.15% (b) £f = 0.4%, weak- 
beam dark field (WBDF). 
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ORIGINAL PAGE 

BLACK AND WHITE PHOTOGRAPH 



Figure 5. WBDF of superlattice intrinsic stacking faults (SISFs) in y precipitate. 




Figure 6. Room temperature cross-slip in hydrogen charged samples. 
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Figure 7. [Ill] micrograph showing superdislocations bowing around Y precipitates. 
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Single crystal nickel base superalloys have been evaluated for potential application in 
the Space Shuttle Main Engine (SSME) high pressure turbopump turbines. Single 
crystals exhibit low cycle and high cycle fatigue life improvements over directionally 
solidified, hafnium modified, MAR-M246 (DS MAR-M246), the current SSME turbine 
blade alloy. The extreme start and stop thermal transients, high rotational speeds and 
high frequency vibrational modes dictate that these properties are life limiting in the 
SSME. In contrast, the long time creep and stress rupture behavior in gas turbine 
engines, for which the turbine blade alloys were developed, are not of concern for the 
SSME application. Although the fatigue life improvements due to a direct material 
substitution are significant, additional gains are possible through the application of 
advanced processing techniques to the single crystal production sequence. The primary 
initiation sites for fracture in single crystal superalloys, especially under cyclic loads, 
are internal casting pores. This casting porosity is inherent due to the dendritic 
solidification of these two phase systems. Pores are caused by shrinkage of the 
interdendritic liquid during final solidification. The primary objective of this program 
was to evaluate the potential Improvements in microstructures and material properties 
due to the reduction in size and density of these casting defects due to the application of 
high thermal gradient casting and hot isostatic pressing. 

PWA 1480 was chosen as the alloy to be evaluated for this program. The alloy had 
previously been chosen for the SSME application after screening of the commercially 
available single crystal alloys. PWA 1480 exhibited superior hydrogen environment 
embrittlement resistance relative to other candidate single crystal superalloys. In 
addition, a significant data base of material properties had been compiled for this alloy. 
Baseline tests for this program were conducted on standard commercial thermal gradient 
cast and standard PWA 1480 heat treated material. High thermal gradient casting was 
evaluated as an avenue for reducing the size of casting porosity. Hot isostatic pressing 
(HIP) was also employed for the elimination of casting pores. An alternate to the 
standard PWA 1480 coating plus diffusion bonding aging heat treatment cycle was also 
evaluated for potential improvements in the properties of interest to the SSME 
application. Microstructural changes associated with the high thermal gradient casting 
process were quantified by measurement of the size and density of the casting porosity, 
amount of retained casting eutectic and dendrite arm spacings. Tensile tests were 
conducted in air at 760°C. Stress rupture tests were conducted in air at 871 °C at an 
initial applied stress of 620 MPa. Statistically determined numbers of low cycle and 
high cycle fatigue tests were conducted, for each material condition, to quantify changes 
in life due to process improvements. Three low cycle fatigue tests were conducted at 
538°C and 2.0% total strain range and eight high cycle fatigue tests were conducted at 
room temperature and at 871 °C, at a stress ratio of 0.47. 


High thermal gradient casting was found to reduce both the size and density of the 
internal casting porosity relative to the standard thermal gradient process. The smaller 
pore size is a result of the decreased dendrite arm spacing afforded by the higher casting 
thermal gradient. Elemental segregation between the dendritic and interdendritic regions 
is also reduced. The reduced pore size provided an increase in both low cycle and high 
cycle fatigue lives relative to the standard gradient material. In addition, the high 
thermal gradient material exhibited a superior combination of tensile strength and 
ductility, though with a reduction in short time stress rupture life. These results are 
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complicated somewhat due to the introduction of the alternate heat treatment schedule 
along with the high gradient casting process. The high cycle fatigue life improvement 
should be more sensitive to the initiating flaw (pore) size, while the low cycle fatigue, 
tensile and stress rupture behaviors are more generally related to the morphology and 
character of the strengthening gamma prime precipitates. 

Test results from the high gradient cast and HIP material were poor due to post-HIP heat 
treatment problems. HIP of the standard thermal gradient material was followed by the 
alternate heat treatment cycle. The combination of strength and ductility was improved 
somewhat, especially at 760°C. Stress rupture life was again found to be reduced 
relative to the standard process material. Low cycle fatigue life was approximately 
doubled, by the application of HIP and the alternate heat treatment relative to the 
baseline material. High cycle fatigue lives were conducted at stress levels different from 
the baseline material so a direct statistical comparison is difficult. An increase in life is 
indicated. 

The results of the program have shown an improvement in material microstructure due 
to high thermal gradient casting. Improved homogeneity of PWA 1480 is advantageous in 
providing an improved solution heat treatment window and, potentially, easier HIP. High 
thermal gradient casting improves fatigue life by reducing casting pore size. The 
alternate heat treatment improves the balance of strength and ductility which appears to 
improve low cycle fatigue life, but with a reduction in short time stress rupture life. 
Based upon the testing from this program, hot isostatic pressing appears to afford 
further improvements in cyclic life, though additional evaluation is suggested. 
Development of the alternate heat treatment is not recommended due to the reduced 
stress rupture capability and the need to develop a new properties data base. High 
thermal gradient casting and HIP are recommended for application to single crystal 
castings. 
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Improving Properties of Rocket Engine Turbine Blades 


Technical Change Property Benefit 

• Substitute single crystal for • Eliminate carbides as fatigue 

columnar grained alloy initiation sites 


• High gradient single crystal 
casting 

• Hot isostatically pressed 
single crystal 

• Alternate heat treatment 


• Crystallographic orientation 
control 


• Improve homogeneity, increase 
fatigue life 

• Eliminate porosity as fatigue 
initiation site 

• Improve hydrogen environment 
embrittlement resistance 


• Improve hydrogen environment 
embrittlement resistance, 
control dynamic response 


• Rocket engine turbine blade 
alloy development 




Rockwell International 

Rockctdf «• Oi.ltmo 


• Hydrogen resistant, fatigue 
capable alloy 

880 - 32*2 
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ADVANCED PROCESSING OF PWA 1480 
PROGRAM OBJECTIVE 


EVALUATE THE INFLUENCE OF HIGH THERMAL GRADIENT CASTING, HOT 
ISOSTATIC PRESSING AND ALTERNATE HEAT TREATMENT ON THE 
MICROSTRUCTURE AND MECHANICAL PROPERTIES OF A SINGLE CRYSTAL 

NICKEL BASE SUPERALLOY 


D 1-89-2-2 
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ADVANCED PROCESSING OF PWA 1480 
PROGRAM LOGIC 



ADVANCED PROCESSING OF PWA 1480 
ALLOY SELECTION CRITERIA 


♦ HYDROGEN ENVIRONMENT EMBRITTLEMENT RESISTANCE 
SUPERIOR TO DS MAR-M246 


• EXISTING DATA BASE 


• PWA 1400 HAD BEEN SELECTED FOR SSME TURBINE BLADE DEVELOPMENT 


D1 -89-2-4 
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Influence of Casting Process on Microstructure of PWA 1480 



ADVANCED PROCESSING OF PWA 1480 
CASTING THERMAL GRADIENT 


CASTING PROCESS AVERAGE DAS CALCULATED GRADIENT 



(MICRONS) 

(CVCM) 

STANDARD 

445 

12 

HIGH 

220 

50 


original page is 

POO* quauS 
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ADVANCED PROCESSING OF PWA 1480 
DENDRITE ARM SPACINGS 



D1 -89-2-7 
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ADVANCED PROCESSING OF PWA 1480 
AVERAGE TENSILE TEST RESULTS 


CASTING 

GRADIENT 

HEAT 

TREATMENT 

TEMPERATURE 

(C) 

YIELD 

STRENGTH 

(MPa) 

ULTIMATE 

STRENGTH 

(MPa) 

REDUCTION 
OF AREA 

(%) 

ELONGATION 

(%) 

STANDARD 

STANDARD 

24 

1024 

1075 

12.5 

11.7 

STANDARD 

H IF/ALT. 

24 

989 

1219 

9.2 

9.8 

HIGH 

ALTERNATE 

24 

1080 

1209 

10.3 

10.3 

HIGH 

HIP/ALT. 

24 

973 

1003 

3.8 

NA 

STANDARD 

STANDARD 

760 

1149 

1273 

7.6 

5.0 

STANDARD 

HIP/ALT. 

760 

1067 

1240 

12.9 

13,2 

HIGH 

ALTERNATE 

760 

1110 

1303 

12.6 

NA 

HIGH 

HIP/ALT. 

760 

972 

1136 

24.8 

12.5 
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ADVANCED PROCESSING OF PWA 1480 
AVERAGE STRESS RUPTURE TEST RESULTS 


CASTING 

GRADIENT 

HEAT 

TREATMENT 

TIME TO 

RUPTURE 

(HOURS) 

ELONGATION 

(%) 

STANDARD 

STANDARD 

14 

NA 

STANDARD 

HIP/ALT. 

4 

22 

HIGH 

ALTERNATE 

9.4 

15 

HIGH 

HIP/ALT. 

2.8 

6.2 
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44-4 
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ADVANCED PROCESSING OF PWA 1480 
STRESS RUPTURE FRACTOGRAPHY 



FRACTURE SURFACE OF NON-HIP, HIGH 
GRADIENT STRESS RUPTURE BAR. 
FRACTURE INITIATES AT INTERNAL CASTING 
POROSITY AT 'A' AND LINKS UP BY DUCTILE 
TEARING AT ’B\ 
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ADVANCED PROCESSING OF PWA 1480 
STRESS RUPTURE FRACTOGRAPHY 


IfSSt SC 3 2 


si82 25KU 


Xlfil 188PS U032 


LONGITUDINAL SECTION THROUGH FAILED STRESS RUPTURE BAR SHOWS 
CRACK INITIATION AT INTERNAL CASTING POROSITY (ARROWS) 


D1 -89-2-11 
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ADVANCED PROCESSING OF PWA 1480 
LOW CYCLE FATIGUE TEST RESULTS 


STANDARD/STANDARD 
STANDARD/HIP/ALT 
HIGH/ALT 
HIGH/HIP/ALT 
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CYCLES TO FAILURE 

538C, 2.0% STRAIN RANGE D1 -89-2-1 2 


HBHI 3 SAMPLES 

3 SAMPLES 

2 SAMPLES 

| 3 SAMPLES, POOR SOLUTION HEAT TREATMENT 
I I I I L 


ADVANCED PROCESSING OF PWA 1480 
LOW CYCLE FATIGUE FRACTOGRAPHY 



CRACK INITIATION IS AT THE SPECIMEN 
SURFACE OR AT NEAR SURFACE DEFECTS 


PROPAGATION IS STAGE I ON {1 1 1 } TYPE 
PLANES 


R = -1 .0, 538C 
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ADVANCED PROCESSING OF PWA 1480 
HIGH CYCLE FATIGUE TEST RESULTS 

(R = 0.47) 


CASTING 

GRADIENT 

HEAT 

TREATMENT 

TEMPERATURE 

(C) 

NUMBER 
OF TESTS 

AVERAGE 
CYCLES TO 
FAILURE 

STANDARD 

DEVIATION 

STANDARD 

STANDARD 

24 

8 

58881 

42528 

STANDARD 

HIP/ALT. 

24 

- 

- 

* 

HIGH 

ALTERNATE 

24 

4 

3.13X10 6 

4.6 X 10 6 

HIGH 

HIP/ALT. 

24 

1 

4.3 X 10 6 

- 

STANDARD 

STANDARD 

871 

8 

74320 

60163 

STANDARD 

HIP/ALT. 

871 

- 

- 

^ « 

HIGH 

ALTERNATE 

871 

5 

1.18 X 106 

594000 

HIGH 

HIP/ALT. 

871 

6 

71246 

72817 


* TESTED AT DIFFERENT STRESS LEVELS 

D 1-89-2- 15 


ADVANCED PROCESSING OF PWA 1480 
STANDARD GRADIENT HCF RESULTS 


H = 0.47 

24C 871 C 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


ORIGINAL PA^p 
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ADVANCED PROCESSING OF PWA 1480 
CONCLUSIONS AND RECOMMENDATIONS 


.HIGH THERMAL GRADIENT CASTING 


* PRODUCES LESS SEGREGATED, FINER MICROSTRUCTURE 

* REDUCES PORE SIZE 

* RECOMMENDED FOR POTENTIAL IMPROVED HEAT TREATMENT AND HIP 


. HOT ISOSTATIC PRESSING 


• REMOVES INTERNAL CASTING POROSITY 

• SIGNIFICANT CYCLIC LIFE IMPROVEMENT 

• RECOMMENDED FOR APPLICATIONS INVOLVING FATIGUE LIMITED LIFE 


• ALTERNATE HEAT TREATMENT 

. SOME TENSILE PROPERTY IMPROVEMENT 

* BENEFIT IS NOT SIGNIFICANT FMOUHII TO WARRANT APPLICATION 


D1 -89-2-4 
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OVERVIEW OF THE INSTRUMENTATION PROGRAM 

William C. Nieberding 
NASA Lewis Research Center 


As you have heard or will hear in the rest of this conference, one of the 
major goals of the durability program is to develop computational models of 
the structural behavior, the fatigue life, and the environmental conditions 
pertinent to advanced reusable space propulsion systems. In many cases, how- 
ever, there is very little experimental data available either to verify these 
models or to use as input conditions to the models. 

The instrumentation portion of this durability program is thus aimed pri- 
marily at developing sensors and measurement systems capable of obtaining the 
data necessary for the verification of computational models. Indeed, some of 
these sensor systems may and probably will find their way into the test bed 
engine and possibly even into flight engines. This is increasingly likely be- 
cause of the recent rise in the importance of instrumentation for health moni- 
toring purposes. The major effort, though, is to develop what is needed for 
bench mark experiments in research rigs to more fully define the actual condi- 
tions experienced in such an engine. In some cases the only research rig 
available will turn out to be the test bed engine simply because the condi- 
tions cannot be adequately simulated in less extensive facilities. 

One of the characteristics of measurement systems needed to verify codes 
is that the sensors must be nonintrusive or at least minimally intrusive so as 
not to significantly perturb the conditions being measured. The reason for 
this is that the computer codes do not take into account the perturbations 
that would be caused by an intrusive probe. This leads to a heavy emphasis on 
laser optical techniques and on thin-film sensors. 

Another characteristic of such instruments is that they must be highly ac- 
curate and produce very high spatial and temporal resolution of the parameter 
being measured. Verifying codes demands this because the codes usually calcu- 
late highly resolved maps of such things as velocity, strain, temperature, etc. 

The measurement systems needed generally fall into a number of broad cate- 
gories. First there are the measurements needed on the surfaces of components 
such as turbine blades and vanes. Some of the desired parameters are tempera- 
ture, strain, and heat flux. In this session the significant progress made in 
these measurement areas since the last durability conference will be discussed. 

Another broad category encompasses those measurements needed in the flow 
environment around these components. Here, the desired results are high reso- 
lution maps of such parameters as flow velocity, temperature, density, pres- 
sure, and species concentration. One of the papers in this session deals with 
this type of measurement. 

The remaining category deals with measurements necessary for monitoring 
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the health of the engine. This category, not reported at the last conference, 
has loomed ever more important since the Challenger disaster. One paper deals 
with an optical method for determining the characteristics of the plume. An- 
other looks at holographic measurement of structural damage. 

In summary, the instrumentation program is aimed at high accuracy, high 
resolution, minimally intrusive measurements on space propulsion system compo- 
nents and on the environment that surrounds these components • The primary 
goal is code verification, but some of the results may apply to later flight 
instruments. A more recently added goal is that of health monitoring. 
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TWO-DIMENSIONAL HIGH TEMPERATURE STRAIN MEASUREMENT SYSTEM 

Christian T. Lant 
Sverdrup Technology, Inc. 

Lewis Research Center Group 

and 

John P. Bar ranger 
NASA Lewis Research Center 


The Instrumentation and Control Technology Division of the Lewis Research 
Center has developed the in-house capability to make two-dimensional optical 
strain measurements on high temperature test specimens. The first phase of 
this effort demonstrated one-dimensional strain measurement at temperatures to 
450 °C, with a resolution of 18 microstrain (ref. 1). The phase II effort ex- 
pands the phase I system to provide a two-dimensional strain measurement capa- 
bility. This two-dimensional capability is implemented through a rotatable 
sensitive strain axis. Three components of surface strain can be measured 
automatically, from which the first and second principal strains are calcula- 
ted. The phase II system demonstrated one- and two-dimensional strain measure- 
ments at temperatures beyond 750 °C with a resolution of 15 microstrain. 

The system is based on the one-dimensional speckle shift technique of 
I. Yamaguchi (ref. 2). Some of the features of the technique include noncon- 
tact measurements, automatic cancellation of rigid body motion, no surface 
preparation requirement, and near-real time results. The Lewis system also 
features a short gauge length (<1 mm) and a programmably rotatable sensitive 
axis. 


The speckle shift technique makes use of the linear relationship between 
surface strain and the differential shift of laser speckle patterns in the dif- 
fraction plane. Laser speckle is a phase effect that occurs when spatially co- 
herent light interacts with an optically rough surface. Since speckle is gen- 
erated by any diffusely reflecting surface, no specimen preparation is needed 
to obtain a good signal. Speckle shift is measured from a laser beam incident 
on a test specimen at 30° from the surface normal. By also measuring the 
speckle shift from a beam incident at an equal and opposite angle, potential 
errors due to rigid body motions are eliminated. 

Figure 1 shows a schematic of the phase II optical system. An argon-ion 
laser beam is directed onto a test specimen at an angle. The Pockels cell and 
beamsplitting cube combination switches the beam from -30° to +30° by rotating 
the polarization of the beam. The acousto-optic modulator controls the expo- 
sure time of the line scan camera. The rotating assembly allows the plane of 
incidence of the beams onto the specimen to rotate, changing the sensitive 
axis of the instrument. A pair of 1/4-wave retardation plates maintains the 
polarization state of the beam going into the rotating assembly. Waist posi- 
tioning lenses locate the laser beam waist at the specimen surface, improving 
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the error cancellation of the system. It is suspected that a planar wave 
front at the specimen surface also stabilizes the speckle patterns in the pres- 
ence of nonuniform heating. The optical pallet can be easily mounted on a 
multiaxis positioning table. 

Figure 2 is a view of the actual test setup. The rotating assembly is 
shown mounted on a goniometer. The optical axis is coincident with the goniom- 
eter’s axis of rotation. The test specimen is located within a thermal enclo- 
sure in order to decrease the temperature gradients near the specimen surface. 

Testing was done on a flat specimen of Inconel 600 mounted in a fatigue 
testing machine. A load cell measured the stress on the specimen before and 
after acquiring the speckle data. Figure 3 shows one-dimensional strain data 
at room temperature. The coincidence of the data loading up and loading down 
demonstrates the repeatability of the instrument. Figure 4 shows a two- 
dimensional run at room temperature. Strain components were measured at 0° 
(parallel to the load axis) and at ±45° , and the plots indicate the calculated 
values of the first and second principal strains. The measured values of 
Young’s modulus and Poisson's ratio are in good agreement with handbook values. 

An RF induction heater was used to heat the specimen in the high-tempera- 
ture runs. Good linearity of the principal strain moduli at high temperatures 
indicate precision and stability of the system. However, a systematic error in 
the high-temperature test setup introduced a scale factor in the slopes of the 
two-dimensional stress-strain curves. The scale factor is due to a gradual 
load relaxation during the retrieval of the three components of strain. The 
load relaxation occurred because, due to the spacing of the induction heating 
coils, the ends of the specimen were hotter than the test section; although the 
test section behaved linearly, the specimen ends were under plastic deforma- 
tion. This result underscored the need for the load to be static during the 
system’s data acquisition time of about 10 sec per two-dimensional strain 
point. No high temperature effects, however, have been observed to degrade 
speckle correlation. 

Future efforts will include reading two-dimensional speckle patterns with 
an area array camera. This will allow the correlation peak to be maximized 
and will extend the range of tolerable rigid body motion. Dedicated correla- 
tion hardware will reduce the correlation times to the point that true real 
time measurements can be made. 


REFERENCES 

1. Lant, C.T.; and Qaqish, W. : Optical Strain Measurement System Development - 

Phase I. (Sverdrup Technology, Inc.; NASA Contract NAS3-24105) NASA 
CR-179619 , 1987. 
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PHASE II OPTICAL DESIGN 
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Figure 1 


PHASE II TEST SETUP 
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Figure 2 
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OPTICAL INSPECTION OF SPACE-PROPULSION COMPONENTS 
USING AN INJECTION SEEDED Nd:YAG LASER SYSTEM 


Arthur J. Decker 

National Aeronautics and Space Administration 
Lewis Research Center 


5 / 3 - 3 ^ 


Cleveland, Ohio 44135 


A dual-beam, injection-seeded, Nd:YAC laser has been demonstrated for detecting 

structural defects. This demonstration was part of an ongoing project^ to use 
dual— reference-beam holographic interferometry to inspect space propulsion 
components for cracks, defects, structural failures, structural changes, and gas leaks. 
Potential subjects for inspection include welds, duct work, casings, turbopumps, 
blades, composites, and ceramics. The maximum dimension of an inspected area 
ranges from a few centimeters to a meter. The entire optical inspection system is 
now ready for a demonstration application 

Dual-reference-beam holographic interferometry overcomes the most significant 
limitation of holographic interferometry: interpretation of the recorded fringe 

patterns. Measurement and analysis of the fringe pattern are fully computerized and 
automated. There is also a factor of 10 to 100 increase in sensitivity and dynamic 
range over other holographic methods. The electronic detection and computerized 
analysis of an interference fringe pattern can be accomplished by optical 
heterodyning, quasiheterodyning, or phase shifting. Optical heterodyning was 
selected as the most sensitive. 

Pulsed lasers are needed for inspections in the field between tests and during tests of 
components or engines. The short pulse of the laser ( 7 nanoseconds) makes the 
hologram recording process insensitive to vibration ana motion of the subject. A 
double exposure iccords the vibration or load-induced deflection of the subject. 
The deflection is encoded in the interference fringe pattern reconstructed from the 
hologram. The term "dual reference beam" means that geometrically distinct 
reference beams are used for each of the two exposures of the double-exposure 
hologram. 

Both pulsed-laser and continuous-wave-laser holograms are interpreted by the same 
computerized, automated, readout facility. The recorded scene is reconstructed in 
that facility with duplicates of the two reference beams originally used to record the 
hologram. In optical heterodyning, an acoustooptic frequency shifter shifts the 
frequency of one duplicate relative to the other. A 125 kHz frequency shift is used in 
our facility. A lens forms an image of the fringe pattern and subject. The intensity 
of the fringe pattern varies sinusoidally at 125 kHz, and the phase of this signal 
contains the interference phase A<p . The interference phase, in turn, depends 
upon the structural-deflection fie! i ! £(r) occurring between the two exposures. The 
signal is detected at each point of the fringe pattern by a photodetector positioned 
by a computer controlled, precision, XYZ, positioning system. The signal is sent 
to a phase sensitive detector which measures the interference phase A <p relative to 
a reference. 
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We note that the interference phase A^ in holographic interferometry is most 
sensitive to deflections along the optical axis or viewing direction. Therefore, we 
orient the surface of the subject as nearly perpendicular to the viewing direction as 
possible for an optical inspection. We also evaluate numerically one or more of the 
second derivatives of interference phase 

(P\ A d^Atz> cP A(£ 
dx dy dxdy 

These quantities depend primarily on the deflection-induced change of curvature 
(bending) of the surface or on the spatial rate of twist (torsion) of the surface. 


The bending distribution is sensitive to structural defects such as cracks. The use of 
heterodyne holographic interferometry to monitor the formation of a fatigue-induced 

crack in a blade of composite material has been reported. * We used the same kind 
of fatigue failure test to evaluate the performance of the pulsed laser system for 
recording dual-reference-beam holograms for heterodyne readout. 

The pulsed laser system, intended for field inspections of propulsion components, 
contains two injection seeded, frequency doubled, Nd:YAG lasers. There is a 
l-to-1000 microsecond adjustable delay between laser firings, so that 
double-exposure holograms with different time separations can be recorded. A 
single injection seeder assures that the Nd:YAG lasers operate in a single 
longitudinal mode. The laser system will produce two pulses of green light (532 nm) 
where each pulse has an energy equal to 500 mJ. The spectroscopic line-width of 

the pulses is 0.005 cm * . Consequently, good holograms can be recorded of subjects 
having a projected size up to about a meter and a surface relief (scene depth) of 
about half a meter. 

The pulsed laser was compared with a continuous wave laser (argon-ion) for 
detecting fatigue cracking in a blade of composite material 8 cm long and 2 cm 
wide. The blade was vibrated at a large amplitude in a shaker to induce cracking. 
Periodically, this destruction was interrupted and optical inspections performed. 
The pulsed-laser optical inspections were performed by recording double-exposure 
holograms of the blade as it vibrated in its first, bending mode at a small amplitude. 
The laser firing was synchronized with an accelerometer attached to the blade. The 
two laser beams were each divided, routed into two separate reference-beam 
channels, andcombined carefully in a single object-beam channel. The argon-ion 
inspections were performed by accurately displacing the blade tip between exposures 
of the double-exposure hologram. 

The Nd:YAG laser and argon-ion laser performed equally well for detecting cracking 
in the blade using the procedures above. Both lasers were used with a small angle 
between reference beams (about 0.015 degree). The small angle makes it easy to 
form detectable fringes from holograms recorded in the field. 

Optical inspections of the bending induced by a carefully excited vibrational mode, 
with the laser accurately synchronized with the phase of an accelerometer output, 
have the advantage of being repeatable. Hence changes in complex structures should 
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bo easier to recognize. This method is suitable between runs of a component or for 
acceptance testing. It may be desirable to inspect a component during a run or 
other test. The Nd:YAG laser is able to record up to 10 double-exposure holograms 
per second of such a time varying component. Using film to record the holograms, 
one could conceive of performing 100 optical inspections of a time varying and 
changing component from holograms recorded during a 10 second run. 

The conclusion is that the combination of automated, computerized fringe 
measurement system and injection seeded Nd;YAG laser is both suitable and ready 
for a field demonstration of an optical inspection of a propulsion component based 
on dual— reference-beam holographic interferometry. The detailed inspection 
procedure will depend on the component and conditions. 


1 Decker, A.J.; Krasowski, M.J.; and Kn igulec, M.A.: Optical Inspection of 

Propulsion System Components Using Heterodyne Holographic Interferometry, in 
Advanced Earth-To-Orbit Propulsion Technology 1988, vol. n, NASA Conference 
Publication 3012. 
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TALK AND PROGRAM OBJECTIVES 


DISCUSS PROGRESS IN USING NEW Nd:YAG LASER 
FOR OPTICAL INSPECTION OF SPACE PROPULSION 
COMPONENTS 

INSPECT FOR CRACKS, DEFECTS, STRUCTURAL FAIL- 
URES, STRUCTURAL CHANGES, GAS LEAKS 

INSPECT WELDS, DUCT WORK, CASINGS, TURBOPUMPS, 
BLADES, COMPOSITES, CERAMICS 


INSPECT STRUCTURES WITH MAXIMUM PROJECTED 
DIMENSIONS PER INSPECTION FROM A FEW CENTI- 
METERS TO A METER 


CD-89-39394 


METHOD OF INSPECTION 

DUAL-REFERENCE-BEAM HOLOGRAPHIC INTERFER- 
OMETRY 

AUTOMATED FRINGE MEASUREMENTS USING OPTICAL 
HETERODYNING 

CW LASER IMPLEMENTATION DISCUSSED AT EARTH- 
TO-ORBIT PROPULSION TECHNOLOGY CONFERENCE, 
MARSHALL, MAY, 1988 

PULSED LASERS NEEDED FOR INSPECTIONS IN FIELD 
BETWEEN TESTS AND DURING TESTS OF COMPONENTS 
OR ENGINES 


COSINE INTERFERENCE-FRINGE PATTERN 
BETWEEN BENT AND UNBENT STATES 
OF CANTILEVER 



CD-89-393 96 


PHASE CHANGE FOR DEFLECTION 
OF STRUCTURE 

A<t>(r) * £i — K - r . ) . 6(f) 

X 2 

CD-89-39397 
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Nd:YAG LASER LAYOUT 



CD-89-39398 


HETERODYNE HOLOGRAPHY 

I ■ Io + Ij Cos [A<P ± <wt] 

CD-89-39399 
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STRUCTURAL PROPERTIES FROM 
HOLOGRAPHY 


STRAIN 


^ 2 \dy dx J 



BENDING 

A K n * cos^a + 2 cos a sin a + ^ sin 2 a 
0 dx 2 dxd V dy 2 


TORSION 
d 2 &. 


dfi 

dS 


dxdy 


z cos 2 a + | 
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OPTICAL INSPECTIONS OF COMPOSITE BLADE 


2.5 
2.0 

1.5 

0 


(a) COMPOSITE, VIBRATED 0 SEC. (t» COMPOSITE, VIBRATED GO SEC AT 50 Hz. 



ce) COMPOSITE. VIBRATED 2400 SEC AT 50 Hz. 
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TESTS IN THE FIELD 

NEED DOUBLE-PULSE, PULSED-LASER SYSTEM 
(INJECTION-SEEDED Nd:YAG) 

CAN USE EXISTING AUTOMATED FRINGE MEASUREMENT 
SYSTEM USED FOR CW LASER APPLICATIONS 

NEED WAY TO CREATE INTEREXPOSURE DISPLACEMENT 
FIELD (VIBRATION) 

NEED PROCEDURES FOR INSURING REPEATABILITY OF 
TEST CONDITIONS 


CD-89- 3940 3 
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ORIGINAL PAGE 

black and white photograph 



ORIGINAL PAGE 

BLACK AND WHITE PHOlOGRAPH 

INJECTION SEEDED Nd:YAG LASER SYSTEM 

• TWO INDEPENDENT, Q-SWITCHED, Nd:YAG, 
OSCILLATOR, AMPLIFIER, FREQUENCY DOUBLER 
COMBINATIONS (QUANTEL YG 580 SERIES LASERS) 

• ONE INJECTION SEEDER 

• 500 mJ PER LASER OF 532 nm LIGHT 

• 0.005 cm" 1 LINE-WIDTH (2 m COHERENCE LENGTH) 

• 1 TO 1000 MICROSECONDS BETWEEN PULSES FOR 
DOUBLE EXPOSURE 

• 10 PULSE PAIRS PER SECOND 

CD-09- 39404 


SAMPLE LAYOUT 

DUEL REFERENCE-BEAM HOLOGRAPHIC INTERFEROMETRY 
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ORIGINAL PAGE IS 
OF POOR QUALITY 



FRINGE PATTERN 

SIZE ABOUT 1 m 
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INTEREXPOSURE DISPLACEMENT 

SYNCHRONIZE LASER TO VIBRATIONAL MODE 

COMPONENT OR ENGINE INSPECTION BETWEEN RUNS: 
EXCITE VIBRATIONAL MODE NONDESTRUCTIVELY 

CD-89-39407 
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ORIGINAL PAGE 

BLACK AND WHITE PHOTOGRAPH 

TWO LASER PULSES SYNCHRONIZED WITH 
ACCELEROMETER OUTPUT 



CD-89-37408 


COMPOSITE AT 21 000 FATIGUE CYCLES 


CW LASER PULSED LASER 



CD-89-39409 
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STATUS 


INJECTION SEEDED Nd.YAC LASERS ARE SUITABLE 
FOR. DUAL-REFERENCE-BEAM HOLOGRAPHIC 
INTERFEROMETRY 

SENSITIVITY OF CW AND PULSED LASER INSPECTIONS 
SIMILAR FOR THE SAME REFERENCE BEAM 
INTERANGLE 


3 STANDARD DEVIATIONS: 

1/300 FRINGE, CW, LARGE ANGLE 
1/100—1/60 FRINGE, CW AND PULSED LASER, 
SMALL ANGLE 

CW AND PULSED LASER METHODS CAN NOW BE 
DELIVERED TO SPECIFIC PROJECTS 
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CALIBRATOR TESTS OF HEAT FLUX GAUGES MOUNTED IN SSME BLADES 

Curt H. Liebert 
NASA Lewis Research Center 


Measurements of heat flux to space shuttle main engine (SSME) turbine blade 
surfaces are being made in the Lewis heat flux calibration facility. This fa- 
cility was described at the 1987 Structural Integrity and Durability of Reusable 
Space Propulsion Systems Conference (ref 1). Surface heat flux information is 
obtained from transient temperature measurements taken at points within the 
gauge (refs. 2 to 4). A 100-kW Vortek arc lamp (ref. 5) is used as a source of 
thermal radiant energy. Figure 1 shows a gauge that has been placed at the ex- 
ternal focus of an elliptical reflector in a lamp head attached to the side of 
a service module in the lamp system. 

The schematic drawing of figure 2 details the miniature plug-type heat flux 
gauges being developed for measurement and study of blade surface heat flux. An 
annulus is electrical discharge machined (trepanned) into MAR M 246 (Hf)(DS) 
blade material creating a post, or thermoplug. The annulus is machined only 
part way through the blade thickness. This fabrication method is advantageous 
because a joining process (welding) is not required between the thermoplug and 
the gauge. Therefore, no seams are formed between the thermoplug and blade sur- 
face; seams can disturb boundary layer flow and heat transfer on the gauge sur- 
face. Figure 3 shows the rear of a gauge with the cover removed. Chromel- 
Alumel thermocouple wires (not shown) with diameters of 0.00254 cm are attached 
to the front surface of the gauge and along the sides and on the rear of the 
thermoplug. These wires are routed through the annulus to the back of the blade 
where they are attached to electrically insulated extension wires. The exten- 
sion wires are mounted in grooves machined into the airfoil pressure surface 
(fig. 3), then routed through holes machined into the platform and finally fas- 
tened to the blade shank. 

Thermoplugs, with diameters of about 0.190 cm and lengths varying from 
about 0.190 to 0.320 cm, are being investigated. The thermoplug is surrounded 
on all surfaces except the active surface by a pocket of air located in the cir- 
cular annulus and under the back cover. Since the thermoplug is insulated, it 
is assumed that heat is conducted in a one-dimensional manner from the hot ac- 
tive surface to the cooler back side of the thermoplug. A one-dimensional rela- 
tionship for obtaining surface heat flux is then (ref 2): 


f (pC 8T/8t)dz, W/m^ 

J o p 


( 1 ) 


where q s is the surface heat flux, L is the thermoplug length, and p and 
Cp are the density and specific heat of MAR-M-246(Hf ) (DS) , 8T/8t is the 
partial derivative of temperature with respect to time, and z is a distance 
along the thermoplug axis. 


Values of specific heat as a function of temperature for MAR-M-246 (Hf ) ( DS ) 
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and a similar alloy (MAR-M-200) are shown in figure 4, These data were obtained 
from reference 6. The density of MAR-M-246 (Hf ) (DS) is 8442 kg/m 3 (ref. 6). 
Unfortunately, no indication of experimental error associated with property 
measurement is given in reference 6. Experimental errors of specific heat and 
density data are assumed as ±20 percent and ±5 percent. These uncertainties of 
measured property data are typical for many alloys. The curve of specific heat 
versus temperature used herein to solve equation (1) is shown as a solid curve 
in figure 4. The solid curve and literature data correspond within experimen- 
tal error. 

The suction surface of an instrumented SSME blade mounted on the calibra- 
tor is shown in figure 5, A heat flux gauge is mounted on both the suction and 
pressure surfaces at midspan and midchord. In figure 5 radiation is incident 
on the active surface of the suction surface gauge, and radiation is also inci- 
dent on the back of the pressure surface gauge. Both this and the reverse ori- 
entation are being investigated. Welds attaching a cover plate (304 stainless 
steel) to the suction surface of the blade and a cover plate welded over exten- 
sion wires buried in grooves machined into the airfoil are discernible in fig- 
ure 5. Because there is no airflow over the airfoil surface tested in the cali- 
brator, these welds were not machined flush to the surface. 

The heat flux gauge mounted on the suction surface below the cover is not 
discernible in figure 5 because the thermoplug and blade at this location are 
one piece. Thermoplug temperature history measured within the suction surface 
gauge is plotted in figure 6, and the corresponding heat flux history is shown 
in figure 7. Blade heat flux data are measured over a temperature range of 
about 100 to 1100 K at a lamp current of 400 A. This data are compared in fig- 
ure 7 with data previously obtained with other gauges: (1) four miniature plug- 

type gauges mounted in MAR-M-246 (Hf) (DS ) coupons, (2) two commercial factory- 
calibrated Hy-Cal gauges, and (3) a photocell. The photocell millivolt data 
show that the lamp intensity rises rapidly during startup (time = 0.3 to 
0.6 sec) and then becomes relatively constant at longer times. This trend was 
also observed when the two Hy-Cal gauges and the four miniature plug-type gauges 
mounted in the coupons were investigated in the calibrator. Furthermore, the 
plug-in-coupon data compare reasonably well (-6%) with the Hy-Cal gauge data. 
Surface heat flux measured on the blade instrumented with a miniature plug-type 
gauge compares well with both the Hy-Cal and gauge-in-coupon data (fig 7), 

Because these comparisons are reasonable, it is concluded that the minia- 
ture plug-type gauge concept is feasible for measurement of blade surface heat 
flux. Reference 7 suggests that it is important to measure heat flux near the 
hub on the suction surface and at the throat on SSME blades rotating in engines 
because stress and heat transfer coefficients are high in this region. Nonro- 
tating tests in the calibrator on gauges mounted at this throat location are 
being planned. Also, tests on gauges with lower thermal conductivity insula- 
tion (compared with air) around the thermoplug are planned in the calibrator. 

Use of insulation with lower thermal conductivity will minimize heat conduction 
between the plug, the surrounding annulus, and the cover plate. Furthermore, 
designs for locating heat flux gauges at the leading edge of SSME blades have 
been started. Finally, plans to test heat flux gauges mounted on nonrotating 
SSME blades inserted into the Marshall Space Flight Center turbine blade tester 
are being formulated. The results, in terms of measured gas-side heat transfer 
coefficients, will be compared with calculations (ref. 7). 
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HEAT FLUX MEASUREMENTS WITH GAUGE 
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summary 

During the past twenty years, laser Doppler anemometry has become an important 
non-intrusive diagnostic tool. It has been used for detailed and accurate measurement of 
fluid flows in wind tunnels, air-breathing engines, rocket motors, shock tubes, etc. A recent 
study by Rocketdyne for NASA identified laser anemometry, using a compact optical head, 
as a feasible diagnostic instrument for the SSME Model Verification experiments. PRi is 
presently under contract from NASA Lewis to develop and deliver such a laser anemometer 
system. For this application, it is desired to place the laser at a remote distance from the 
engine, and use single mode polarization preserving fiber optics for the transmission of the 
laser light to and from the measurement head. Other requirements for the instrument 
include: 

1. Two simultaneous components 

2. Flow velocity: -300 to 300 m/s 

3. Flow angle: 0 to 360" 

4. Severe vibration and temperature environments as specified 

5. Measurement resolution: 1 mm (normal to the flow) 

6. Optical probe dimensions: 25 mm diameter; 15 cm length 

7. Total run time: 500 seconds 

8. Accuracy: 5% velocity magnitude, 3 degree in angle 
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Analytical and experimental tools are being used to develop the technologies 
required for the laser anemometer. T hese include finite element analysis of the optical 
head and vibration tests for various optical and mechanical components. Design of the 
optical head and the fiber optic connectors are driven by the temperature and vibration 
requirements for the measurement environment. Results of the finite element analysis and 
the vibration tests of the components are included. Conceptual design of the fiber optic 
launcher and the optical probe has also bee complete. Detailed design of the probe as 
well as the fabrication and assembly of the components is in progress. 

It is anticipated that the anemometer will be fabricated and its performance 
evaluated through a series of vibration and environmental tests before the delivery of the 
optical system to NASA for further evaluation. 
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Finite Element Model of the Probe, Including the lenses and Fiber Holders 
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ABSTRACT 


The purpose of this analysis is to determine the effect of transient heat flux on heat 
flux sensor response and calibration. 

We have begun this work by studying a one dimensional case. This was done to 
elucidate the key parameters and trends for the problem. It has the added advantage that 
the solutions to the governing equations can be obtained by analytical means. The. 
analytical results obtained to date indicate that the transient response of a heat flux sensor 
depends on the thermal boundary conditions, the geometry and the thermal properties of 
the sensor. In particular if was shown that if the thermal diffusivity of the sensor is small, 
sav for a material such as Nichrome, then the transient behavior must be taken info 
account. 
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INTRODUCTION 


The purpose of this analysis was to determine the effect of transient heat flux on 
sensor response and calibration. 

Consider the case of a heat flux sensor mouifted in a wall subjected to both radiative 
and convective heat transfer on the top side and purely convective heat transfer on the 
bottom side as shown in Fig. I, where x and r are the axial and radial coordinates, 
respectively; L and f, D and d are the thicknesses and diameters of the two layers; k is the 
thermal conductivity, a is the diffusivity, T is temperature, Ip and h 2 are the convective 
heat transfer coefficients, q is the radiant incident heat flux, and T is the ambient 
temperature. The subscripts i and 2 refer to the different materials of the sensor. 


In this analysis we make the following simplifying assumptions 

1 11 neglect re— radiation; 

2) constant thermal properties; 

31 constant, uniform ambient temperatures; 

41 constant, uniform convection coefficients; 

51 constant, uniform radiant incident, heat flux; 

6) perfect thermal cont act between t he regions of different k; 

(7) the equations arc one-dijnensional; 

Steady - State Solution To determine a base line case we first, solve for the 
steady— state case. 


For this case, t he governing equations and corresponding boundary condit ions are, 
under the above assumptions. 


d 2 T, _ 
d^ 2 "-° 

d 2 T 2 _ n 

HP "- 0 


0 < x > C 
( < x > L 


(1) 

( 2 ) 


-k 


t 


dT, 

dx 


+ l ',(T|-V-q 0 =0 


at x = 0 


(V 
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(3) and (5) are rearranged as: 



T hi l'i ~ li|l| 
+ T, = h 2 f 2 


a! x — 0 
at x = L 


where f| 




f . = To 

~00 


(3-a) 
(5 — a) 


The solutions to equs. ( I ) and (2) arc 

T] = Ai x + B| 

To = Ao x + Bo 

The coefficients A|. B|, Ao, I3o can he determined in terms of the Iroiindaiy conditions 
(3— a). (1), and (- r )-a). 

The sensor output is based on the temperature difference between two fixed points 
where the two alumel wires are attached. Assume t hat the coordinates of these two points 
are xj and x-> respectively (they are usually at different layers), then 
T, (x,) = A,x, + n, 

To (xo) = AoXo + Bo 

Hence, 

AT S = T, (x,)-T,,(xo) (fi) 
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wll<T0 


A| - jrj A . , B| — ( k | A i + It 1 1 1 )/li i 
n, - B| + f (A, - A,) 

A. — ( {2 ~ f I ) 

k 1 koh 1 + kilijhoL + kohjho^ + kjl^ha ~ kjli 1 li 2 ^ ^ 

S 

Note t hat AT’ is independent of time. 

Unsteady — State Solution 
l or unsteady state ease, we have 


, d 2, L't 
ai W 2 ~ 

<9T| 

0 < x < l, 



-n •• * 1 

~lh 

t >0 

(8) 

0 ^ T i 

a-2^72- 

= <9T> 
TTt~ 

0 < x < L, 

t > 0 



with 


1 . 0 1, 
-k * 7Tx 

h 1) | T| — It 1 1 

1 1 at x — 0 , 

t. > 0 

T, = 

T 2 at 

. x = f , t > 0 


dTi 

1 PTo . 

x = 7, t > 0 


li <sr 

~ k;i 


, (7V‘> 

Ko “T — 11 — 

“ ox 

+ lloTo = 1 1 2 f 2 

at x = L, t > 0 



and 

Tj (x,0) = T, (x,D) = T 0 


( 10 ) 


I he solutions to cqus. (7) and (8) with nonhomogeneous boundary conditions (9) and (10) 
can be expressed in t he form 

T. (x,t) = 0 ] (x,t) + <p. (x) f, + (x) f, (I!) 

i = 1 , 2 


where the functions v*(x), <A(x), and ft(x, t) are the solutions of the following subproblems: 
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(.1) £l J 


(Ix- 


=-- 0 


(I y 2 n 

<lx - ” - 0 


will) 



lWi 


= 1'ifi 


V>) = V 1 - 1 


ill X — 0 
at x = ( 


Therefore, 



at x = i 
at. x = L 


(b) 



V?. (x) = C ,x + D, 
V>> (x) = C'ox + Dj 


djj h 

dx 2 


= 0 


with 


i. s> v 

_Ki ar 

- + h|i/>i = 0 

at > 

i = 0 

ih = h 


at > 

; = f 

A , V 

i. d t - ‘» 

^ k > ■] — ~ — 
- (lx 

at > 

: = f 

(I I- f | 

1<2 dx 

+ h_>X'2 = h 2 l*2 

at > 

: = L 


we have 


V'i (x) = Ri x + F| 
4h (x) = E 2 x + V-> 


In this rase, the sensor outpui will he dependent on lime. To determine the sensoi 
output we calculate the temperature difference to be 
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_vr u (i.) = T| (x,i)-t,(x,,i ) 

= 0 [ (xt,t) — 0 >( Xj,l.)+ [ ^l(xi) — V>i(X2)]f| 
+ [Vi(xj) - V l -i(x>)]fo 


Computational results yliow that 


bi(xi) - v> 2 (x 2 )]fi + (V’i(xi) - &(x 2 )Jf 2 = AT S 
Therefore, 

AT u (t) = 0j(xi,t) - #>(xal) + AT S (12) 


(c) 0,(x,t) and 0,(x,t) are determined from the following homogeneous 
problem: 

(Pjh dO x 
a '0x 2 = Ot 

(P 02 00 , 

nj Ox 2 = 01 


dOx 


Subject to — k| dx + h|0| = 0 


at x — 0 

Ox = 0> 

at x = f 


00 1 00, 

kj Ox = bo Ox 

a t. x = f 


00, 

k _> Ox + 1 1 2 0 > = 0 

at x = I, 



and fA(x,0) = T q -^(x)I, - ^(x)!'. 

Jet ^(x,t) = x.(x) r(i) 

=> F(t) = e ~' j n '' 


118 


win to /in are eigenvalues which will Ik* determined later. 

f p 

X (x) = A Sin i — — x) + B Cos ( — - — x) 

• n > n fir ,n fir 
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x ,„ w = X,, Shl ( 


a. A, , 

(-p= x ) + B <7=-*) 

V a “ vi ® 


Without, loss of generality, we choose Ai» = I. R n , A (|J , and R (i can he det.ermined with 
the corresponding boundary conditionsfand several pages ol algebra). 

Finally, 


AT"(l) = S C n [X |H (x) - X, n (x 2 )]o^n l + AT S 

11 = 1 


(id) 


where C is determined by the initial conditions. The equation for the determination of 


the eigenvalues p is 

sin 7 + Hi 7 sin 7 —sin fj 

K cos 7— KHi7sin 7 cos // 

0 II>//cos ijV) + sin (jii) cos(^//)-U:*//sin(j//) 


-cos // 

—sin // =0 


where 7 = ^n ( , tj = f, Hi = 77 


f~f 


id? 


k| O) 


' V - yj n -’ 


1 2 

~ rrx 


Some additional alegcbra will lead to the following transcendtal equation lor the 
determination of the eigenvalues /in: 
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l-si •»'"J 

/n , 


h-, f , ,],-f 


irr 


ift(- 




/in) 


(M) 


Hewrii o equal ion ( 1 .1) as: 

AT u (t) = E C e — /Jn 1 + AT S 
n = I 

whore 


(15) 


C n = C n [X Ip (x i )-X 2 n*( x )] 

Numerical example shows that pr « ih 2 << A 2 << , so we can just lake the 

first term to approximate the result. In other words, we have: 

AT 11 = cie“ /,|2t + AT S 


Inil 


AT" = 0 


at t = 0 


Therefore, C| = — AT S 


Then, 


A'l 


ui 


A'l 


,s 



(Ih) 


where -i\ can be determined from Equation (II). The following observations can be made 
by examining Eqs. (14) and (Hi): 

(1) AT u (t) = AT S as t^ao; 


(2) The transient process depends on I lie boundary conditions (I 11 J 12 ), 
the geometry (f,L) and the properties ( k 1 - k > n 1 , n ) of the sensor. 

{:]) The transient term (e - ^ 1 *') may be very important if ii\ is small. 
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I epical Fxample 


To illustrate t.Iio impact of t hose results consider a typical example, Assume that, D 
- 0.8 cm, f = 0.02 cm, L = 0.11 cm, choose h>/hi = 1.0. See Fig. 1. The material of the 
sublayer is alumel whose properties are k > = V17 w/m-k, a-> = 0.7 x 10 ' nr/s. If the 
materials of the foil are pure nickel and nichrome (80% Ni , 20% Cr), respectively. The 
result s arc shown in Table I 




Table l 


Foil Material 

Oi (nr/s) 

Ih i 

_/j9| 

? ’ (at t = 0.05 second) 

Nickel 

(pure) 

2.3 x I0~ 5 

16.745 

8.15 x I0- 7 

Nichrome 3.1 

(80% Ni, 20% Cr) 

x 10 -6 

6.438 

1.26 x I0 _l 


As can be seen from Table I, if we choose Ni— Cr alloy as the foil material, t he 

a 2 , 

transient term "e 1 " would he important. 


SUMMARY 

YVe have init ialed a study of the t ransient heat flux response of heat flux sensors, 
i his paper presents a one dimensional analysis of the problem. This was done to elucidate 
the key parameters and trend for the problem. The results of a heat flux sensor depends on 
the thermal boundary conditions, the geometry and the thermal properties of the sensor. 

In particular, it was shown that if the thermal diffusivitv of the sensor is small, say for a 
material like Nichrome, then the transient behavior must, be taken into account. 
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The goal of this work is to develop and test thin-film thermocouples for 
space shuttle main engine (SSME) components. Thin-film thermocouples have been 
developed for aircraft gas turbine engines and are in use for temperature meas- 
urement on turbine blades up to 1800 °F (refs. 1 to 3). Established aircraft 
engine gas turbine technology is currently being adapted to turbine engine 
blade materials and the environment encountered in the SSME, especially severe 
thermal shock from cryogenic fuel to combustion temperatures. Initial results 
using coupons of MAR M-246(+Hf) and PWA 1480 have been followed by fabrication 
of thin-film thermocouples on SSME turbine blades. Current efforts are focused 
on preparing for testing in the Turbine Blade Tester at the NASA Marshall Space 
Flight Center (MSFC). Future work will include testing of thin-film thermo- 
couples on SSME blades of single crystal PWA 1480 at MSFC. 

Successful fabrication and testing of thin-film thermocouples on a coupon 
of SSME high-pressure fuel turbopump (HPFTP) material (DS MAR M-246(+Hf)) was 
discussed previously (ref. 4) . The same procedure, adapted from aircraft tech- 
nology (ref. 3), was used on an HPFTP blade. The test specimen was coated 
with 120 pm of NiCoCrAlY (PWA 270 Spec. - Chromalloy Corp.). An adherent, 
electrically insulating A 1 2 C >3 film was then grown for 50 hr. This thermally 
grown film was followed by a sputtered AI 2 O 3 film* Finally, the thin-film 
thermocouple of Pt/Pt(13% Rh) was deposited. A typical fabrication process on 
a plasma jet sprayed NiCrAlY-coated SSME blade is shown along with that for a 
jet aircraft gas turbine blade in figure 1. An example of a recently fabri- 
cated device is shown in figure 2. 

Currently work is underway to prepare for testing of the SSME blade in 
MSFC ' s turbine blade tester (shown schematically in fig. 3). We anticipate 
completing at least one round of testing during the next 5 months. The instru- 
mented blade will be inserted into a blade holder, supplied by MSFC (fig. 4), 
in the middle position and will be surrounded by two other SSME blades with or 
without instrumentation, depending on the test requirements. The test will be 
conducted with the instrumented blades in a wired holder installed in the tur- 
bine blade tester at either blade position A or B. Since only the middle 
blade experiences close to SSME conditions, only one blade per holder assembly 
will model SSME conditions during a run on the turbine blade tester. A mockup 
is being prepared with wiring to the installed instrumentation on the blades. 
Final qualification at MSFC will precede final design and building of a holder 
to accommodate the wiring and instrumented blades. 

In a related development, we have successfully instrumented coupons of 
single-crystal PWA 1480. It is anticipated that this material will eventually 
replace the current MAR M-246(+Hf). Thermal cycling data (fig. 5) show the 
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thin-film thermocouple on single-crystal PWA 1480 to be of comparable surviva- 
bility with thin-film thermocouple on coupons of MAR M-246 (ref. 4). Success- 
ful application of this technology to this blade material will allow rapid 

response to testing needs in the future for single-crystal PWA 1480 blades. A 
further enhancement of our ability to respond to future testing needs is our 

recent improvement of our in-house capability to fabricate, test, and charac- 

terize devices. We have expanded our analytical capability (ref. 4) with an 
energy dispersive x-ray spectroscopy instrument (which allows elemental identi- 
fication of materials). 
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THIN-FILM THERMOCOUPLE ON SSME BLADE 



Figure 2 
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HOLDER FOR TURBINE BLADE TESTER 
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Figure 4 


SIGNAL OUTPUT OF THIN-FILM THERMOCOUPLE 
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ABSTRACT 

Visible artifacts have been seen by observers of Space 
Shuttle Main Engine (SSME) test firings, either directly or 
on film, prior to several failures, some of which were 
catastrophic. A program was undertaken to attempt to 
identify, in the plume of an SSME, anomalous atomic — or 
molecular- -species resulting from wear, normal or excessive, 
of internal parts, allowing real time monitoring of engine 
condition or detection of failure precursors. To this end, 
measurements were taken during test firings on stand A1 at 
Stennis Space Center and on SSME, and other engines, at the 
Santa Suzanna Facility in California- The results indicated 
that a system having high spectral resolution, a fast time 
response, and a wide spectral range was required to meet all 
requirements, thus two special systems have been designed and 
built. One is the Optical Plume Anomaly Detector (OPAD) 
polychromator--a 16 channel spectroradiometer having 
adjustable optical bandwidth and center wavelength and a time 
response of lmS. This instrument has been described 
elsewhere. The other instrument is the superspectrometer, an 
optical multichannel analyzer having 8,192 channels covering 
the spectral band 250 to 1,000 nm. This instrument is 
described in this paper. 

INTRODUCTION 

During the process of testing Space Shuttle Main Engines 
(SSME), movie film surveillance is routinely employed. These 
films typically form the major source of information 
available after unusual occurrences. Film is used because 
rates such as 400 frames-per-second are not uncommon and high 
frame rates are needed to capture many of the anomalies 
observed. It is not unusual for "streaking” to occur for one 
frame only. "Streaking" is a visible (on film) phenomenon 
thought to result from extraneous material being injected, by 
some method, into the propellant streams within the engine. 
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The most commonly observed "streaking" appears to be at the 
periphery of the plume. It is assumed that anomalies 
occurring upstream of the Main Combustion Chamber (MCC) will 
be diffused throughout the plume and, therefore, also in the 
"mach disc." Anomalies occurring downstream of^the injector 
face plate may appear only at the plume edge. ' 

BACKGROUND 

It has been theorized that if such events can be seen on 
film, and indeed occasionally with the naked eye, then it 
should be possible to use electronic-optical detection to 
devise a system allowing automatic engine health monitoring. 
Such a system might permit the prevention of a catastrophic 
event. Much effort has been expended in background work in 
the attempt to build a data base on which such a detection 
system could be based. Multi spectral photometric 
measurements (Fig 1) were acquired over about a 2-year period 
at test stands at Stennis Space Center (SSC) in Southern 
Mississippi ' and at the SSME development stand (A3) -at' 
the Santa Suzanna Test Facility (SSTF) in Ca li fornia . 

Some gata was also obtained f ro m tests of other than 
SSME. Using relatively low spectral resolution equipment, 
several elemental signatures were found in the spectra 
acquired from tests. The base data acquired showed a very 
strong continuum from OH emission. That background is 
sufficient to cover most of the elemental spectra since those 
lines are very narrow and normally of low amplitude. The 
nature of the background places demands upon the equipment 
such that commercially available instruments can be used only 
to detect the most significant incidences. It would be 
preferable to be able to follow over a period of time several 
of the individual elemental or molecular spectral lines so 
that detectable variations from the norm may be properly 
interpreted . 


RESULTS AND OBSERVATIONS 

Calcium (Ca) was identified some time before a bearing 
cage failure occurred on an SSME; Ca was also seen during the 
sudden failure of a small H^/0„ engine having a compound 
combustion cycle. The increased level of Ca was detected as 
an erratic and significant increase in the normally observed 
level (believed to come from the ) . During another SSME 
test at SSTF an injector face plate carrying intentionally 
plugged LOX tubes suffered some damage and that damage 
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resulted in a "seeded" plume. Spectra indicated the presence 
of iron (Fe), nickel (Ni), chromium (Cr), and some copper 
(Cu). Other species were also noted in trace quantities. In 
truth, this event was a massive failure--the sort which would 
be extremely serious were the engine a flight unit. 

Another unplanned incident occurred at SSTF . Certain 
shipping procedures for SSMEs at SSC require that minor 
cracks in the MCC be sealed with a copper foil tape to allow 
maintenance of an inert dry purge. The tape was not removed 
at SSTF and when the engine was fired, the copper tape rather 
handily "seeded" the plume. This event was also detected 
easily. 


INSTRUMENT DESIGN 

The basic instrument for the Optical Plume Anomaly 
Detector (OPAD) program is called a Polychromator (PCM) ; 
for further base data a Superspectrometer (SSM) was 
designed. The PCM provides 16 individually adjustable 
spectrophotometer channels; the wavelengths are adjustable, 
the slit-widths are selectable, the amplifiers provide 
high-speed autoranging, the computers permit setting alarm 
limits or multichannel complex algorithyms for 
decision-making, and a second computer allows archival 
storage on removable magnetic disc. This system employs a 
quartz fiber-optic input system, allowing mounting of the 
optics and electronics inside the hard-core of the test 
stand, where the environment is friendly to non-ruggedized 
equipment. It is this instrument which is the prototype of 
the safety monitoring system. All the spectral lines 
accruing from previous failures are monitored as well as OH 
at 308 nM and two broad banded background channels to observe 
such things as spurious light and variable occlusions (such 
as vapor (H„0). The spectral band covered ranges from 250 
nM to 1,000 nM. This is a continuous device with a digital 
data system; the sample time is about 1.5 ms. 

The system speed is chosen with the fact in mind that the 
films indicate some events occur within 2.5 ms (400 
frame-per-sec film). It is not assumed that these events are 
individually significant, but rather may be indicative of 
normal wear. Confounding the ability to see these "streaks" 
is the fact that they cover a relatively small percentage of 
the total emitting area of the plume. Edge phenomena are 
therefore difficult to capture; mach disc inclusions might 
well be much easier to find, since the higher temperatures of 
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the disc cause elemental emission strength to increase. 
Experience with the system will tell more about the area of 
the plume which should be observed. 

The other instrument (the SSM) is intended to provide 
overall surveillance of the plume perchance some species for 
which the PCM has not been set appears. It had to have a 
wide spectral range, high data capability, and high spectral 
resolution. Since no commercial instrument provides all 
these features, a special one was designed. It Is based upon 
a 1/2 meter grating spectrometer, uses an input system 
similar to that of the PCM- -telescope and fiber optic 
cable — , and has at the exit slit folding mirrors directing 
the dispersed spectrum onto four 2,048 element self-scanning 
linear silicon diode arrays, for a total of 8,192 elements. 
The outputs are digitized and processed by computer for 
display and storage ( Fig 3 ) ; storage for archiving is by 
optical disc. 

The input system is composed of two telescopes ( Fig 2 ) , 
two quartz multifilament optical fiber cables. The numerical 
aperture of the fibers, projected onto the plume, does not 
quite cover the whole plume. The mounting point on the Test 
Bed Stand is about 63 feet from the "mach disc, " and at that 
distance the viewed field is about 6-1/2 feet across. Two 
input sets are used because one cable will illuminate only 
half the entry slit of the spectrograph. Using a direct 
viewing optic system would eliminate these difficulties of 
the fiber optics but would also require that the spectrograph 
and detector set be mounted on the engine deck in some way. 
All In all, the cable route was considered the most 
practical . The cables used are round at the telescope end 
and rectangular at the spectrograph entry slit in order to 
maximize light transfer at both ends; that scheme also 
isolates the field of view from the slit shape. 

The range of coverage is similar to that of the PCM--250 
nM to 1,000 nM, yielding an ultimate resolution of less than 
1A (0.1 nM). The sampling technique used, combined with the 
factory provided service boards ( for the arrays ) , allows at 
least 10 samples per second. The system normally will 
display each sample as a line spectrum; several modes would 
allow detailed examination of selected narrow ranges. 

The computer for this system is an "AT" style machine 
whose prime purpose is to store the data and display it. The 
archival data is placed on WORM (write once, read many) 
optical disc; because of the data rates the WORM drive 
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cannot accept the continuous stream of data, but the internal 
high capacity hard disc can. Therefore, the processed data 
is loaded into the hard disc, then on completion of the 
acquisition run, the hard disc is dumped to the WORM drive. 
Intensive data processing is done after-the-fact; too much is 
required of the computer to achieve real-time processing of 
this data. 


SUMMARY 

These two instruments should be sufficiently flexible and 
capable of acquiring spectral data to allow an in-depth study 
of the plume for most any reason. These systems will be used 
in an ongoing program wherein the proper algorithyms will be 
determined through experience. At first "simple" redline 
cuts are envisioned as the operational mode for the PCM, but 
more complex computations may be necessary. Each channel is 
presently corrected for background conditions, but 
multichannel dependence could drive a change to the basic 
philosophy. 

Also, the PCM could be considered a conceptual (from a 
systems, or scientific, standpoint) prototype for a flight 
unit. This system, as such, cannot fly — it is much too 
heavy, bulky, and fragile to put on a shuttle--but the 
concept is valid, whatever the configuration of a final 
system. Also, the problem of viewing the proper place in the 
plume affects the design; the varying nature of the plume 
with respect to altitude forces considerable compromises on 
the system. For example, it may be desirable to view the MCC 
(Main Combustion Chamber) directly by some method. The 
Technology Test Bed operation should indicate what will be 
required and how to affect it. 
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. Maximum dock rate » 100 KHZ. 
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PROBABILISTIC STRUCTURAL ANALYSIS - 
INTRODUCTORY REMARKS 

Christos C. Chamis 
NASA Lewis Research Center 


It is becoming increasingly evident that deterministic structural analy- 
sis methods will not be sufficient to properly design critical components in 
hot engine structures. These structural components are subjected to a variety 
of complex and severe cyclic loading conditions, including high temperatures 
and high temperature gradients. Most of these are quantifiable only as best 
engineering estimates. These complex loading conditions subject the material 
to coupled nonlinear behavior which depends on stress, temperature, and time. 
Coupled nonlinear material behavior is nonuniform, is very difficult to deter- 
mine experimentally, and is perhaps impossible to describe deterministically. 

In addition, hot rotating structural components are relatively small. Fabrica- 
tion tolerances on these components, which in essence are small thickness vari- 
ations, can have significant effects on the component structural response. 
Fabrication tolerances by their very nature are statistical. Furthermore, the 
attachment of components in the structural system generally differs by some 
indeterminate degree from that assumed for designing the component. In sum- 
mary, all four fundamental aspects - loading conditions, material behavior, 
geometric configuration, and supports - on which structural analyses are based, 
are of a random nature. The direct way to formally account for all these un- 
certain aspects is to develop probabilistic structural analysis methods where 
all participating variables are described by appropriate probability functions 
(fig. 1) . 

The development of probabilistic structural analysis methodology consists 
of the following program elements: (1) composite load spectra models, (2) com- 

putational probabilistic structural analysis methods, and (3) probabilistic 
constitutive relationships. The development of the probabilistic structural 
analysis methodology is a joint program of NASA Lewis in-house and sponsored 
research. The objective of this session is to illustrate recent progress on 
the application of this methodology to determine the reliability of structural 
components for rocket propulsion systems. 

The session contains descriptions of and progress reports on the follow- 
ing specific activities: (1) The NESSUS computer code, (2) approximate meth- 

ods, (3) advanced methods, (4) composite load spectra applications, (5) proba- 
bilistic fracture mechanisms, and (6) probability of failure and risk analysis. 
Collectively, the progress to date demonstrates that the structural durability 
of hot engine structural components can be effectively evaluated in a formal 
probabilistic/reliability framework . 


PRECEDING PAGE BLANK NOT FILMED 
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UNCERTAINTIES IN LOAD CONDITIONS AND STRUCTURAL 
PRIMITIVE VARIABLES ARE INTEGRATED BY PERFORMING 
MULTIPLE FINITE ELEMENT ANALYSES 


GLOBAL: AT ALL NODAL POINTS 

[M][uj + [c]fuj + [k](uj = [F(t)j 


LOCAL: AT EACH REGION BOUNDED BY 

[M] = J v [B M ] T [D][B M ]dv 

[C] = J v [B c ] T [q[D][B c ]dv 

[K] = J v [B k ] T [D][B]dv 


THE UNCERTAINTIES IN THE 
STRUCTURAL RESPONSE VARIABLES 
(u, w, o, e, G) ARE DETERMINED FROM 
THE RESULTS OF THE MULTIPLE FINITE 
ELEMENT ANALYSES 


^-NODAL LINE 
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Figure 1 


STRUCTURES RESEARCH ACTIVITIES ON 
SSME STRUCTURAL DURABILITY 



RESEARCH ACTIVITY 

TYPE 

INSTITUTION 

B16 

COMPOSITE LOAD SPECTRA 

CONTRACT 

ROCKETDYNE (PRIME) 

B15(A) 

PROBABILISTIC STRUCTURAL 
ANALYSIS METHODS (PSAM) 

CONTRACT 

SWRI (PRIME) 

B15(G) 

PROBABILISTIC FINITE ELEMENTS 
FOR FRACTURE MECHANICS 

GRANT 

NORTHWESTERN UNIV. 

B15(C) 

PROBABILISTIC STRUCTURAL 
ANALYSIS 

IN HOUSE 

SVERDRUP SUPPORTED 
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Figure 2 
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PROBABILISTIC STRUCTURAL ANALYSIS METHODS FOR 


SELECT SPACE PROPULSION SYSTEM COMPONENTS 

H.R. Millwater and T.A. Cruse 
Southwest Research Institute 
San Antonio, Texas 
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The Probabilistic Structural Analysis Methods (PSAM) project developed at SwRI integrates 
state-of-the-art structural analysis techniques with probability theory for the design and analysis of 
complex large-scale engineering structures. PSAM is currently in the fifth year of a five year 
contract. An advanced efficient software system (NESSUS) capable of performing complex 
probabilistic analysis has been developed. NESSUS contains a number of software components to 
perform probabilistic analysis of structures. These components include: an expert system, a 
probabilistic finite element code, a probabilistic boundary element code and a fast probability 
integrator. The NESSUS software system is shown in figure 1. 

An expert system is included to capture and utilize PSAM knowledge and experience. 
NESSUS/EXPERT is an interactive menu-driven expert system that provides information to assist 
in the use of the probabilistic finie element code NESSUS/FEM and the fast probability integrator 
(FPI). Figure 2 summarizes the expert system menu structure. 

The NESSUS system contains a state-of-the-art nonlinear probabilistic finite element code, 
NESSUS/FEM, to determine the structural response and sensitivities. A broad range of analysis 
capabilities and an extensive element library is present. A summary of the current analysis 
capabilities include static, eigenvalue, harmonic, random vibration, and buckling analyses. A variety 
of random variables are available such as geometry, loads, material properties, temepratures, 
boundary conditions, harmonic excitation, power spectral density input, etc. In addition, special 
algorithms have been developed and implemented to efficiently perform sensitivity computations 
required for probabilistic structural analysis. A summary of the current code capabilities is presented 
in figure 3. 

A fast probability integrator (FPI) has been developed to efficiently perform probabilistic 
analysis. Fast probability integration is an approximate technique which offers significant 
advantages over traditional Monte Carlo methods particularly in the low probabilities of failure 
typically found in structural engineering. FPI contains a number of features and options which are 
summarized in figure 4. 

A probabilistic boundary element analysis program (PBEM) is available in the NESSUS 
system as an alternative to the finite element code. For many types of analyses, the boundary 
element method offers significant efficiency and accuracy advantages over finite elements. PBEM 
is a robust three-dimensional nonlinear boundary element code. The current capabilities of the code 
are: linear elastic and plastic analyses with temperature dependent materials, time dependent loading, 
and normal mode vibration. 
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NESSUS has been validated on a number of small-scale problems where the exact solution 
can be determied. Figure 5 is a finite element model for a notched plate. The primary random 
variable is the radius of the notch and the response variable is the probabilistic stress concentration. 
Figure 6 demonstrates the random vibration capabilities of NESSUS. The random variables include 
the Power Spectral Density acceleration, the structural damping and the PSD cutoff frequency. The 
response variable is the root-mean-square tip displacement. In both cases, excellent agreement is 
found between the NESSUS solutions and exact solutions. In general, NESSUS has been validated 
for a number of capabilities including: static stress analysis, eigenvalue, buckling, harmonic 
excitation and random vibration. 

The application of NESSUS to SSME components is being conducted under the verification 
phase of PSAM. The four selected components are shown in figure 7. Each component exercises 
adifferent capability of NESSUS. Presently, verification studies of the turbine blade and the HPOD 
have been completed. The turbine blade model exercises the static stress analysis and natural 
frequency capabilities of NESSUS. The HPOD model exercises the random vibration capabilities. 

Graphical interfaces between NESSUS and pre- and post-processors have been developed to 
enhance data interpretation. Figure 8 is a contour plot of the probability of exceeding 80,000 psi 
stress for the turbine blade. This plot can be used to quickly determine critical areas which can 
then be examined more accurately with further computations. 

NESSUS is currently being extended into transient and nonlinear material and geometric 
problems. Random variables will include initial accelerations, displacements, velocities and 
nonlinear parameters such as yield stress and hardening parameters. 


Overview of the NESSUS System 



NESSUS/EXPERT Menu System 
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Verification Analyses 
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AN APPROXIMATE METHODS APPROACH TO PROBABILISTIC STRUCTURAL ANALYSIS 

R. C. Me Clung , H. R. Mlllwater, Y.-T. Wu, B. H. Thacker, and 0. H. Burnside 

Southwest Research Institute 
San Antonio, Texas 

A major research and technology program in Probabilistic Structural Analysis 
Methods (PSAM) is currently being sponsored by the NASA Lewis Research Center 
with Southwest Research Institute as the prime contractor [1,2]. This program 
is motivated by the need to accurately predict structural response in an environment 
where the loadings, the material properties, and even the structure may be 
considered random. The heart of PSAM is a software package which combines 
advanced structural analysis codes with a fast probability integration (FPI) 
algorithm for the efficient calculation of stochastic structural response. 

The most common structural analysis techniques in use today are finite 
element and boundary element methods. These techniques permit highly accurate 
three-dimensional modeling of structural geometry, thermal and mechanical load 
environment, and material properties. Unfortunately, these techniques can also 
be expensive and time-consuming. In the early stages of design, when geometries, 
loads, and materials are only tentative, it is not practical to assemble an 
exact, comprehensive three-dimensional model for every critical component. What 
Is needed, typically, is some means of estimating the nature and rough magnitude 
of stresses, displacements, natural frequencies, etc. Of further value is the 
identification of those design parameters which exert the most influence on the 
total system performance, so that further design evolution is more efficient. 

These needs are addressed in the PSAM software by the Probabilistic 
Approximate Analysis Methods (PAAM) module. The basic idea of PAAM is simple: 
make an approximate calculation of system response, Including calculation of the 
associated probabilities, with minimal computation time and cost, based on a 
simplified representation of the geometry, loads, and material. The deterministic 
solution resulting should give a reasonable and realistic description of 
performance-limiting system responses, although some error will be inevitable. 
If the simple model has correctly captured the basic mechanics of the system, 
however, including the proper functional dependence of stress, frequency, etc. 
on design parameters, then the response sensitivities calculated may be of 
significantly higher accuracy. In other words, the calculated probabilistic 
distribution of the response variable may be in significant error only by some 
offset of the mean value. 

Three factors make up the ’’approximate" analysis approach. The first is a 
simplified representation of the part geometry. Complex three-dimensional shapes 
are replaced by simple beams, shells, etc., with relatively few descriptive 
parameters. The second factor is a similar simplification of the applied 
mechanical and thermal loads. Point loads, uniform fields, and linear or parabolic 
distributions are used to describe these quantities. Often a complex load 
environment can be approximated through linear superposition of the simple 
descriptors. The third factor in the PAAM methodology is the solution technique 
itself. Complex numerical techniques with many degrees of freedom are replaced 
by simpler computational schemes. In some cases, it is possible to construct 
a simple mechanics -of -materials model, often in two dimensions. In other cases, 
more sophisticated closed-form solutions can be derived or (preferably) adapted 
from previous research results in the literature. These may be based on elasticity 
or may employ approximate energy methods. 
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The. PAAM code employs the fast probability integration (FPI) algorithm of 
Wu and Virsching [3], which has been demonstrated to consistently provide fast 
and accurate estimates of point probabilities for typical engineering response 
functions. This algorithm identifies a most-probable -point (or design point) 
on the response surface, establishes a quadratic polynomial approximation to the 
response function at that point, and then transforms the quadratic form into a 
linear one. Input variables are defined using probability distributions. An 
optimization routine is employed to approximate non-normal variates as equivalent 
three - parameter normals, thereby approximating the limit state as linear in 
normally distributed design factors. Solution of this simplified problem is 
straightforward. The FPI Is repeatedly applied to a number of limit states to 
compute the cumulative distribution function (CDF). 

The PSAM code is being demonstrated and validated by analyzing four 
representative critical components In the current Space Shuttle Main Engine 
(SSME). The same four components are considered in the PAAM code, which by its 
very nature is component- specific . It should be noted, however, that the 
architecture of the PAAM software permits new closed- form expressions for these 
or other components to be installed and evaluated quickly. In this presentation 
one of these demonstration components, the liquid oxygen (LOX) post, Is used to 
illustrate typical PAAM solution strategies and results. 

Two different simple models are used to analyze the LOX post. The steady 
temperature and pressure gradients between the internal and external radii are 
evaluated by modeling the post as an axisymmetric thick cylinder. This is a 
classic elasticity problem, and exact solutions are readily available In standard 
textbooks. Other loadings, including the transverse fluid flow, are considered 
by modeling the post as a beam with hollow circular cross-section and elastically 
restrained ends (i.e., a mechanics-of-materials approach). Note how this model 
Is used to solve a vibration problem. The basic solutions for free vibration 
of fixed- fixed and fixed-pinned uniform beams are available in handbooks. 
Solutions for a limited number of intermediate end conditions are available In 
tabular form. These results can be fitted with a simple empirical relationship 
to estimate the frequency factor for any mode and any end condition. Subsequent 
calculations of mode shapes and stresses for forced vibration follow directly 
from simple equations. 

Several sample PAAM/FPI calculations have been performed in order to validate 
and further demonstrate the approximate methods approach. Results are given In 
terms of cumulative distribution functions for critical response variables. FPI 
and Monte Carlo calculations are shown to agree closely. Also determined are 
the sensitivity factors of the random variables. These sensitivity factors are 
a function not only of the mathematical form of the equation, but also of the 
distribution types and coefficients of variation of each random variable. They 
are of particular value for guiding improvements in the analytical scheme or 
component design. 
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AN OVERVIEW OF PROBABILISTIC STRUCTURAL ANALYSIS METHODS 


0 PSAM LINKS STRUCTURAL ANALYSIS 
WITH FAST PROBABILITY INTEGRATION 
FOR EFFICIENT CALCULATION OF 
STOCHASTIC STRUCTURAL RESPONSE 
WHEN LOADS , MATERIALS, AND 
GEOMETRY ARE RANDOM 

0 FEM AND BEH ARE MOST COMMON 
STRUCTURAL ANALYSIS TECHNIQUES 

- DETAILED MODELS 

- ACCURATE 3-D SOLUTIONS 

- EXPENSIVE, TIME-CONSUMING 

- NOT APPROPRIATE FOR 
EARLY DESIGN/ANALYSIS 

o Need a rapid, inexpensive estimate 
OF structural response which 

INCLUDES PROBABILISTIC INFORMATION 



o "Probabilistic Approximate Analysis Methods" (PAAM) 


WHAT IS AN "APPROXIMATE METHOD" ? 


o Simplified representations of 

- component geometry 

- APPLIED LOADS 

- material properties 

o Rapid, efficient solution scheme 

- mechanics-of-materials 

- simple energy methods 

- closed form or tabulated 

literature solutions 

o Calculated mean value of response 
function contains some error 

o If basic mechanics are correct, 

RESPONSE DISTRIBUTIONS AND 
PARAMETER SENSITIVITIES WILL 
BE ACCEPTABLY ACCURATE 
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PAAM CODE USES FAST PROBABILITY INTEGRATION ALGORITHM 


0 Wu-WlRSCHING FPI ALGORITHM 

o Shared by other PSAM modules 
o Identify MPP on response surface 

FOR GIVEN PROBABILITY LEVEL 

o Establish quadratic approximation 

TO THE RESPONSE FUNCTION AT MPP 
AND TRANSFORM TO A LINEAR FORM 

o Input variables defined using 

PROBABILITY DISTRIBUTIONS 

o Approximate non-normal variates 

AS EQUIVALENT 3-PARAMETER NORMALS 

o Solve the resulting simple problem 
(linear with normal design factors) 


JOIOT PROB. DENSITY IfX) 



PROB I T « WOJ . J^f »(X)X {PROEL OF FAILURE) 


PAAM CODE DEMONSTRATED WITH SSME LOX POST 


o One of four SSME 

DEMONSTRATION COMPONENTS 

0 A SLENDER, HOLLOW TUBE 
FIXED AT BOTH ENDS 

o Important loadings include 

- AT AND A p ACROSS TUBE WALL 

- STEADY AND DYNAMIC LOADS DUE 

TO TRANSVERSE FLUID FLOW 

0 TWO DIFFERENT PAAM MODELS SHOWN 

- THICK CYLINDER MODEL 

- BEAM MODEL 

o Although these equations 

ARE COMPONENT-SPECIFIC, 

THE BASIC PAAM ALGORITHM 
IS COMPONENT-INDEPENDENT 


ROW 
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THICK CYLINDER MODEL USES KNOWN ELASTICITY SOLUTIONS 


o Thick cylinder subjected to 

DIFFERENT INTERNAL AND EXTERNAL 
PRESSURES AND TEMPERATURES 

o Exact closed-form solutions 

AVAILABLE IN STANDARD TEXTBOOKS 

o Consider the hoop stress: 

- PRESSURE SOLUTION 

(i o/' ) ? (p 0 - p,) p, - p„k ? 

°° k ? - 1 + k ? ~\ 

- THERMAL SOLUTION 

f5 (r„//) ? +l 1 - 1 n ( r „ / r ) 

° 0 = " P F-i + P “ mot — 

a T 7 (7~, - 7~ a ) 

J?d-v) 
k = r a /r, 



BEAM MODEL USES A MECHANICS-OF-MATERIALS APPROACH 


o Slender beam with hollow cross-section 

- FIXED AT ONE END 

- ELASTICALLY RESTRAINED AT THE OTHER 

o Consider the vibration problem 
o Basic equations in handbooks 

IN TERMS OF A FREQUENCY FACTOR \ 

o Lambda depends on end fixity and mode number 

- closed form expressions for fixed, pinned 

- tabulated solutions for elastic restraint 

- FIT THESE WITH AN EMPIRICAL EXPRESSION 

o Simple to calculate natural frequencies, 

DISPLACEMENTS AND STRESSES DUE TO 
HARMONIC AND RANDOM FORCED VIBRATION 


\\VYv\Vi 
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PAAM INPUT INCLUDES 

STATISTICAL DESCRIPTIONS OF RANDOM VARIABLES 


Input Variables for Lox Post Thick Cylinder Model 


Variable 

Distribution 

Mean 

COV 

r ( 

truncated Normal 

0.094 In. 

1.06% 


(*0.003) 



>0 

truncated Normal 

0.110 in. 

4.55% 


(-0.002, 40.01) 



F 

Normal 

3.40E47 psl 

2% 

V 

Normal 

0.3594 

2% 

a 

Normal 

5.65F.-6 /R 

5% 

n> 

Lo gnormal 

3077 psl 

4% 


T.ognormal 

3232 psi 

4% 


Lognormal 

194 R 

1.55% 

LjL- 

l.o gnonnal 

1444 R 

1.55% 


o User supplies mean value, 

DISTRIBUTION TYPE, AND COV 
FOR EACH RANDOM VARIABLE 

o Many distribution types 

POSSIBLE (INCLUDING TRUNCATED) 

o Set C0V=0 for 

DETERMINISTIC VARIABLES 

0 If STATISTICAL INFORMATION 
UNKNOWN, MAKE ESTIMATES 
FROM DEFAULT SPECIFICATIONS 


PAAM OUTPUT INCLUDES CDF OF RESPONSE FUNCTION 


O FPI CALCULATES CUMULATIVE 
DISTRIBUTION FUNCTION (CDF) 
OF RESPONSE VARIABLE 


0 CDF BASED ON EQUIVALENT 
NORMAL DISTRIBUTION 
OF RESPONSE VARIABLE 


o Excellent agreement between 
FPI and Monte Carlo 


■H 'i'Jt.Hl'. r — — 





U 



ii ( « > p stress inn 
- MONT! t'AHU) 


.’WiS 


O FPI MUCH FASTER THAN MONTE CARLO 


GR.GJi'JAL PAGE IS 

OF POOR QUALITY 
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PAAM CALCULATES SENSITIVITY FACTORS OF RANDOM VARIABLES 


LOX POST (THICK CYLINDER) 



ri no i nu alpha pi po n to 


RAJHXJM VARIAJU 


o Shows relative influence of 

RANDOM VARIABLES ON THE TOTAL 
UNCERTAINTY OF THE RESPONSE 

o Function of both 

- MATHEMATICS OF THE EQUATION 

- DISTRIBUTION TYPE AND COV 

OF THE RANDOM VARIABLES 

o Helps the designer/analyst to 

DETERMINE REQUIREMENTS FOR 
FURTHER DESIGN/ANALYSIS 

- MORE ACCURATE INPUT STATISTICS 

- MORE EXACT ANALYSIS MODEL 

- TIGHTER TOLERANCES IN DESIGN 

- BEST OPPORTUNITIES FOR 

DESIGN IMPROVEMENTS 


SUMMARY 


O PAAM MAKES A RAPID, INEXPENSIVE ESTIMATE OF THE 
PROBABILISTIC STRUCTURAL RESPONSE OF A SYSTEM 

o Requirements: 

- SIMPLIFIED REPRESENTATION OF COMPONENT GEOMETRY AND LOADS 

- CLOSED- FORM SOLUTIONS DERIVED OR ADAPTED FROM THE LITERATURE 

- STATISTICAL DISTRIBUTIONS OF INPUT VARIABLES 

o Probabilistic calculation performed with FPI algorithm 
o Output includes: 

- CDF OF RESPONSE VARIABLE 

- SENSITIVITY FACTORS FOR RANDOM VARIABLES 

0 PAAM PARTICULARLY USEFUL IN EARLY STAGES OF DESIGN/ANALYSIS 
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PROBABILISTIC BOUNDARY ELEMENT METHOD 

T.A. Cruse and S.T. Raveendra 
Southwest Research Institute 
San Antonio, Texas 
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The purpose of the Probabilistic Structural Analysis Methods (PSAM) project is to develop 
structural analysis capabilities for the design analysis of advanced space propulsion system 
hardware. The PSAM effort consists of three major technical thrusts: probabilistic finite element 
methods (PFEM), probabilistic approximate analysis methods (PAAM), and probabilistic advanced 
analysis methods (PAD AM). The boundary element method is used as the basis of the Probabilistic 
Advanced Analysis Methods (PADAM) this is discussed in this text. In particular, the BEST3D 
code developed under NASA/HOST program is modified for inclusion in the PSAM module as 
NESSUS/BEM. 

The probabilistic boundary element method code (PBEM) is used to obtain the structural 
response and sensitivity results to a set of random variables. As such, PBEM performs analogous 
to other structural analysis codes such as finite elements in the PSAM system. The probabilistic 
analysis is performed by coupling PBEM and the Fast Probability Integrator (FPI) using a highly 
efficient mean based algorithm. Preliminary validation studies have shown PBEM to be an accurate 
tool for probabilistic analysis. 

For linear problems, unlike the finite element method (FEM), the BEM governing equations 
are written at the boundary of the body only, thus, the method eliminates the need to model the 
volume of the body. However, for general body force problems, a direct condensation of the 
governing equations to the boundary of the body is not possible and therefore volume modelling is 
generally required. Since such volume modeling mostly eliminates the advantage of the BEM 
procedure, a surface transformation technique based on particular integrals is used to replace the 
volume integral by equivalent surface integrals in the current analysis. 

To illustrate the particular integral procedure, consider the solution of thermoelastic problem. 
The surface transformation technique requires the evaluation of the particular integrals of the 
inhomogeneous Navier’s equilibrium equations. However, in general, the temperature field is not 
known as a continuous functions, instead, the values of the temperatures are known at selected 
points at the surface and volume of the body. The procedure assumes that the temperature field 
then can be expressed in terms of interpolation functions and unknown densities associated with 
each collocation point. Since this procedure replaces the actual temperature field by an approximate 
field, the resulting integrals equations, while satisfying the equilibrium conditions, are only 
approximate. A plot of the interpolated temperature field on the surface of a cube subjected to 
uniform temperature shown in figure 4 indicates that the global interpolation results is in error 
except at collocation points. Nevertheless, the particular integral solution procedure uses nodal 
values interpolated by regular isoparametric shape functions in their computations which seems to 
minimize the interpolation error effect. 


preceding page blank not filmed 
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As a first example, consider the free vibration analysis of a cantilever beam. The first five 
normal modes shown in figure 5 using two BEM models agree well with FEM results, thus, validating 
the free-vibration analysis procedure. To validate the particular integral procedure for the 
thermoelastic problem, we considered the solution of hollow sphere with an external radius to 
internal radius ratio of 2 subjected to linear radial temperature variation. The normalized hoop stress 
compared to the theoretical results in figure 6 validates the particular integral procedure for thermal 
problems. To further investigate the accuracy of the procedure for higher temperature variations, 
we considered the same sphere subjected to cubic radial temperature variation. Again the agreement 
between the computational the theoretical results for hoop stresses is excellent as shown in figure 
7. This shows that while the interpolated temperature field is in error, the solution of the 
corresponding boundary integral equations even for a higher order temperature field is still very 
accurate. As a final example, we considered a hollow cylinder subjected to linear temperature 
variation. We further assumed that the Young’s modulus also varied linearly for this case. The 
resulting hoop stress compared to FEM in figure 8 validates the material inhomogeneity procedure. 


Overview of NESSUS/BEM System 
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BEM Formulation for General Body Force Analysis 
Does Not Use Domain Modeling 


• Regular Domain BIE for Thermoelasticity 

cu + J Tuds = J Utds+f S6dv 

__ j s __ j s __ j v _ 

• BIE for Particular Solution Replaces Domain Integral 

cu p +f T u p ds - J U t p ds = f S9dv 
s __ s~ J v~ 

• Surface-Only BIE Replaces Thermoelastic BIE 

A f A f A 

cu + J Tuds=J Utds 
s s 

A n 

— u M (0) 
t = t - t P (0) 

Particular Solutions Obtained for Approximating 

Body Forces 

• Navier’s Equilibrium Equation (Thermoelastic) 

N • u p = mVQ 

• Collocate fi with Boundary Interpolation Functions 

N 

K(q m ,Q n ) * *0 (Qri^ Q n rS 

n=1 

• Particular Solution for Collocated Temperatures Replaces Domain 
Integral 

cu p +f T u p ds - f llPds^J SOdV 

__ j s __ — j v _ 

• Surface-Only BIE for Interpolated Temperatures 



Interpolation Error Within the Domain 


Error in Global Interpolation 
Function for Constant Temperature 

• Global Interpolation Functions 


• Collocated at User-Defined 
Points (Surface and Interior) 

• Particular Solution Method 
Minimizes Interpolation 
Error Effect 

• Particular Solution 
Approximation Uses BIE 
Shape Functions 



Collocation Points 



Error 


Normal Vibration Modes of BEM and FEM 
Models Are Compared 



a) BEM Map 1 




b) BEM Map 2 



b) FEM Map 2 


First Five Modes of Vibration for the Cantilever Models 



Mode 1 

Mode 2 

Mode 3 

Mode 4 

Mode 5 

Nodes/ 

Elements 

CPU (sec) 
VAX 8700 

BEM-1 

2686 

5441 

13485 

16810 

25382 

44/14 

352 

BEM-2 

2810 

5333 

12736 

16255 

24523 

86/28 

1676 

FEM-1 

3033 

5449 

13191 

17227 

25141 

195/96 

158 

FEM-2 

_ . ... 

2885 

5354 

12778 

16224 

24522 

1125/768 

801 
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Results Are Accurate 
Even For Cubic Temperature Variation 



v = 0.3 
a = 10“ 4 


T = 100 + lOr + r 2 + 0.2r 3 


Temperature-Dependent Thermal 
Analysis Procedure Validated 


Hollow Cylinder 



b/a = 2 

a = 10 

E 0 = 30 x 10 6 

v = 0.3 
a = 10 -4 



T = lOr 

E = E 0 - 10 4 T 
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Current NESSUS/BEM Code Capabilities 


• Linear Elastic Stress Analysis 

• Centrifugal Loading Analysis 

• Free-Vibration Analysis 

• Thermal Analysis 

• Thermal Analysis with Temperature-Dependent 
Material Properties 


NESSUS/FEM Future Developments 


• Additional Particular Integrals 

(e.g.. Steep Thermal Gradient, Linear Temperature Field) 

• Elastoplastic Analysis 

• Large Deformation Analysis 

• Transient Analysis 
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Composite load Spectra 

[or Select Space Propulsion Structural Components'* 


H. Ho arid J. I. Newell 

Rockwell International Corporation, Rocketdync Division 

R. f . Kurth 

Bat tell Columbus laboratory 
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the objective of the Composite load Spectra (CIS) project is to build a 
knowledge based system to synthcsi/e probabilistic loads for selected space 
propulsion engine components, the knowledge based system has a load expert 
syslem module and a load calculation module, the load expert system provide- 
load information and the load calculation module generates the probabilistic 
load distributions. 

The engine loads can be divided into four broad classes. the engiw 
independent loads, the engine system dependent load, the component local 
independent loads arid the component loads. The engine independent loads 
include engine operating conditions and the engine hardware parameters. Ihe 
engine operating conditions are the 0/1 mixture ratio, and the engine fuel and 
oxidizer inlet conditions etc. the engine hardware parameters .ire for example 
the HPMP pump efficiency, the HPOIP turbine efficiency, etc. the engine 
system dependent loads are the operating conditions of the engine subsystems 
such as the 11PPIP turbine inlet and discharge pressures. The component loads 
including the local independent loads are loads internal to the engine 
subsystems, the component loads are evaluated with a multi level engine model 
implemented on CI.S. The multi-level engine model includes an engine 
probabilistic influence model and many component load models. the engine 
subsystem dependent loads are evaluated by the probabilistic influence model 
and they are in turn fed into the component load models lo generate the 
component loads, figures 1 and ?. illustrate the loads and the engine models. 

the engine probabilistic influence model calculates the engine system 
dependent load variations from a set of engine influence cocflicienfs as shown 
i n I i <;j u r c 3 , If i o i n I 1 u o ri r; c c o c C I i c i o n l s < i r c d c j.) <: n 1 1 e n l of i f i c c n y i n o c onirnd n d e d 

power level, these influence coefficients were derived from the engine (light 
and test data. The engine model implemented on CIS has f/i independent loads 
and 99 system dependent loads. I he load mean values, coefficient oT 

variations and the influence coclicicnts are stored in the CIS knowledge 
base, they ran be retrieved by the load expert system when needed. Iberc are 
many component load models as listed in figure 'I, each of which is for a 
certain component load. lor example, the pressure scaling model is for 
evaluation the component static pressure for the turbine blade, the Irons!, r 
duct and other components. Ihe probabilistic thermal load model is fi r 
evaluating the temperatures of the components. tri many cases, each compere a l 
has its own geometry model and related parameters. I hey are stored in a 
database format in the CIS knowledge base. 


* Work performed under NASA Contract NASS -?/\MYS 
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The four engine components selected for implementation are the turbine blade 
the transfer duels, the LOX posts and the HPOIP discharge duct. The loads for 

these components important to the probabilistic structural analyses arid 

reliability calculations are listed in figure 5. 

The load calculation module provides three probabilistic methods for 

generating probabilistic loads. They arc the Gaussian moment method, the 
RASCAL (Random Sampling Condensation Algorithm) method and the Monte’ Carlo 
Method. The Gaussian moment method assumes all input loads arc normally 
distributed. The dependent loads are evaluated analytically, the RASCAL is a 
discrete probability distribution method. It reduces computational time by 
restricting the sample space as it combines loads, tiic Monte Carlo method is 
a standard version of the Monte Carlo sampling method. The probabilistic 

models implemented on CI.S are shown in figure 6. they include models for 
calculating the steady state load and the transient stale loads. At present, 
the steady state and the quasi-steady slate models arc available for duty 

cycle mission calculations. The other probabilistic models are available for 
stand alone calculations. figure 7 shows a comparison of the results of the 
Gaussian method arid the RASCAL method with the HPI-IP discharge pressure data 
at 100% power level. The calculations match reasonably well with the data, 
figures 8 and 9 show the confident interval calculations for the HPflP turbine 
torque using RASCAL method with various bin numbers and different sampling 
options. The confident intervals compare well with the Monte Carlo method. 
Notice that the bin numbers used iri RASCAL are much smaller than the Monte 
Carlo sample size. 

Ihe Composite load Spectra knowledge based system as delineated in figures 10 
and II are developed for load syntheses of selected space propulsion engine 
components. Its domain knowledge includes engine load information and 

probabilistic engine load models. Ihe knowledge based system has a built-in 
database system that can handle the database operations. It provides three 
basic probabilistic methods and numerous probability models for building 
models for component load calculations. Ihe load expert system can retrieve 
the required load data and prepare input files for the load calculations. In 

an application study, presented separately in this conference, we have 
demonstrated that using the marginal distribution method the CIS can provide 
the composite load spectra to a structural analysis program for probabilistic 
analysis arid thus eliminates the need of modeling the loads inside the 
structural analysis program. 


ORIGINAL PAGE IS 
OF POOR QUALITY 
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SYSTEM CLASS OF LOADS 


• STEADY STATE 

• DIRECT VARIABLES — THRUST, INLET CONDITIONS 

• RANDOM VARIABLES - HARDWARE. TEST VARIATIONS 

• TRANSIENTS 

• START ANO CUTOFF CONTROLLED 

• LOCAL EFFECTS 


DIRECT VARIABLES- 



INJECTOR - RANDOM VARIABLES 
• FLOW LOSSES 
* LOXDOME 
v • HOT GAS 



DUCT - RANDOM VARIABLES 
FLOW LOSSES 


TURBOPUMP RANDOM 









C0MPD5ITE LOAD SPECTRA 
WULTr-L£ V€L ENGINE MODEL 


ENGINE SYSTEM MODEL 


PROBABILISTIC INFLUENCE MODEL 


Jo 


I (IC) 
1 


J1 


lo 


WHERE X^s ARE ENGINE INDEPENDENT LOADS 


f-g. ENGINE MIXTURE 

FUEL AND OXIDIZER INLET PRESSURES AMO INLET TEMPERATURES 

HPFTP 4 HPOTP PUMP EFFICIENCIES 

HPFTP 4 HPOTR TURBINE EFFICIENCIES, ETC. 

Yj's ARE ENGINE SrSTEH DEPENDENT LOADS 

e.g. HPFTP TURBINE INLET PRESSURE AND DISCHARGE PRESSURE 
HPOTP TURBINE TORQUE, ETC. 

(IC)^ ARE THE INFLUENCE COEFFICIENTS 

THEY ARE FUNCTIONS OF THE COMMANDED POWER LEVEL 


KNOWLEDGE BASE 

• ENGINE LDAO KNOWLEDGE - MEAN, CV, DISTRIBUTION TYPE 

■ INFLUENCE MODEL - INFLUENCE COEFFICIENT SET. GAIN INFORMATION, DITTY-CYCLE DATA 


Figure 3. 


COMPOSITE LOAD SPECTRA 
MULTI-LEVEL ENGINE MODEL 

COMPONENT LOAD MODELS 

• COMPONENT STATIC PRESSURE SCALING MODEL 

• TURBINE BLADE DYNAMIC PRESSURE LOAD MODEL 

• TRANSFER DUCT FLUCTUATION PRESSURE MODEL 

• PROBABILISTIC THERMAL LOAD MODEL 
KNOWLEDGE BASE 

• COMPONENT LOAD MODELS 

• DEPENDENCY ON ENGINE LOADS 

• LOAD PARAMETER DATABASE 

• INFLUENCE AND SCALING COEFFICIENTS 

• COMPONENT GEOMETRY DATABASE 

Figure 4. 
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SUMMARY HA7RU OF INDIVIDUAL LOADS 


VS COMPONENTS 


Individual 

load 


TURBINE 

LOAD 


TRANSFER 

DUCT 


r 

i 

j 

• 

STATIC PRESSURE | 

X 

1 

X 

I 

| X 

1 

1 

X 

i DUTY CYCLE 

1 

1 

i 

I 

* 

DYNAMIC PRESSURE I 


1 

1 


1 

I 

1 

1 



1 

I 

i 

i 


• CHUGGING (TRANSIENT) J 

- 

1 

1 

X 

I 

) 

I 

- 

1 AMS, STATOS 

1 

1 

i 

i 


• TURBULENCE 1 


1 

1 


1 

1 

1 

I 



1 

1 

i 


• SINUSOIDAL 1 


1 

1 


I 

1 

1 

1 



I 

[ 

i 

i 


(REPEATED PULSE) J 

X 

1 

j 


1 

1 

1 

1 


! AMS, PSD. STATOS 

1 

i 

i 


• RANDOM | 

- 

1 

1 

X 

1 

| X 

1 

1 

X 

I AMS. PSDS 

I 

1 

i 

i 

i 

• 

CENTRIFUGAL | 

X 

1 

1 

1 

- 

1 

1 

1 

1 

1 

- 

| DUTY CYCLE 

1 

1 

i 

i 

i 

• 

TEMPERATURE j 

X 

1 

X 

1 

[ X 

1 

1 

X 

1 DUTY CYCLE 

1 

I 

i 

i 

i 

• 

STRUCTURAL VIBRATION j 


1 

1 


[ 

I 

1 

1 



1 

I 

i 

i 

■ 


• transient 1 

- 

1 

1 

1 

X 

I 

I X 
j 

1 

1 

1 

X 

I AMS. STATOS 

1 

1 

i 

i 

i 

j_ 


* STEADY STATE j 

- 

1 

1 

j 

X 

l 

| X 
1 

j 

1 

I 

1 

- J- 

X 

1 AMS, PSD. STATOS 

1 

1 

l 


LOX 

POST 


HPOTPDD 

FORMAT 


LOAD FORM/ 
TEST DATA 
FORHAT 


Figure 5. 


COMPOSITE LOAD SPECTRA 
ANLOAD = THE LOAD CALCULATION MODULE 

• PROBABILISTIC METHODS 

• GAUSSIAN MOMENT METHOD 

• RASCAL (A VARIANCE OF DPD METHOD) 

• MONTE CARLO 

• PROBABILISTIC MODEL 

• STEADY STATE 

• TRANSIENT STATE 

• QUASI “STEADY STATE 

• POISSON ARRIVAL 

• RARE EVENT 

• PERIODIC MODEL AND PSD 

Figure 6. 
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Pr««aur» (p»la) 

Ctfnf1d«nc« Interval (ff-lb) (Thou*and*) 
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Conf1d«nc« interval (ft-lb) 


Initial Interval Effect On Confidence 

2.6 

2.4 
TJX 

2 

1*8 

1.6 

1.4 

1.2 

1 

0.8 
0.6 
0.4 
CL2 
0 

4.678 4.682 4.686 4.69 4.694 4.698 4.702 4.706 4.71 

(Thousands) 

Tonqus (ft— lb) 

□ Unoquab 25 4 Equal: 25 O Monts Carlo 


HPFTP Turbins Tarqus 



CONFIDENCE BAND CALCULATIONS 

0 HONTE CARLO BANDS ARE CALCULATED USING CLASSICAL STATISTICAL RESULTS 

0 MONTE CARLO CONFIDENCE BAND INTERVALS BECOME SMALLER {I.E. ACCURACY 
INCREASES) WITH INCREASING NUMBER OF SLHJUTICKS 

O RASCAL CALCULATIONS INCLUDE THE OCNTTNUOUS PORTION OF THE PROBABILITY 
SPACE IN THE CONFIDENCE BAND CALCULATIONS 

O THE CCNT LNTXXB PROBABILITY SPACE THAT IS NOT REPRESENTED IN THE 
INDIVIDUAL DISCRETE DISTRIBUTIONS CAUSES A DECREASE IN THE POINT 
ACCURACY BY RASCAL 

0 THEREFORE, THE CONFIDENCE BAND INTERVALS BECOME SMALLER WITH AN 
INCREASING NUMBER OF INITIAL DISCRETE INTERVALS 


Figure 9. 



LDEXPT LOAD EXPERT SYSTEM 


• EXPERT SYSTEM DRIVER 

• DECISION TREE INFERENCE 

• QUERY ON THE DATABASE KEY VARIABLES 

• LOAD DATABASE SYSTEM 

• STAND-ALONE DATABASE SYSTEM 

• EXPERT SYSTEM INTERFACE 

• KEY VARIABLES ARE ATTRIBUTES OF THE EXPERT SYSTEM 

• USER/EXPERT SYSTEM SELECT OPTIONS ON KEY VALUES 

• SIMPLE WORKING MEMORY MODEL 

• PASSING INFORMATION BETWEEN RULE MODULES 

• LDEXPT RULE MODULES 

• IMPLEMENTING PROCESS AND CONTROL KNOWLEDGE 

• e.g., RETRIEVING LOAD INFORMATION 

• IMPLEMENTING PROBLEM-SOLVING KNOWLEDGE 

• e.g.. SELECTING INDEPENDENT LOADS BASED ON GAINS 

• PROBABILISTIC LOAD MODULE - ANLOAD (BATTELLE) 


Figure 10. 


LDEXPT: LOAD EXPERT SYSTEM 



Figure 11. 
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PROBABILISTIC DESIGN ANALYSIS USING COMPOSITE LOADS SPECTRA (CLS) 

COUPLED WITH PROBABILISTIC STRUCTURAL ANALYSIS METHODOLOGIES (PSAM)* 

J. F. NEWELL, K. R. RAJAGOPAL and H. HO J3 3$ / 

ROCKWELL INTERNATIONAL CORPORATION, ROCKETDYNE DIVISION 

This paper discusses the application of the CLS to generate probabilistic 
loads for use in the PSAM nonlinear evaluation of stochastic structures under 
stress (NESSUS) finite element code. The CLS approach allows for quantifying 
loads as mean values and distributions around a central value rather than 
maximum or enveloped values typically used in deterministic analysis. NESSUS 
uses these loads to determine mean and perturbation responses. These results 
are probabilistically evaluated with the distributional information from CLS 
using a fast probabilistic integration (FPI) technique to define response 
distributions. The main example discussed describes a method of obtaining 
probabilistic (dependent) pressure, temperature and centrifugal steady state 
load descriptions and stress response of the second-stage turbine blade of the 
Space Shuttle main engine (SSME) high-pressure fuel turbopump ( HPFTP) . 
Additional information is presented on the on-going analysis of the high 
pressure oxidizer turbopump discharge duct (HPOTP) where probabilistic dynamic 
loads have been generated and are in the process of being used for dynamic 
analysis. Example comparisons of load analysis and engine data are furnished 
for partial verification and/or justification for the methodology. 

Figure 1 depicts the component loading analysis. The description of the 
dependent random variables are obtained from a multilevel physical model that 
includes an engine system influence model and component load models of the 
turbine to describe the local conditions on the turbine blade. The 
independent loads are quantified at the appropriate model - engine or 
component - to best define the basic variables in the system, as well as 
maintain the required correlation between the dependent loads on the turbine 
blade. Both duty cycle and related probabilistic variations are accounted for 
in these models. The specific turbine blade model. Figure 2, is used both for 
part of the load modeling-thermal and pressure distribution - as well as the 
structural analysis. The structural model combines the loads with a 

probabilistic description of geometry, material property, and material 
orientations. The work herein discussed is an extension of the analysis 
performed under the PSAM contract where probabilistic loads were not 
considered. 

Figure 3 is a cross section of the HPFTP turbine showing the overall 
configuration and the hot gas coolant network. The hot gas flow enters from 
the fuel preburner at 11 and discharges at 15. The second stage blade and 
rotor are between 14 and 15. Two geometric variations play a key role in the 
distribution of the coolant flow around the turbine blade, the interstage seal 
at 10 and the aft platform seal at 8. 

The HPFTP turbine blade loads are dependent to the engine system level major 
component (turbine) dependent loads. The major component dependent loads are 
in turn dependent to a set of engine independent loads, including the engine 
operating parameter loads and the engine hardware parameter loads. A summary 
and categorization of these loads are listed in figure 4. The independent 
system level loads and hardware variations are in the first category listed. 
The major component dependent loads that are calculated in the engine model 
are the next level. The component level dependent loads are pressure, 
temperature and centrifugal. Additional local independent loads from the 
seals are internal to the turbine. 

*Work performed under NASA Contract NAS3-24382 
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Probabilistic methods are built into these models and are described in the CLS 
discussion in this document. For this analysis, the Gaussian moment method 
was employed to illustrate the application of the probabilistic methodology. 
The pressure, temperature and centrifugal loads are correlated with higher 
level system independent loads. Theoretically, one could obtain correlation 
fields between pairs of the loads and then decompose them as part of the PSAM 
solution scheme. These fields are difficult to obtain and decompose. The 
marginal distribution method is used as a better solution to the load 
description and interface to the PSAM codes. Figure 5. Figure 6 outlines the 
process of generating the HPFTP turbine component loads marginal distribution. 

At this time, validation of these loads is limited to the standard 
instrumentation available from engine tests. Pump speed and turbine discharge 
temperatures are regularly recorded as well as the independent loads at the 
pump interfaces. Figure 7 has typical plots of test data from a series of 
tests that can be used to estimate 2 sigma variations of independent loads. 
Figure 8 has typical dependent pump speed variations. A sample from a flight 
is also furnished of the 3 pump speeds at 104% power level. The variation in 
speeds are caused by a combination of hardware and duty cycle variations of 
primarily the low pressure oxidizer pump inlet pressure and the resulting 
cavitation of the HPOTP. Table 1 furnishes a correlation of calculated and 
measured variations of turbine speed and discharge temperature for both high 
pressure turbopumps. There is excellent correlation with speed and slightly 
poorer correlation with temperature. Part of the temperature error is 

attributed to instrumentation error that is not readily accounted for. 

Table 2 summarizes the set of random variables, how they affect the structural 
model and their mean and standard deviation values. Figure 9 shows two points 
on the structural model where cumulative distributions of effective stress 
were calculated using NESSUS and FP1 codes. 

PSAM is performing a series of verification analysis on SSME hardware. CLS is 
being used to develop the loads for these analyses to furnish realistic input 
to PSAM and to demonstrate the integrated use of CLS and PSAM. Vibration 
loads for the HPOTP discharge duct. Figure 10, have been completed. This duct 
is attached to the HPOTP at two points. Zone G, and at the main injector (MI) 
at Zone A of the engine. The MI is part of a combustor that generates random 
vibrations. The pumps generate both random and sinusoidal loads. These 
vibration loads are dependent on power level of the engine and pump power and 
speed. The structure of the engine couples the responses of the components 
such that sinusoidals from the high pressure pumps are found throughout the 
engine at attenuated levels compared to the point where they are generated. 
Table 3 summarizes the 38 individual load components that are being applied to 
the duct. Figure 11 shows one of the six normalized PSD shapes used for the 
random loading. For this analysis, the random load composite has been modeled 
as a normalized shape that is scaled by G£ ms . The sinusoids are depicted 

in Figure 12. The frequencies generated by a pump are correlated to the 
synchronous speed, IN. The IN, 2N and 3N response levels are typically caused 
by rotor unbalance and rubbing or nonlinear effects. The 4N is caused by 
fluid loads on the four primary blades of the inducer. The variations are 
based on engine test data. The sinusoid amplitude mean levels are correlated 
with pump speed, but the standard deviation magnitudes are uncorrelated. The 
sinusoid level distribution is modeled as lognormal, see Figure 13. The 
variations of the sinusoidal magnitudes have a coefficient of variation of .5 
to over 1.0. The high variations are roughly equally caused by component 
hardware variations and test to test variations. Figure 14 is a recent study 
of random composite levels coefficients of variations for several environments 
on the SSME. This information will be used in the development of generic 
dynamic physical load models. 
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FIGURE 1. COMPONENT LOADS 



FIGURE 2. COMPONENT MODELS FOR LOADS AND STRUCTURAL ANALYSIS 
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FIGURE 3. SSME HPFTP TURBINE HOT GAS AND COOLANT FLOW NETWORK 
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FIGURE 4. COMPOSITE LOADS OF THE HPFTP TURBINE BLADE 
FOR THE PROBABILISTIC STRUCTURAL ANALYSIS 
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PURPOSE: 


• PREPARING COMPOSITE LOADS FOR STRUCTURAL ANALYSIS. 

• MARGINAL DISTRIBUTION LOADS USED DIRECTLY AS INPUT TO PROBABILISTIC STRUCTURAL 
ANALYSIS CODES, E.G. NESSUS & FPI . 

ADVANTAGES: 

• MUCH SIMPLER THAN CORRELATION FIELD APPROACH. 

• CORRELATIONS TO THE ENGINE INDEPENDENT LOAD LEVEL AS REQUIRED. 

• DIRECTLY USABLE FOR PERTURBATION LOADING IN PSAM. 

DEFINITION: 

• MARGINAL DISTRIBUTION FOR A COMPONENT LOAD IS A DISTRIBUTION FOR THE COMPONENT 
LOADS WITH ONLY ONE INDEPENDENT LOAD VARYING AND THE REST OF THE INDEPENDENT 
LOADS STAYING CONSTANT. 

• FOR A LINEAR SYSTEM 


AY = XB jAX j 

MARGINAL 

° Y, = 1 0 1 1 ° X , 

MARGINAL 

*Y_ . B IL 

3Xj B I I B , I * %i 

TIGURE 5. MARGINAL DISTRIBUTION METHOD 
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TOTAL RANDOM 
LOW NPSP - DET. 
HIGH NPSP - DET. 

360 

620 

-317 

— 

165 

225 

-219 

RANGE - RANDOM + DET. 

1660 

1500 

475 

MAX 

29090 

29000 

1630 

NOM 

28100 

— 

1374 

MIN 

27430 

27500 

1155 


— 

388 

396 

:::: 

65 

20 

— — 

— 

554 

56 

-94 

— 

70 

52 

-62 

:::: 



— - 

400 

1260 

1000 

114 

150 

1650 

35742 

35750 

1740 

1760 

— 

35130 

— 

1688 

— 

1250 

34482 

34750 

1625 

1610 


MEASURED — 
HARDWARE — 
TEST 
DET. 
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TABLE I. CALCULATED VS. MEASURED VARIATIONS IN PUMP SPEED AND TEMPERATURE 



NO 

RANDOM VARIABLE 
DESCRIPTION 

i 

JL 

TYPE 

1 FEH | 
(QUANTITIES ( 
1 AFFECTEO 1 

MEAN 

STANDARD 

DEVIATION 



1 

MATERIAL AXIS 

1 



-0.087?66 

0.067544 




ABOUT Z 

1 


I MATERIAL | 

RADIANS 

RADIANS 





1 

MATERIAL 






2 

MATERIAL AXIS 

1 


IORIENTATIONI 

-0.034907 

0.067544 




ABOUT Y 

1 

AXIS 


RADIANS 

RADIANS 





1 


| ANGLES I 





3 

MATERIAL AXIS 

1 

VARIATIONS 


+0.052360 

0.067544 




ABOUT X 

JL 


1 L 

RADIANS 

RADIANS 



4 

ELASTIC MODULUS 

1 



18.38E6 

0.4595E6 





1 

ELASTIC 


KSI 

KSI 





1 


I ELASTIC | 





5 

POISSON'S RATIO 

1 

PROPERTY 


0.386 

0.00965 





1 


| CONSTANTS | 





6 

SHEAR MODULUS 

1 

VARIATIONS 


18.63E6 

0 . 46575E6 





1 

j i 

KSI 

KSI 



7 

GEOMETRIC LEAN 

1 



0.0 

0.14 




ANGLE ABOUT X 

1 




DEGREES 





1 

GEOMETRY 

1 NODAL I 





8 

GEOMETRIC TILT 

1 



0.0 

0.14 




ANGLE ABOUT Y 

1 

1 

VARIATIONS 

(COORDINATES! 


DEGREES 



9 

GEOMETRIC TWIST 

1 

1 



0.0 

0.30 




ANGLE ABOUT Z 

JL 




DEGREES 



10 

MIXTURE RATIO 

1 








LIQUID HYDROGEN/ 

1 

INDEPENDENT 


6.00 

0.02 




LIQUID OXYGEN 

1 









1 

LOAD 






n 

FUEL INLET 

1 



30.00 psia 

5.00 




PRESSURE 

1 

i 

VARIATIONS 

| PRESSURE | 





12 

OXIDIZER INLET 

I 

1 

DEPENDENT 


100.00 psia 

26.00 




PRESSURE 

1 


(TEMPERATURE) 







1 

LOADS ARE 






13 

FUEL INLET ^ 

1 



38.5* R 

0.5 




TEMPERATURE 

(TURBINE BLADE 

i 

ICENTRIFUGALI 





14 

OXIDIZER INLET 

i 

i 

PRESSURE 


167.0* R 

1.33 




TEMPERATURE 

i 


| LOAD | 







i 

TEMPERATURE 






15 

HPFP EFFICIENCY 

i 



1.00 

0.008 





i 

AND SPEED 






16 

HPFP HEAD 

i 



1.0237 

0.008 




COEFFICIENT 

JL 







17 

COOLANT SEAL 

1 

LOCAL 


1 .00 

0.1 




LEAKAGE FACTOR 

1 









1 

GEOMETRY 

(TEMPERATURE) 





18 

HOT GAS SEAL 

1 



1.0 

0.5 




LEAKAGE FACTOR 

L 

FACTORS 


-] 




TABLE 2. LIST OF 2nd STAGE BLADE RANDOM VARIABLES 
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Location A 



Mean = 61,086 (PSI), SIGMA = 1,280 (PSI) 



(Thousands) 

Effective Stress (PSI) D5051-7 



(Thousands) 

Effective Sires: (PSI) D5051-6 


FIGURE 9. STRESS RESPONSE FOR HPFTP TURBINE BLADE 
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HPOTP 
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MAIN 

INJECTOR 


X 


X 


N 


FIGURE 10. HPOTP DISCHARGE DUCT 
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SEQUENTIAL 

NO. 

1 

_L 

RANDOM VARIABLE 
DESCRIPTION 


| MEAN | 

I 1 

ST. D. 

1 

1 

cov 

1) 

1 

ZONE G - X AXIS, 

PSD 

| 222.0 | 

163.0 

| 

0.73 

2) 

1 

ZONE G - Y AXIS, 

PSD 

1 73.5 | 

59.4 

| 

0.808 

3) 

1 

ZONE G - Z AXIS, 

PSD 

1 73.5 | 

59.4 

| 

0.808 

4) 

1 

ZONE A - X AXIS, 

PSD 

1 22.5 | 

4.5 

| 

0.20 

5) 

1 

ZONE A - Y AXIS, 

PSD 

1 54.0 | 

10.8 

| 

0.20 

6) 

1 

ZONE A - Z AXIS, 

PSD 

I 69.5 | 

13.9 

| 

0.2 

7) 

1 

OXIDIZER PUMP SPEED 

I 2940.53 | 

41 .167 

| 

0.014 

8) 

1 

FUEL PUMP SPEED 


I 3707.08 | 

37.07 

| 

0.01 

9) 

1 

1 

1 

DAMPING 

ZONE A OXIDIZER TURBOPUMP 
OXIDIZER PUMP SINE AMPLITUDES 

1 0.033 | 

1 1 

1 | 

.005 

1 

1 

1 

0.15 

10) 

1 

X DIRECTION 

IN 

1 .30 | 

.120 

| 

0.4 

11) 

1 


2N 

I 0.30 | 

.045 

| 

0.15 

12) 

1 


4N 

1 1.5 | 

.45 

| 

0.3 

13) 

1 

Y DIRECTION 

IN 

I 0.60 | 

0.30 

| 

0.5 

14) 

1 


2N 

I 0.70 | 

0.28 

| 

.40 

15) 

1 


4N 

1 2.6 | 

.78 

| 

.3 

16) 

1 

Z DIRECTION 

IN 

1 0.5 | 

.225 

| 

0.45 

17) 

1 


2N 

I 0.70 | 

.140 

| 

0.20 

18) 

1 

1 

4N 

FUEL PUMP SINE AMPLITUDES 

I 0.70 | 

1 | 

.140 

1 

1 

.20 

19) 

1 

X DIRECTION 

IN 

1 .35 | 

1 .050 

1 

| 

.3 

20) 

1 

Y DIRECTION 

IN 

I 0.80 | 

0.280 

| 

0.35 

21) 

1 

1 

1 

1 

Z DIRECTION IN 

ZONE G - MAIN INJECTOR 
OXIDIZER TURBOPUMP 
SINE AMPLITUDES 

1 1.20 | 
I 1 

1 1 

I | 

0.96 

1 

1 

1 

1 

0.3 

22) 

1 

X DIRECTION 

IN 

I 1.35 | 

1.35 

i 

| 

1 .0 

23) 

1 


2N 

I 1.50 | 

0.75 

| 

0.5 

24) 

1 


3N 

1 1.10 | 

0.495 

| 

0.45 

25) 

1 


4N 

1 n.o | 

2.75 

| 

0.25 

26, 27) 

1 

X-Y DIRECTION 

IN 

1 1.9 | 

1 .71 

1 

0.9 

28, 29) 

1 


2N 

1 1.6 | 

0.96 

1 

| 

0.6 

30, 31) 

1 


3N 

| 0.75 | 

.225 

| 

0.3 

32, 33) 

1 

1 

4N 

FUEL PUMP SINE AMPLITUDES 

1 5.5 | 

I | 

3.30 

1 

1 

0.6 

34) 

1 

X DIRECTION 

IN 

| 0.65 | 

0.2275 

V 

1 

0.35 

35, 36) 

1 

X-Y DIRECTION 

IN 

I 0.45 | 

1.35 

i 

.3 

37. 38) 

_L_ 

Y-Z DIRECTION 

2N 

1 0.45 1 

0.180 

1 

.4 


Note: 1) P.S.O. units are in G^ 

2) Pump Speed units are in Radians/SCE 

3) Sinusoidal amplitude units are in GRMS 


TABLE 3. BASIC RANDOM VARIABLES AND THEIR STATISTICS (INPUT) 
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BHS POWER SPECTRAL DENSITY 


ZONE B. X-AXIS 



FIGURE 11. TYPICAL NORMALIZED PSD 
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FIGURE 12. RANDOM AND SINUSOIDAL VIBRATION 
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PROBABILISTIC FINITE ELEMENTS FOR FRACTURE AND FATIGUE ANALYSIS 

W.K. Liu, T. Belytschko, M. Lawrence, and G.H. Besterfield 
Department of Mechanical Engineering 
Northwestern University 
Evanston, Illinois 


The fusion of the probabilistic finite element method (PFEM) and reliability' analysis 
for probabilistic fracture mechanics (PFM) is presented. The authors have developed a 
comprehensive method for determining the probability of fatigue failure for curved crack 
growth. The criterion for failure or performance function is stated as: the fatigue life of a 
component must exceed the service life of the component; otherwise failure will occur. An 
enriched element that has the near crack-tip singular strain field embedded in the element is 
used to formulate the equilibrium equation and solve for the stress intensity factors at the 
crack-tip. The crack growth law (e.g., Paris law) and crack direction law (e.g., direction 
which gives rise to the maximum tangential stress at the crack-tip) are also incorporated into 
the formulation. The loading is mixed-mode (i.e., opening mode and in-plane shear mode 
which are designated mode I and mode II, respectively) with randomness in the applied 
loads, material properties, component geometry, crack geometry, crack growth law, and 
the crack direction law. Examples of random crack geometry are the initial and final crack 
lengths, initial crack angle, and initial crack position. The methodology consists of 
calculating the reliability index, which is used to calculate the first-order probability of 
failure, by solving a constrained optimization problem. Included in the optimization for 
reliability is the equilibrium equation for stress intensity factors, the crack growth law, and 
the crack direction law. 

Performance and accuracy of the method is demonstrated on a classical mode I 
fatigue problem. The initial and final crack lengths, material parameters in the Paris law, 
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and the applied stress are modeled with random variables. The results imply that the 
service life of a component must be well below the deterministic fatigue life of the 
component when many parameters are considered random. The fusion of PFEM and 
reliability for fatigue crack growth is computationally quite efficient and provides a 
powerful tool for the design engineer. 
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SCOPE OF WORK 


Probabilistic Finite Element Method fPFEMl 

■ Uncertainties 

Random Load 

Random Material Properties 
Random Geometry 

• Statics and Dynamics 

• Elasticity and Nonlinear Elasticity 


Brittle Fracture bv PFEM 

• Probability of the Stress Intensity Factor exceeding 
the Fracture Toughness of the Material 


P ( K I > Kj c ) 

P(C>G c ) 


Fatigue Crack Growth bv PFF.M 

• Probability of the Desired Service Life exceeding the 
Fatigue Life 


P ( T s > T ) 


Introduction 


Probabilistic Fracture Mechanics fPFM't 

• Model the crack-tip singularity in stress. 

• Determine statistics of stress intensity factors. 

• Probability of fracture and fatigue failure. 


Uncertainties which need to be considered in PFM are: 

• Random load, material properties, and geometry; 

• Random crack geometry; 

length, 

location, 

orientation, 

• Random crack growth laws; 

• Random crack direction laws; 

• Random correlation among detectable cracks and 
micro-cracks. 
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Reliability of Components 
Subject to 

Fatigue Crack Growth 
by 

Probabilistic Finite Elements 


Methodology and Objective 

• Given : 


• Component with Crack Subject to Cyclic Stress 

* Method for Determining Stress Intensity Factors 

• Crack Growth Law 

* Crack Direction Law 

• With Randomness in : 

• Applied Loads 

■ Component Geometry 

• Crack Geometry 

* Crack Direction Parameters 

* Crack Growth Parameters 

• Material Properties 

• Determine Probability of Fatigue Failure 

* Probability that the desired service life is greater 
than the fatigue life. 

P(T s >T) 


Stress Intensity Factor 



Kj = mode I stress intensity factor 
K II = mo ^ e H stress intensity factor 
K m = mode III stress intensity factor 


Consider only Modes I and II: 


x 


K I 

K II _ 
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Random Variables 


System Uncertainties 

• Applied Loads 

• Component Geometry 

length, height, width, etc. 

• Crack Geometry' 

length, orientation, position 

• Crack Direction Parameters 

• Crack Growth Parameters 

• Material Properties 

elastic modulus, poisson's ratio, 
fracture toughness, etc., 


Modeled by a correlated random variable vector 
b (qxl) 


Transform b to a Independent Standard Normal Random 
Variable Vector 

r (qxl) 


* Eigenvalue Orthogonalization (Normal R.V.’s) 

Covariance Matrix for b 

* Rosenblatt Transformation (Non-normal R.V/s) 

Joint Probability Distribution for b 


Reliability Analysis 

Performance Function or Margin of Safety 

g > 0 g = 0 g<0 

safe state limit- state failure state 

Probability of Failure 
P f = P(g < 0) 

First-Order Estimate of the Probability of Failure 
P fl - (-P) 

<J> is the standard normal cumulative distribution, 

(3 is the reliability index. 

Reliability Index 

Minimum Distance from Limit-State Surface to 
Origin in Independent Standard Normal Space 

P 2 = Vr T r 
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Log (Ak) 


Region A: A sufficient level of activity at the crack-tip has 
not been attained resulting in very slow crack 
growth. 

Region C: Corresponds to rapid unstable crack growth. 

Region B: Represents the zone where most fatigue crack 
growth is analyzed. 


Paris Law: ^ = D (AK eq ) n 

D is a material parameter 
n is a material parameter 
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Integrating from the initial crack length to the final crack 
length gives in the fatigue life 


/■ da 

. D(A* eq ) n 

a i 


Equivalent Stress Intensity Factor 
* Maximum Stress Theory 


eq 


© T x 


e 



COS 5- 

, . e 

3 sin j 


CRACK DIRECTION LAW 


• Maximum Tangential Stress Criterion 

• Minimum Strain Energy Density Criterion 

• Maximum Energy Release Rate Criterion 


Maximum Tangential Stress Criterion 

• Crack grows in a radial direction from the crack 

* Crack grows in the plane perpendicular to the 
direction of maximum tension 


Crack grows in the direction where the tangential 
stress is a maximum and the shear stress is zero. 

Z(k, 6) s O t r = 0 

f sin 0 I 


O = 


L 3 cos 0 - 1 J 

For mode I fracture, 9 = 0° and is constant 
For mode II is fracture, 0 = 70.54° and changes 
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PERFORMANCE FUNCTION FOR FATIGUE 


Basis for determining the reliability and/or probability of 
failure depends on the choice of a performance function. 

p f = P(g < 0) 


Performance function for fatigue crack growth 


T s is the desired service life 
T is the fatigue life 

The component fails when the fatigue life is less than 
the desired service life. 


DISCRETIZATION 


Discretize the Crack Path into "nint” Integration Points 

* Crack Length 

a i = \ [(a f - a t ) + (a f + ap] i = 1 nint 

• Crack Growth Law 


T = 


+1 

f 

- 1 


D ( K eq) 


¥ 


J 


da 


• Crack Direction Law 

T 

Z k = ^k x k = 0 k = 1 nint 

Crack direction is calculated at each integration point 
Crack grows to the next integration point 

• Equilibrium 

Kj = fj i = 1, nint 
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RELIABILITY AS AN OPTIMIZATION PROBLEM 
FOR FATIGUE CRACK GROWTH 


First-Order Probability of Failure 
— > Reliability Index 

— » Constrained Optimization Problem 

Minimize the Reliability Index 

P 2 = VT7 

Subject to Equality Constraints 
Crack growth direction law 

2^ = 0 k = I, nint 

Equilibrium 

Ki ^i = Fi i = l, nint 

Subject to Inequality Constraint 

Limit-State Surface and Failure State 
T - T s < 0 
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Define a Lagrange Functional 


& (b, \l | , H n i nP $ 1 . ^nint’ a > 

‘Pl ( Pnint’ 0 l e nim> 

nint 

= (3 2 + X |i T [ ^ . K j ] 

1=1 

nint 

+ X <Pk t z k = 0 ] 

k=l 

+ M T - T s + a2 } 


Necessary Conditions 

(derivatives w.r.t. independent variables) 


Probabilistic Finite Element Method 
(PFEM) 


Eliminate Lagrange Multipliers 


Final Form of Optimization for Reliability 


+ \ L g = 0 


L P 



L.S 



X 

k=l 

nint 


- X 

k=l 


( T, 


6i- 


(z k ), 

V T >c 


(Zk)' 


(Z k )' Kj 

-(z k )'b 


+ T 


+ T.v 


db 


Solve for b via an Optimization Algorithm 

* Lagrange Multipliers 

* Gradient Projection Algorithm 

* HL-RF Algorithm (Tangent Matching) 

* Modified HL-RF Algorithm 
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CRACK GROWTH MOVING MESH 


• Variation of ALE formulation 

Pseudo-Static (neglect transport term) 

• Moving Element Mesh with Constraints 

Aspect ratio of elements must remain one 
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EXAMPLE 


Single-Edged Cracked Plate with an Applied Load 
b = [ a 1 , a f , D , n , a ] T 

a 3 is the initial crack length 
a f is the final crack length 

D is a material parameter 
n is a material parameter 
a is the applied stress 



Mean 

WKSSSSsBEESBSM 

Percent 

Lengtn (Lj 

10.0 in 

0.0 

0.0 

Width (W) 

4.0 in 

0.0 

0.0 

Thickness (t) 

1.0 in 

0.0 

0.0 

Young's Modulus (E) 

30,000.0 ksi 

0.0 

0.0 

Poisson's Ratio (v) 

0.3 

0.0 

0.0 

Applied Stress (o) 

12.0 ksi 

3.0 ksi 

25.0 

Initial Crack Length (ap 

0.01 in 

0.01 in 

100.0 

Final crack Length (ap 

0.1 in 

0.01 in 

10.0 

Fatigue Parameter (D) 

1.0 x 10' 10 

3.0 x 10' 11 

30.0 

Fatigue Parameter (n) 

3.25 

0.08 

2.5 


• Variations in the parameters are reasonable. 


• nint = 50 


• Deterministic fatigue life, T = 7,500,000 cycles 
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R I R * R i R 1 



Reference Solution 
Kj = ^ a Va 

^ = 1-99 - 0 . 41 ^ + 18 . 7 (|r ) 2 - 38 . 48(^) 3 + 53 . 85(^) 4 


Deterministic Stress Intensity Factor for Reference and FEM Solutions 
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Reliability Index Dclcrminislic Fatigue Life (cycles) 









Probability of Failure for Combined Randomness 
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PROBABILITY OF FAILURE AND RISK ASSESSMENT 
OF PROPULSION STRUCTURAL COMPONENTS 


^2 £-3 f 
/^3SS 


Michael C. Shiao 

NASA Lewis Research Center Group 
Sverdrup Technology, Inc. 
Cleveland, Ohio 


ABSTRACT 


The probabilistic structural analysis method (PSAM) has been developed to 
analyze the cfTects of fluctuating loads, variable material properties, and uncertain 
analytical models especially for high performance structures such as SSME 
turbopump blades. In the deterministic approach, uncertainties in the responses 
were not quantified and the actual safety margin remains unknown. Risk is calcu- 
lated after expensive service experience. However, probabilistic structural analysis 
provides a rational alternative method to quantify uncertainties in the structural 
performance and the durability. NESSUS (Numerical Evaluation of Stochastic 
Structures Under Stress) has been developed by SwRI under the PSAM research 
project sponsored by NASA Lewis. It is a probabilistic structural analysis computer 
code which integrates finite element methods and reliability algorithms, capable of 
predicting the probability distributions of structural response variables such as 
stress, displacement, natural frequencies, and buckling loads etc. In the design and 
analysis practice, it is also important to assess the risk associated with a new or 
existing structure for safety or serviceability consideration. Probable failure criteria 
should be identified for different structures due to various structural functions. For 
example, failure events such as stress greater than strength or the displacement ex- 
ceeds maximum allowable are often used for the reliability assessment. Probability 
of occurrence of those failure events can be determined once the probability dis- 
tributions of the stress or displacement are calculated by NESSUS. 

In undertaking a rcliability/risk analysis, all suspected sources in uncertainties 
should be taken into account in order to control the probability of failure in service 
environments within an acceptable range. Reliability and risk obtained by a prob- 
abilistic structural analysis can be useful in evaluating the traditional design, setting 
quality control requirements and inspection intervals. It can also be used to iden- 
tify candidate material and design concepts in the absence of technology base. In 
this presentation, a risk/cost assessment and a reliability analysis using NESSUS 
with a generic probabilistic material model is described. 

The risk assessment includes the initial cost and the cost due to failure. The 
initial cost is defined as the cost for component service readiness which can be a 
function of several key design variables. The consequential cost is the cost incurred 
due to failure. Total cost is the sum of initial cost and a fraction of consequential 
cost. The fraction is weighted by the probability of failure. Since the lower initial 
cost is often associated with higher risk for the structural failure, higher initial cost 
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will normally reduce the risk. IIow to minimize the total cost for an acceptable 
structural rcliability/risk is also described in this presentation. 

In previous studies, probability distributions of material properties such as 
Young's modulus, thermal expansion coefficient, and material strength were as- 
sumed. In the present study, a generic probabilistic material model is incorporated 
in NESSUS. The assumption for this model is that the probability distributions of 
material properties arc simulated by using primitive random variables. These vari- 
ables can be temperature, stress, fatigue cycles, etc. The model is especially at- 
tractive for the structures under high temperature environments. A risk analysis of 
an SSME blade is performed under this concept. Since the material properties are 
functions of stress, vet stress is a function of material properties in an implicit way, 
an iterative procedure is necessary to obtain a convergent solution for stress and 
material properties. During the iteration process, the joint cumulative distribution 
functions of nodal stresses are required in order to apply the generic probabilistic 
material model. It is found that only a few iterations are needed for a converged 
result. From this analysis, a Risk-Fatigue cycle curve is developed for critical lo- 
cations. This curve is useful for assessing the risk of existing structures. For in- 
stance, at a given acceptable risk level, the number of fatigue cycles to initial local 
failure can be determined. With this information available, criteria can be set for 
quality control, inspection intervals and retirement for cause. 

In summary, a rcliability/risk cost methodology has been developed. It consists 
of a probabilistic structural analysis by NESSUS and a generic probabilistic mate- 
rial model. The methodology is versatile and particularly applicable to high tem- 
perature structures where data is difficult to obtain. It is demonstrated by using it 
to assess the risk associated with fatigue cycles to initiate local failure. 
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OBJECTIVE 


DEVELOP A METHODOLOGY TO EVALUATE PROBABILITY OF FAILURE 
AND PERFORM A RISK ASSESSMENT USING NESSUS AND A GENERIC 
PROBABILISTIC MATERIAL MODEL OF COMPLEX STRUCTURES 
OPERATING IN HIGH TEMPERATURE SERVICE ENVIRONMENTS 


CD-89-39414 


UNCERTAINTIES IN THE PROBABILISTIC 
STRUCTURAL ANALYSIS 

STOCHASTIC 

THERMOMECHANICAL 
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PROBABILISTIC STRUCTURAL ANALYSIS 

BY NESSUS & A GENERIC PROBABILISTIC MATERIAL MODEL 
ITERATION 0 

STEP 1: CALCULATE THE PROBABILITY DISTRIBUTIONS OF 

MATERIAL PROPERTIES USING THE GENERIC 
MATERIAL MODEL WITHOUT STRESS AND FATIGUE 
CYCLES 

STEP 2: CALCULATE THE NODAL STRESS BY NESSUS 

ITERATION 1 

STEP 1 : CALCULATE THE PROBABILITY DISTRIBUTIONS OF 

MATERIAL PROPERTIES USING COMPLETE GENERIC 
MATERIAL MODEL 

STEP 2: CALCULATE THE NODAL STRESS BY NESSUS WITH 

UPDATED MATERIAL PROPERTIES 

ITERATION 2 

(REPEAT THE ITERATIONS UNTIL THE PROBABILITY 
DISTRIBUTIONS OF MATERIAL PROPERTIES AND 
STRESS HAVE CONVERGED) 


CD-89'39422 


JOINT CUMULATIVE DISTRIBUTION FUNCTION 
OF NODAL STRESSES 

Oj & Oj CALCULATED BY FIRST-ORDER SECOND-MOMENT METHOD 


o, = NODAL STRESS 
y, = REAL VALUE 
P/ = RELIABILITY INDEX 
<t> = NORMAL CDF 
<p = NORMAL PDF 
P ij = cos V,j 



NORMALIZED COORDINATE SYSTEM 


P (Oj < y, and o, < Yj ) = Of- pjOf- fy) + g«f <p (- - p j;Z )dz 

CD-89-39423 
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GENERIC PROBABILISTIC MATERIAL PROPERTY 
MODEL IN TERMS OF PRIMITIVE VARIABLES 


1 

1 

n 

Sp — o 

p 

log N mf - log N m 

1 

K° 

1 

1 


S F -a 0 


,0 9 n mf ~ ,0 9 n mo 


PRIMITIVE VARIABLES 
M p = MATERIAL PROPERTY 
T - TEMPERATURE 
S = STRENGTH 
a = STRESS 

N m = MECHANICAL CYCLES 


SUBSCRIPTS 

F = FINAL CHARACTERISTIC VALUE 
0 = REFERENCE PROPERTY 
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PRIMITIVE VARIABLE PROBABILITY DISTRIBUTIONS 
FOR PROBABILISTIC MATERIAL PROPERTY MODEL 


VARIABLE 

DISTRIBUTION 

MEAN 

STANDARD DEVIATION 


TYPE 


(VALUE) 

(% OF MEAN) 

Tf 

NORMAL 

2750 °F 

51.4 °F 

2.0 

T 0 

NORMAL 

68 °F 

2.04 °F 

3.0 

S F 

NORMAL 

212.0 ksi 

10.6 ksi 

5.0 

°0 

CONSTANT 

0 

0 

0 

n mf 

LOGNORMAL 

10 8 

5x10 6 

5.0 

n mo 

LOGNORMAL 

10 3 

50 

5.0 

n 

NORMAL 

0.25 


3.0 

P 

NORMAL 

0.25 


3.0 

q 

NORMAL 

0.25 


3.0 

CD-89-3^ 
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RISK ASSESSMENT 


fo'f S 


7 

\ 

\ 

\ 


u PROBABILITY 
OF FAILURE 


a > 

I 

I 

I 

I 

l stress 


s 

I 

I 

I 


STRENGTH 


(J?.'sfc)<k)Wdx 


r STRESS-a 
(FROM NESSUS) 


r STRENGTH-S 

\ (FROM GENERIC PROBABILISTIC 
\ MATERIAL PROPERTY MODEL) 
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PROBABILITY OF LOCAL FAILURE 
CAN BE QUANTIFIED 

(STRESS EXCEEDS STRENGTH) 



PROBABILITY OF FAILURE 


. 0869 . 

.0790 

. 0711 . 

. 0632 . 

. 0553 . 

. 0474 . 

.0395 

.0316 

.0237 

.0158 

.0079 


PROBABILITY 9 

CD-89-39419 
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HYPOTHESIS FOR QUANTIFYING 
ACCEPTANCE OF RISK 


Cj INITIAL COST = COST OF COMPONENT SERVICE 

READINESS 

C f CONSEQUENTIAL COST = COST DUE TO FAILURE 

OCCURRENCE 

PROBABILITY OF FAILURE 

TOTAL COST 

c t = Cj + PfCf 

CD-87-39424 


PROBABILITY OF FAILURE CAN BE QUANTIFIED 
IN TERMS OF IMPROVED MEAN STRENGTH 



STRESS/STRENGTH MEAN STRENGTH 

CD-87-39425 


Pf 

C t 
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THE INITIAL COST TO IMPROVE THE STRUCTURAL 
RELIABILITY CAN BE QUANTIFIED IN TERMS OF 
MEAN STRENGTH (GIVEN QUALITY) 


INITIAL 

COST 



PROBABILITY 

OF 

FAILURE 


CD-8?- 39426 


THE TOTAL COST TO IMPROVE THE STRUCTURAL 
RELIABILITY CAN BE QUANTIFIED IN TERMS OF 
MEAN STRENGTH (GIVEN QUALITY) 


TOTAL 

COST 



PROBABILITY 

OF 

FAILURE 


CD-89-39427 
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SSME BLADE SHOWING LOCATIONS WHERE 
PROBABILISTIC STRUCTURAL RESPONSE 
WAS EVALUATED 


LOCAL 

STRESSES 



GLOBAL 

DISPLACEMENTS 


CD-89-39420 


PROBABILISTIC MODULUS SIMULATED BY USING 
THE GENERIC PROBABILISTIC MATERIAL 
PROPERTY MODEL 


MEAN STRESS 
STRENGTH RATIO = 0.83 


MEAN STRESS 
STRENGTH RATIO = 0.82 


PROBABILITY 

OF 

OCCURRENCE 



CD-09-39429 
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PROBABILISTIC FATIGUE STRESS/STRENGTH 
SIMULATED BY USING THE GENERIC 
PROBABILSTIC MATERIAL PROPERTY MODEL 


MEAN STRESS 
STRENGTH RATIO = 0.83 


MEAN STRESS 
STRENGTH RATIO = 0.82 


PROBABILITY 

OF 

OCCURRENCE 
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PROBABILISTIC DISPLACEMENTS 
CALCULATED BY NESSUS 


PROBABILITY 

OF 

OCCURRENCE 



DISPLACEMENT AT DISPLACEMENT AT 

NODE C, In. NODE D, In. 
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PROBABILISTIC CYCLIC LOAD 
EFFECTS ON MODULUS 



CD-89-39432 


PROBABILISTIC CYCLIC LOAD 
EFFECTS ON STRENGTH 



CD-89-39433 
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PROBABILITY OF LOCAL FAILURE DUE 
TO FATIGUE CYCLES 


PROBABILITY 

OF 

FAILURE 



CD-89-39434 


CONCLUDING REMARKS 

• A METHODOLOGY HAS BEEN DEVELOPED FOR PERFORMING 
RELIABILITY/RISK COST BASED ASSESSMENTS OF AEROSPACE 
PROPULSION STRUCTURES 

• THE METHODOLOGY CONSISTS OF PROBABILISTIC STRUCTURAL 
ANALYSIS (NESSUS) AND A PROBABILISTIC GENERIC MATERIAL 
PROPERTY MODEL 

• THE METHODOLOGY IS ILLUSTRATED BY USING IT TO ASSESS 
THE RISK FOR QUALITY CONTROL, STRENGTH IMPROVEMENT, 
AND FATIGUE CYCLES TO INITIATE LOCAL FAILURE 


CD-89-39435 
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OVERVIEW OF AEROTHERMODYNAMIC LOADS DEFINITION STUDY 

53 7 - 3 ? 

/f#f¥ 


Raymond E. Gaugler 
NASA Lewis Research Center 


For a number of years now, NASA has been conducting the Advanced Earth- 
to-Orbit (AETO) Propulsion Technology Program, now part of the Civil Space 
Technology Initiative (CSTI) Program. The program objective is to provide the 
knowledge, understanding, and design methodology that will allow the develop- 
ment of advanced Earth- to-orbit propulsion systems with high performance, ex- 
tended service life, automated operations, and diagnostics for in-flight 
health monitoring. This program supports both current and future engine and 
vehicle developments. The technology needs of the current Space Shuttle Main 
Engine (SSME) have thus far provided the main focus for the AETO efforts. 
Organizationally, the program has been divided into technology working groups 
made up of people from both the Marshall Space Flight Center and the Lewis 
Research Center. The focus of this report is the Lewis effort under the Fluid 
and Gas Dynamics working group, the Aerothermodynamic Loads Definition Study. 

The objective of the Aerothermodynamic Loads Definition Study is to 
develop methods to more accurately predict the operating environment in AETO 
propulsion systems, such as the Space Shuttle Main Engine (SSME) powerhead. 
Development of time-averaged and time-dependent, three-dimensional viscous com- 
puter codes as well as experimental verification and engine diagnostic testing 
are considered to be essential in achieving that objective. Time-averaged, 
nonsteady, an- i transient operating loads must all be well defined in order to 
accurately predict powerhead life. Improvements in the structural durability 
of the SSME turbine drive systems will depend on our knowledge of the aerother- 
modynamic behavior of the flow through the preburner, turbine, turnaround 
duct, gas manifold, and injector post regions. 

The approach taken under this study consists of two parts: (1) to modi- 

fy, apply, and disseminate existing computational fluid dynamics (CFD) tools 
in response to current needs and (2) to develop new technology that will ena- 
ble more accurate computation of the time-averaged and unsteady aerothermody- 
namic loads in the SSME powerhead. The new technology development effort is 
further divided into three parts: (1) new computer code development, (2) ex- 

periments to provide data for physical models required by the codes, and (3) 
experiments to provide data for validating the codes. With the more accurate 
aerothermodynamic load predictions providing boundary conditions to improved 
structural and fatigue life analyses, the goal of improved durability will be 
met . 


The Aerothermodynamic Loads Definition Study was begun in October 1983. 
The initial effort involved the use of existing CFD tools to analyze problems 
in the fuel and oxidizer turbopump turbines, the fuel turbine turnaround duct, 
the fuel-side preburner, and the main injector liquid oxygen posts. Results 
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of those efforts have been presented previously and will not be discussed 
here. What follows is a description of the currently active tasks and their 
status . 


H4 ( a) . UNSTEADY VISCOUS MULTI BLADE -ROW TURBINE ANALYSIS 

In this task, the objective is to develop a numerical simulation capabili- 
ty for unsteady viscous flow and heat transfer in multiblade-row turbines. 

The numerical simulation of unsteady viscous flows in multiblade-row tur- 
bines requires a means of reducing the three-dimensional, nonperiodic, multi- 
row problem to one of more manageable proportions. This work will develop a 
two-dimensional test bed to be used in exploration for ways of solving the 
full problem. The approach is to use an existing quasi-3D viscous cascade 
code as the basis of a 2D multipassage rotor-stator interaction code. 

A detailed description of the latest efforts under this task will be pre- 
sented later in this session. 


H4(b) . COMPUTATION OF THREE-DIMENSIONAL FLOW WITH HEAT TRANSFER 

In this task the objective is to develop an accurate, efficient numerical 
analysis of steady 3D viscous flow and heat transfer in turbine blade rows. 

The approach being taken is to modify a 3D viscous turbomachinery code 
( RVC3D) to include an advanced turbulence model and increased resolution of 
near-wall gradients and to validate and upgrade the code by application to 
suitable 3D data sets. 


H4(e). THREE-DIMENSIONAL STATOR-ROTOR BLADE FLOW EXPERIMENT 

The objective of this task is to investigate the steady and unsteady 3D 
viscous flow field in an axial turbine stage, including stator-rotor interac- 
tion, mixing effects, blade and end wall boundary layers, and temperature gra- 
dient effects on secondary flow development. This work is being done in the 
Pennsylvania State University Axial Flow Turbine Research Facility. A 
progress report will be presented later in this session. 


H4 (m) . TURBULENCE MODELING FOR PROTEUS TURBOMACHINERY CODE 

The objective of this task is to develop the capability to predict turbu- 
lent flow losses and heat transfer in the SSME turbines using the Proteus 3D 
Navier-Stokes code. Existing turbulence models are inadequate or require em- 
pirical tuning to generate reasonable results. Efforts have begun, working 
first with a 2D version of the code. 


H5( a) . IMPROVED MODELING, PREBURNER FLOW 

In this task the objective is to obtain steady-state reacting hydrogen 
data of mean and fluctuating velocities and concentrations and to compare 
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benchmark data with predictions from computer models and evaluate turbulence- 
chemistry interactions. 

Computations of the 3D flow in the SSME turbopump turbines requires knowl- 
edge of the condition of the gases exiting the preburner. Current models for 
predicting those conditions are inadequate. Fluctuating concentration and ve- 
locity data can now be obtained with laser diagnostics. These measurements 
are necessary to determine the unmixedness of the reacting flow and the reac- 
tion rates. The measurements are also necessary for differentiating the vari- 
ous computer models. 

For this task a team of computer modelers and experimentalists have de- 
fined the measurements required to verify turbulent-reacting flow models. 
Existing optical techniques will be used to make measurements within the shear 
layer: two-component laser Doppler velocimetry (LDV), fluorescence for OH and 

temperature, and TiCl 4 for seeding and H 2 0 concentration. 

Hardware fabrication is nearly complete. Figure 1 is a schematic view of 
the test rig, and figure 2 is a photograph of the test section, mounted in an 
inspection facility. 


H5(f). HIGH REYNOLDS NUMBER AND UNSTEADY HEAT TRANSFER EXPERIMENT 

In this task the objective is to develop a fundamental understanding of 
the effects of high Reynolds numbers and Reynolds number variations, and also 
upstream unsteadiness, on heat transfer in turbulent boundary layers typical 
of turbine airfoils in the SSME turbopump turbines. 

In order to predict the surface pressure and temperature distributions 
within engine flow passages, an accurate turbulence and unsteadiness model is 
needed. Data in the proper range of Reynolds number and disturbance frequency 
are necessary for accurate modeling. A new in-house experimental rig, shown 
schematically in figure 3, is being built in which the major features of an 
SSME turbine airfoil will be simulated, particularly the Reynolds number. A 
particular feature of the rig will be the ability to span the Reynolds number 
range from conventional gas-turbine levels up to SSME levels. Figure 4 is a 
photograph of the test section, instrumented and awaiting installation in the 
facility. Testing will be underway in the near future. 


H5(g). UNSTEADY HEAT TRANSFER ANALYSIS 

In this task the objective is to formulate analytic models to aid in the 
prediction of the unsteady aerothermal loads expected in the SSME turbopump 
turbines . 

The ability to predict the aerothermal loads in the turbine unsteady flow 
environment is an essential input to the ability to predict life and durabili- 
ty. There are two aspects to the unsteadiness which occurs on turbine blades. 
The first is due to free-stream turbulence and wakes shed from upstream vanes 
and blades. The second is due to laminar to turbulent boundary layer transi- 
tion which occurs some distance downstream from the blade leading edge. Both 
of these issues are being addressed in this research program. 
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For the unsteady wake interaction problem, a viscous-inviscid boundary 
layer theory has been used to model the effects of passing wakes and free- 
streain turbulence on the heat transfer to the SSME turbine airfoils, particu- 
larly near the stagnation region. For the boundary layer transition experi- 
ments, an experimental approach is being used to study the effects of concave 
curvature, which occurs on the pressure side of the SSME turbine airfoil, dur- 
ing the transition process. 


H5(h). IMPROVED AERO AND HEAT TRANSFER PREDICTION 

In this task the objective is to improve the calculation of the aerody- 
namic environment in low-aspect-ratio turbines and to experimentally obtain 
the unsteady surface heat transfer response to a temperature disturbance. 

Knowledge of the local flow environment of turbine vanes and blades is es- 
sential for the accurate prediction of turbine temperatures and heat-transfer 
rates. In addition, the response of the turbine blade surfaces to rapid chang- 
es in stagnation temperature is unknown. Improved steady and unsteady heat 
transfer prediction is required for improved life and for turbine efficiency. 

An existing quasi-3D computer code will be activated at Calspan and used to 
compute pressure distributions over stator and rotor surfaces, substantially 
reducing the uncertainty in the aerodynamic environment and providing a relia- 
ble base for heat transfer predictions. The experimental phase will use a 
low-aspect-ratio turbine instrumented with thin-film heat flux sensors located 
on the vanes and rotor blades. The temperature level, disturbance, blade pass- 
ing frequency, and Reynolds number will be varied. A detailed description of 
the latest results under this task will be presented in another paper at this 
meeting. 


SUMMARY 

Significant progress is being made in the Aerothermodynamic Loads Defini- 
tion Study, particularly in the turbomachinery area, where there is a distinct 
overlap between the AETO efforts and research in the aeronautical gas turbine 
field. The continuing work under this program will provide the knowledge, 
understanding, and design methodology that will allow the development of ad- 
vanced Earth-to-orbit propulsion systems with high performance, extended serv- 
ice life, automated operations, and diagnostics for in-flight health 
monitoring. 
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PLANAR REACTING SHEAR LAYER 



SSME AERO-THERMODYNAMIC LOADS 
FLUID AND GAS DYNAMICS WORKING GROUP 
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Figure 2 
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HEAT-TRANSFER AND SURFACE-PRESSURE MEASUREMENTS 
FOR THE SSME FUEL -SIDE TURBOPUMP TURBINE* 

by 43^ -5/’ 

Michael G. Dunn / 3 3 3 

Calspan-UB Research Center 
Buffalo, New York 14225 


INTRODUCTION 

An experimental research program is currently underway at the Calspan-UB 
Research Center (CUBRC) that is designed to obtain detailed heat-flux and surface- 
pressure distribution on the first-stage vane and first-stage blade of the SSME fuel- 
side turbo pump turbine. This particular turbine is a two-stage machine and both stages 
w 11 be in use. However, at the present time, the first-stage vane, the first-stage 
bade, and the second sta;,e vane will be instrumented. The specific turbine being 
utilized is a combination of actual engine hardware and reproduced hardware consistent 
with that being used at NASA Marshall Space Flight Center for the initial measurements 
in their newly constructed blow-down turbine test facility. 

The facility being used at CUBRC is also a blow-down type facility, but it is of 
the short-duration shock-tunnel variety. The short-duration nature of the shock-tunnel 
facility permits use of thin-film thermometers which are used to measure the surface 
temperature histories at prescribed locations on the turbine component parts. Heat- 
flux 'alues are then inferred from the temperature histories using standard data reduction 
procedures. Miniature surface mounted pressure transducers are also used on both the 
vane and the rotating blade to obtain the desired pressure distributions. 


*The research described in this paper was supported by NASA Lewis Research Center 
under Grant No. NAG 3-581. 
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EXPERIMENTAL APPARATUS 


The experimental apparatus sketched in Figure 1 consists of a 0.47- m (18. 5-in) 
i.d. helium-driven shock tube with a 12.2-m (40-ft.) long driver tube and an 18.2-rn 
(60-ft.) long driven tube as a short-duration source of heated air, supplying the test 
section device located near the exit of the primary shock-tunnel nozzle. The 2.7-m 
(9-ft.) i.d. receiver tank is initially evacuated to a pressure on the order of I torr in 
ordfc ' to minimize the flow establishment characteristics of the model. The useful test 
time in this facility for the test conditions to be used for the S5ME measurements is 
on t le order of 35 rns, which is nearly two times greater than the test time available 
for the facility used in the previous CUBRC full-stage turbine studies (References 1-6). 
A more detailed photograph of the test section device housing the two-stage turbine 
is shown n Figure 2. This device consists of a forward transition section with a 
circular opening facing the supersonic primary nozzle flow. This transition section is 
followed by the 360 deg. annular passage containing the turbine stages. Downstream 
of the second rotor exit is an annular passage with a contoured nozzle at the exit. 
This nozzle is used to establish the desired pressure ratio across the turbine. The 
inten al model configuration duplicates the SSME turbopump configuration including the 
twelve upstream struts and the protruding bolt head on the dome. In our configuration, 
this forward dome houses an air motor that is used to accelerate the turbine from test 
to me desired speed just prior to initiation of the experiment and prior to arrival of 
the test gas. The aft internal cavity houses a 200-channel slip ring that is used to 
trai sfer the pressure and heat-transfer data from the rotor blade to the laboratory 
data recorders. 

Figures 3 and 4 are front and rear photographs of the first stage nozzle guide 
vane ring. A contoured leading-edge insert containing ten thin-film heat-flux gages 
will be installed at midspan on vane No. 32. Button-type thin-fim gages will be installed 
on vane No. 16 (pressure surface) and vane No. 17 (suction surface) at 10%, 50%, and 
90% s[ an. Miniature surface mounted pressure transducers will be installed at 10% 

span on vane No. 13 (pressure surface) and vane No. 14 (suction surface), at 50% on 

vane No. 23 (pressure surface) and vane No. 24 (suction surface), and at 90%span on 

vane No. 29 (press ire surface) and vane No. 30 (suction surface). The specific vanes 

and instrumenttion locations on those vanes at which pressure measurements will be 
performed are consistent with the Marshall instrumentation plan. The heat-flux 
measurements will be performed at locations such that direct comparisons can be made 
between the heat-flux and surface-pressure resjlts. 
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ORIGINAL PAGE IS 
OF POOR QUALITY 

Figure: 4 illustrates removal of a portion of the vane near the hub wail. This cut 
in vane profile is consistent with flight engine hardware. Several heat-flux gag^s have 
been placed near this cut on both the pressure and suction surfaces and on the hub 
endwall in order to investigate the influence of the cut on the flow field. 


Figure 5 is a photograph of one of the first-stage rotor blades that will be 
instrumented with flush-mounted miniature pressure transducers and thin-fiim heat-flux 
gages. Both heat-flux and pressure instrumentation will be placed at selected chordwise 
locations at 10%, 0% and 90% spans. The heat- flux instrumentation will be composed 

of button-type gages and a leading-edge insert. 

The measurements in and along the flowpath are important for determining the 
gas dynamics into and out of the turbine stage. The turbine model shown in Figure 2 
has approximately twenty static pressure measurements distributed between the inner 
..nd outer portions of the model flow path. Static pressure measurements of particular 
mportance are those at the inlet to the guide vane row, the intra stage measurements, 
and the static pressure measurements downstream of the second rotor. A spanwise 
total pressure rake is placed downstream of the second rotor and an attempt is being 
made to fit such a rake into the flow just downstream of the first rotor. Further, an 
attempt is being made to install a spanwise rake of total temperature probes downstream 
0 f the first rotor. The rakes just described may result in excessive flow blockage at 
the location of interest, and if so would be moved to a station downstream of the 
second rotor. 

EXPERIMENTAL CONDITIONS 


The SSME turbopump operates with a test gas composed of hydrogen plus steam. 
The facility to be used at CUBRC will operate with air as the test gas, but the 
appropriate scaling has b< en done to achieve the corrected conditions at the turbine 
inlet. It is import int that the experiment duplicate the design flow function ), 

the corrected speed (Np^y/V^), the wall to total temperature ratio (T w /T 0 ), and the 
stage total to static pressure ratio ( p 0 ^in^ P s^out^ ** * s a * s0 im P ortant to produce a 
turbine inlet Reynolds number sufficiently close to the operating value that the inlet 
boundary-layer characteristics are representative. 
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The SSME turbopump turbine measurements will be performed in a large turbine 
test facility at the Calspan Corporation. Because of the large dimensions of this 
facility it is possible to plac^ the inlet to the model housing the two-stage turbine at 
a couple of different locations in the nozzle expansion (see Figure 1), and in this 
manner, change the turbine inlet total pressure and Reynolds number while holding the 
turbine inlet total temperature constant. 

The nominal flow conditions for the SSME fuel side turbopump are; T 0 = 1890°R, 
P 0 = 5526 psia, inlet Mach No. = 0.14, inlet density = 4.29 lbm/ft 3 (H 2 + steam), Flow 
Function = 2.28, and overall total to static pressure ratio = 1.52. 

For the proposed measurements, air would be used as the test gas instead of 
hydrogen plus steam. Two experimental test conditions will be run, both at a total 
temperature of the turbine inlet gas equal to 1000°R and both at the appropriate flow 
function, overall total to static pressure ratio, and corrected speed. The two test 
conditions would provide two different turbine inlet total pressure values and Reynolds 
number values as noted in Table i. 


Table 1 

TEST CONDITIONS WITH AIR AS TEST GAS 


Test 

Condition 

Shock-Tube Conditions 

Model Conditions 
Ahead of Vane 

Turbine Conditions 


Po 

psia 

T 

*0 

°R 

^*)throat 

Ibm/sec 

P o 

psia 

T 

*0 

°R 

Re 

/ft 

w 

Ibm/sec 

Rotor 

Speed 

rpm 

T w / T o 

//I 

noo 

1000 

523 

270 

1000 

IxlO 7 



0.53 

//2 

2500 

1000 

1188 

615 

' 1000 

2xl0 7 



0.53 


CONCLUSIONS 

The current status of the measurement program described herein is as follows: 
(a) The model to house the two stage turbine has been designed and detailed engineering 
drawings have been prepared. Construction of the model has been initiated and is 
planned for completion in July 1989, (b) the pressure and heat-flux gages have been 
constructed and calibrated. Installation of this instrumentation in the components shown 
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in Figures 3 and 5 is scheduled for completion in June 1989, and (c) upon delivery of 
the model components, the instrumented vanes and blades will be installed and the 
measurement program will begin in the Summer of 1989. 
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Figure 3 FRONT-VIEW PHOTOGRAPH OF FIRST-STAGE VANE 



Figure 4 REAR-VIEW PHOTOGRAPH OF FIRST-STAGE VANE 
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Figure 5 PHOTOGRAPH OF FIRST-STAGE BLADE 
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A description of the Axial Flow Turbine Research Facility (AFTRF) 

>eing built at the Turbomachinery Laboratory of the Pennsylvania State 
University is presented in this report. The purpose of the research to be 
performed in this facility is to obtain a better understandir g of 
the rotor/stator interaction, three dimensional viscous flow field in 
nozzle and rotor blade passages, spanwise mixing and losses in these blade 
rows, transport of wake through the rotor passage, and unsteady 
aerodynamics and heat transfer of rotor blade row. The experimental 
results will directly feed and support the analytical and the computational 
tool development. This large scale low spared facility Is heavily 
instrumented with pressure and temperature probes and has provision for 
flow visualization and laser doppler anemometer measurement. The facility 
design permits extensive use of the high frequency response instrumentat ion 
on the station, ry vanes and more Importantly on the rotating blades. 
Furthermore, i facilitates detailed nozzle wake, rotor wake and boundary 
layer sw veys. The large size of the rig also has the advantage of 
operating at Reynolds numbers representat i ve of the engine environment. 


^39-3 V s 
/f88(p 


1. 1 RIEF OUTLINE OF RESEARCH PROGRAM 

The objective of the proposed research is to gain a thorough 
understanding of the fl)w field in a turbine stage, including three 
dimensional inviscld ani viscid effects, unsteady flow field and heat 
transfer in the rotor passage, three dimensional mixing in inter and 
ntra-blade passages, tip clearance effects, effects of film and convective 
cooling on flow and thermal field. A brief outline of the research 
>rograms proposed for the near future is given below: 

Recent investigations have shown that three dimensional effects have 
considerable influence on the performance, and heat transfer in turbine 
stages. A knowledge of the character istics of the three dimensional flow 
field nside the turbine- rotors, including the blade boundary layers, 
endwali flows and the inviscld velocities, is essential for the analysis 
and the design of turbines used in military and civilian aircraft. Viscous 
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and turbulence effects play a very significant role in the study of 
improved design, better efficiency, cooling air requirements, mechanical 
reliability, flow induced vibration, and flutter of the turbine blades* 
Because of the complicated nature of the problem, the three dimensional 
viscous, turbulent, unsteady flow field inside the rotors Is one of the 
least understood phenomena In turbines. 

The Pennsylvania State University research group has recently 
completed a detailed experimental and computational study of the three 
dimensional flow field In the axial flow compressors and this has provided 
valuable information on the endwall flow characteristics, blade boundary 
layers, tip leakage flow and the inviscid flow field (including blade 
pressure distribution) (example Refs. 1-3)- The boundary layers developing 
on turbines are also very complex; these are three dimensional with 
laminar, transitional, turbulent and separated zones. The flow field is a 
function of several parameters such as the solidity, incidence, camber, 
radial and cho^dwise pressure gradients, speed of rotation, blade geometry 
and the Mach n amber. Because of these complications, very few have 
attempted measurement of the rotor flow field. The objective of this 
Investigation Is to measure the entire flow field Inside a turbine rotor 
and nozzle at the design condition. An additional objective Is to develop 
analysis and computational techniques for the prediction of this flow field. 
The laser doppler velocimeter system ( non- intrusive) Is most desirable for 
flow measurements inside the rotor of the proposed facility. The wall and 
blade shear layers will be measured using hotwire sensors. The measurement 
will be taken at sufficient chordwise, radial, and tangential locations so 
as to derive a continuous variation of the flow field (including the blade 
and the annulus wall boundary layers) across the entire passage (hub to 
tip, leading to trailing edge, and pressure to suction surfaces). These 
measurements will provide information on the transition, turbulence 
properties, three dimensional flow velocities, and blade pressure 
distribution. 

In addition, the unsteadiness into rotor inlet will be resolved 
through detailed measurement of nozzle wake characteristics. This will be 
carried out 50 and 20 percent (of nozzle chord) spacing between the nozzle 
and the rotor. The nozzle wake survey at several radial and axial 
locations, including stations close to the nozzle trailing and rotor 
leading edge will be carried out. 

All the Investigations carried out hitherto concern an Independent 
study of steady and unsteady aerodynamics, or steady and unsteady heat 
transfer. These are usually from different set of facilities and operating 
conditions. Since heat transfer and aerodynamic phenomena are closely 
coupled, this approach has not led to a basic understanding of the flow and 
thermal field. One of the major objectives of the Penn State program is to 
carry out both the aerodynamic and heat transfer measurements in the same 
facility with identical operating conditions. The program for the 
measurement of steady aerodynamics was described earlier. The unsteady 
aerodynamics, including unsteady boundary layer properties will be measured 
using the hot wire techniques developed at the Pennsylvania State 
University. It is also planned to carry out unsteady convective heat 
transfer research in the Axial Flow Turbine Research Facility at a later 
date. The objective Is to model the turbulent heat fluxes near the solid 
walls in complicated viscous flow zones such as the areas affected by wake 
passing. During the Investigations, unsteady features of wall heat flux, 
stagnation temperature and wall temperatures will be investigated 
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experimentally and interpreted. The vane-blade interaction phenomenon will 
be captured by using thin film platinum wall heat flux sensors located in 
the areas where the wakes from the upstream rows impinge. These sensors 
are excellent indicators of phenomena such as laminar to turbulent 
transition, separation, reattachment etc. Their time response is 
sufficiently high to capture the unsteady features of the vane-blade 
interaction process. A new method of obtaining unsteady heat transfer 
information from the large scale cold turbine rig is discussed in detail in 
ref. 4. The measurement techniques are described in refs. 5 and 6. 


2. FACILITY DESCRIPTION 

A complete description of the facility and instrumentation is given in 
ref. 4. Only a brief summary is provided below. The Axial Flow Turbine 
Research Facility of the Pennsylvania State University is an open circuit 
facility 91. A cm in diameter and a hub to tip radius ratio of 0.73, with 
advanced axial turbine blading configurations. The facility consists of a 
large bellmouth inlet, a turbulence generating grid section, followed by a 
test section with a nozzle vane row and a rotor as shown in figure 1 . 

There are 23 nozzle guide vanes and 29 rotor blades followed by outlet 
guide vanes. Provisions exist for changing the vane-blade axial spacing 
from 20 to 50 percent of chord. The belmouth inlet is housed in an 
enclosure (not shown) covered with wire mesh and a thin layer of rubber 
foam to filter the air prior to entry to the inlet. 

A variable through flow is provided by two auxiliary, adjustable 
pitch, axial flow fans and an aerodynamically designed throttle. This 
system allows accurate control of the mass flow through the experimental 
stage up to a maximum of 22,000 cfm. The two fans in series produce a 
pressure rise of 74.7 mm Hg (40" of water) with a mass flow of 10.4 m3 per 
second under nominal operating conditions. The power generated by the 
experimental turbine rotor assembly Is absorbed by an eddy-current brake 
which is capable of absorbing up to 90 Hp. The speed of the rotor can be 
varied between 175 and 1695 RPM with the "dynamatic-adjusto speed" control 
system and can be held constant to + 1 RPM, with normal fluctuations in 
line voltage. The eddy current brake is cooled by a closed loop chilled 
water cooling system. 

The rotor and nozzle vane passages are instrumented with high 
frequency response instrumentation to measure steady (time averaged) and 
unsteady pressures and shear stresses. The details of the instrumentation 
used on the nozzle vane, rotor blade, nozzle casing, rotor hub, nozzle hub 
will be described in the next section. Provision has been made for a laser 
window (item B in Fig. 1) for LDV measurement of flow field upstream of the 
nozzle, nozzle passage, spacing between the rotor and the nozzle, rotor 
passage and downstream of the rotor passage. 

The facility is equipped with two traversing mechanisms. One of the 
probe traverse unit is mounted directly behind the rotor disk and has 
provisions for the radial and circumf erential traverses in the rotating 
frame. It is driven by a stepping motor driven by a computer controlled 
Indexer at tangential increments of 0.019 degrees/step to allow accurate 
measurement of c.he rotor wakes. The intra blade radial and tangential 
surveys are accomplished with the use of the L.C. Smith stationary probe 
traverse unit. This mechanism is also driven by a computer controlled 
indexer and stepping motor In 0.05 mm steps and 0.0234 degree steps in the 
radial and circumferential directions respectively. 
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The rotating to stationary interface data transmission system, 
attached to the rotor shaft ahead of the nose cone, is an integral part of 
the facility. It consists of a 150 ring mechanical (brush/coin type) slip 
ring unit, and a specialized ten-channel low noise/signal ratio mercury 
slip ring unit. A 32 channel electronic pressure scanner unit is located 
in the rotating drum downstream of the turbine rotor. The electrical 
signals carrying the pressure information is carried to the stationary 
frame through the slip ring assembly. The rotor frequency will be 
accurately determined by using an infrared emitter/receiver circuit located 
on the casing of the turbine rotor. This device generates the reflected 
infrared emissions from the tip of a selected rotor blade and directly 
provides the rotor frequency (once per revolution pulse-OPR). 

A data processing system is built around a microcomputer with a clock 
rate of 10MHz. The system consists of a 16 bit computer with 1MB random 
access memory, a disk operating system, ^0 Mb hard disk storage space, 
printer and plccter. All of the data from both stationary and rotating 
i nstrumer tatior can be processed on-line. One of the long range goals of 
the turbine research is to acquire unsteady heat transfer and aerodynamics 
data simultaneously. 

A unique 'eature of Axial Flow Turbine Research Facility is its 
ability to generate the aerodynamic and unsteady heat transfer data 
simultaneously. This feature allows the generation of benchmark quality 
lata to validate three dimensional viscous codes, which includes heat 
transfer. It is expected that this integrated approach will improve the 
neat transfer predictic l capabilities, especially in the rotor stator 
interaction region. 

The van* and blade design was carried out by General Electric 
Company's Aircraft Engine group. The aerodynamic design, while not 
representing any specific current or future GE product, does embody modern 
design philosophy. Stage loading, flow coefficient, reaction, aspect 
ratio, blade turning angles are all within the range of current design 
practice. State-of-the-art quasi three dimensional design methods are used 
to design airfoil shapes. The details of the design are incorporated in 
reference [7]. The blade geometry is shown in figure 2. Absolute Mach 
number at the inlet of the rotor changes from 0.30 to 0.23 from hub to tip. 
Exit absolute Mach number is in a range from 0.11 to 0.10 from hub to tip. 
Figure 3 shows a rotor blade and a nozzle guide vane after final machining. 

The design features of the blading and facility are given in Table 1. 
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TABLE 1 


OVERALL PERFORMANCE PARAMETERS 


Total temperature at inlet 

TTO 

520 °R 

Total pressure at inlet 

PTO 

I4.7psia 

Mass flow rate 

W 

24.32 lbm/s 

Specific work output 

AH 

0 Btu.s 

2 ' 36 TbiT 

Rotational speed 

N 

1300 rpm 

Total-to-total isentropic 
efficiency 


0.8930 

Total pressure ratio 

PT0/PT2 

1 .0778 

Pressure drop 

PT02-PT01 

30" of water 
(56.04 mm Hg) 

Pitchline reaction 

RXP 

0.3820 

Power 

HP 

81 . 2hp 

Torque 

T 

328.1 lbf.ft 

Stator efficiency 

n 3 

0.9421 

Rotor efficiency 

n r 

0.8815 

Mean swirl angle at exit of rotor 
(meas. from tangential dir.) 

“2 

29.97° 

No. of stators 

n v 

23 

No. of rotors 

n R 

29 


^ actual exit k.e. 

where n = — — : 7T~, 

s Ideal exit k.e. 

actual relati/e exit k.e. 


n r ideal relative exit k.e. 


, . ideal 

static enthalpy drop across 

rotor 

reaction — . _ 

ideal 

static enthalpy drop across 

stage 
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• 3. FACILITY INSTRUMENTATION 

The Axial Flow Turbine Research facility is equipped with large number 
of static pressure holes (nearly 500) at carefully selected locations on 
the blade, vane, casing and hub walls. Static pressure measurements in the 
rotating frame will be '.aken by a 32 channel electronic pressure scanner 
located in the rotor fr ime. The signal output, after multiplexing in the 
rotor frame, will be passed through a slip ring unit for processing in a 
computer controlled data acquisition system. 

To measure the shear stress on the walls and the blades, shear stress 
sensors are used in the turbine rig. These are deposited thin film, (V 
array) "glue on dual probes" manufactured by Dantec company. The 
construction is similar to those used by McCroskey (Ref. 8) and is shown in 
Fig. i). A total of 23 dual element sensors are mounted around the rotor 
blade at mid span, with close spacing near the leading edge of the blade. 

In addition, five sensors will be located on the hub (endwall) of the rotor 
blade passage. The data from these shear stress sensors will oe acquired 
using a P. c I 156 channel data acquisition system (20 KHz). Unsteady shear 
stresses, which are extremely important quantities In the study of 
vane-blade interaction, will be recorded iy a Kinetic Systems (CAMAC) high 
speed data acquisition system (1 MHz). 

The implementation of the dynamic pressure transducers in the turbine 
rig was driven by space limitations. The pressure transducers are inserted 
Into chambers, which in turn are connected to the turbine airfoil surfaces 
through 0.5 - 0.8 mm diameter holes. The design objective is to achieve a 
frequency response of 40 KHz. The sensors used are Kulite model XCS-093 
with a pressure range of maximum 5 psia. They are capable of measuring 
pressure fluctuations to an accuracy of 0.01 psia. Details of a dynamic 
pressure transducer used in the Turbine Facility are given in Figure 5. 
Sixteen Kulite transducers are located along the chord at the mid span of 
the rotor blade. Seven transducers will be on the pressure side and the 
remaining nine will be located on the suction surface of the next blade in 
the sane passage. A schematic showing the location of dynamic pressure 
transducers on the blade at midspan (sixteen) and on the hub (endwall) 
surfa e (five) is shown in Fig. 6. The low level signals from the dynanic 
press nre transducers are amplified In the rotating frame by using miniature 
Instrumentation amplifiers. The amplifiers will rotate in the rotor frame 
and provide a high level signal output before the signal reaches the slip 
ring unit. These amplifiers are located inside the instrumentation drum 
shown in Fig. 7. It is expected that this amplification may help to 
improve the quality of output signal from the dynamic pressure transducers. 

A preliminary static calibration study of an XCd-093 transducer was 
performed for implementation in the turbine rig. 

The LbV system available in the Turbomachinery Laboratory is a two 
channel fr nge type laser anemometer (LDA ) with on-axis back scatter light 
collection (TSI model 9100-6). The system operates with a 4 W Lexel 
argon-ion aser tuned for the green line (514.5 nm). The entire system 
comprising the laser, the transmitting and receiving optics is mounted on a 
table. Th ; table can be moved along and perpendicular to the optical axis 
and can al ;o be tilted in the vertical plane. The assembly is fitted on a 
hydraulic bench capable of moving in the vertical direction. The three 
degrees of freedom with the possibility to tilt the table helps to position 
the bisector of the beam in the radial direction. The LDV system will be 
used to measure must of the flow field, except the nearwall regions. 
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Provision has been made for hot wire measurement in both rotating and 
stationary system. Both two sensor and three sensor hot wire probes will 
be used in a stationary frame of reference to measure the boundary layers 
on stationary walls (on both the casing and hub walls) and the nozzle vane 
surf ace. 

The mechanism used for traversing the rotating hot wire probe consists 
of both the circumferential and the radial traversing units. The hot wire 
probe location can be incrementally changed by using 0.019 degree steps, in 
circumferential direction. This motion is also checked and corrected for 
possible errors by an angular encoder device. The planetary gearbox 
providing this precision motion is driven by a stepping motor having a 
25000 steps per revolution. Both the radial and circumferential traversing 
units will be simultaneously controlled by a computer based system. 


I). DATA ACQUISITION AND PROCESSING 

The Kulite, skin friction, heat transfer, hotwire and static pressure 
data from the rotor will be transmitted through the rotating drum to a slip 
ring unit (figure 1). The rotating drum consists of two probe traversing 
mechanisms for both the radial and the circumferential traverse, 4 constant 
temperature anemometers, a 32 channel electronic pressure 
scanner/transducer unit and a 30 channel data acquisition system for 
voltage inputs. The data acquisition system is shown in Fig. 8. 

The slip ring unit is of brush type and has 150 channels. The slip 
ring unit is housed in a cowling in front of the facility (Fig. 1). The 
unit has a rigid one piece housing eliminating bearing alignment problems. 
Each ring carries four brushes made of silver graphite. The rings are made 
up of coin silver which withstands up to current levels of 5 Amps. The 
brushes are individually removable and replaceable. The contact resistance 
is about 5 mill i. ohms maximum. 

A medium speed data acquisition system with 96 channels of pressure 
transducers having a 20,000 samples per second capacity is ava lable In the 
laboratory. The PS I model 780-B pressure measurement system is a fully 
integrated test instrument which offers a transducer per port for 
multipressure measurement applications. The present system is also capable 
of using an additional 60 channels which may be described as "electrical 
input" channels for signal types such as thermocouples, strain gauges, 

RTD’s etc. This is a complete data acquisition system with a special 
emphasis on pressure measurement. Total system inaccuracy for pressure 
measurements including sensors is in the order of ± 0.10 percent full scale 
for the worst case. However, a typical value Is in the order of ± 0.07 
percent. The system is interfaced to an IBM compatible personal computer 
through an IEEE -488 (GPIB) with input output rates up to 30,000 bytes per 
second. 

For high speed data acquisition, a 16 channel, 12 bit transient 
recorder with a 1 Megabyte solid state memory is planned (Kinetic Systems 
Int., 4022-DIA-4050-23A) . This system will be built in a CAMAC crate. The 
acquisition system will be controlled by a parallel bus crate controller 
interfaced to an IBM PC-AT computer. With the current configuration the 
system will be able to store 500,000 measurements for a maximum duration of 
one second either In its own memory or in the computer storage medium. The 
main use of this system will be in the area of hotwire measurements, wall 
shear stress measurements, Kulite dynamic pressure measurements. 
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Figure 3 

A Nozzle Guide Vane and a Rotor Blade 
After Final Machining 
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Figure 8 


Measurement Chain and Data Acquisition System 
in the Stationary and Rotor Frame of the 
Pennsylvania State University AFTRF 
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A quasi-three-dimensional analysis has been developed for unsteady rotor- 
stator interaction in turbomachinery. The analysis solves the unsteady Euler 
or thin-layer Navier-Stokes equations in a body-fitted coordinate system. It 
accounts for the effects of rotation, radius change, and stream-surface thick- 
ness. The Baldwin-Lomax eddy-viscosity model is used for turbulent flows. 

The equations are integrated in time using an explicit four-stage Runge-Kutta 
scheme with a constant time step. Implicit residual smoothing is used to in- 
crease the stability limit of the time-accurate computations. The scheme is 
described, and stability and accuracy analyses are given. 

Results are shown for the stage of the space shuttle main engine high 
pressure fuel turbopump. Two stators and three rotors were used to model the 
41:63 blade ratio of the actual machine. Implicit residual smoothing was used 
to increase the time step limit of the explicit scheme by a factor of six with 
negligible differences in the unsteady results. We feel that the implicitly 
smoothed Runge-Kutta scheme is easily competitive with implicit schemes for 
unsteady flows whi.e retaining the simplicity of an explicit scheme. 
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UNSTEADY ROTOR-STATOR INTERACTION CODE 

DEVELOPED BY P. C. E. JORGENSON R. V. CIIIMA 


DESCRIPTION 

• UNSTEADY THIN-LAYER NAVIEIl-STOKES SOLVER 
FOR ROTOR-STATOR INTERACTION 


FEATURES 

• QUASI-3-D FORMULATION 

INCLUDES ROTATION, RADIUS CHANGE, & STREAM SURFACE THICKNESS 

• C-TYPE GRIDS OVERLAP AT INTERFACE 

• MULTI PASSAGE CAPABILITY 

• EXPLICIT RUNGE-KUTTA SCHEME + IMPLICIT RESIDUAL SMOOTHING 
UNCONDITIONALLY STABLE (PRACTICAL LIMIT CFL = 0(20)) 

2nd ORDER TIME ACCURATE 

FIRST USE OF IMPLICIT SMOOTHING FOR UNSTEADY PROBLEMS 
RESULTS 

• SPACE SHUTTLE MAIN ENGINE TURBINE 


CD-89-39438 


QUASI-THREE-DIMENSIONAL STREAM SURFACE 


ROTATION 
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QUASI-3-D THIN LAYER NAVIER-STOKES EQUATIONS 
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- CONTRAVARIANT VELOCITY COMPONENTS IN RELATIVE SYSTEM 


w e - v 0 - rn - RELATIVE TANGENTIAL VELOCITY 

0 - BLADE ROTATION SPEED 
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K 2 - (pVg + p-Re o 2 2 ) ^ m r/r + (p - Re O 33 ) d m h/h 

- CENTRIFUGAL FORCE TERM + p 6 m (AREA) 


Eq - Eg/r 

; 0-METRICS SCALED BY 1/r 

Hq - r^Jr 

1/J-rh/J : JACOBIAN SCALED BY AREA 
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ROTOR-STATOR INTERACTION CODE DETAILS 

INITIAL CONDITIONS 

ANALYTIC 1-D SOLUTION OF FLOW EQUATIONS WITH AREA C HANGE 

BOUND A RY CONDITIONS 

INLET 

P 0 , T 0 , AND v 0 SPECIFIED 

R' -V m -— EXTRAPOLATED 
y-1 

EXIT 

P SPECIFIED 

3 CONSERVATION VARIABLES EXTRAPOLATED 
WALLS 

INVISC1D FLOW - TANGENCY, p EXTRAPOLATED 
VISCOUS FLOW - NO-SLIP, T w SPECIFIED 
NORMAL MOMENTUM EQUATION FOR PRESSURE 

PERIODIC BOUNDARIES 

SOLVED LIKE INTERIOR POINTS 
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MULTISTAGE RUNGE-KUTTA ALGORITHM 


GOVERNING EQUATIONS 
d t q = -J[R, - {R v +D)\ 

n, = INVISCID RESIDUAL 
Rv = VISCOUS RESIDUAL 
D = ARTIFICIAL DISSIPATION TERM 
MULTISTAGE SCHEME 


<7o — Qn 


< 7 i — ?o “ a i J [^/ <70 ~ [Rv + D) go] 


Qk = <?o “ ct k J&t j Rj q k -i - [Rv + D) g 0 ] 

<7n+l = qk 

R v kD EVALUATED AT FIRST STAGE ONLY 
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ARTIFICIAL DISSIPATION 

NONCONSERVATIVE VERSION OF JAMESON FORMULATION 
Dq = ( D i -f Dr,) q 

^ -DIRECTION OPERATOR 
D^q = C (V 2 g^ - V 4 qmt) 

WHERE: 

C = — - — 

JAlij 

(LOCAL A t iti MINIMIZES DISSIPATION) 

V 2 ~ fi 2 max (i/,* + v { „ i ) 

I 7 4 = max (0,/X4 — V2) 

AND 

_ \Jj±Uj Z 2P, >; + p i-i,A 
t,J lA+i,/ + 2f»,y + 

;x 2 = 0(1) 

;i4 =s O(^) 
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IMPLICIT RESIDUAL SMOOTHING 


USE A TIME STEP GREATER THAN THE STABILITY LIMIT 
MAINTAIN STABILITY BY SMOOTHING THE RESIDUAL IMPLICITLY 

REWRITE STAGE (A') OF MULTISTAGE SCHEME AS 
= q W - = - ak JM - ZV°>] 


IMPLICIT SMOOTHING STEP 


(1 - e<5«)(l - 


+ A «/*) 

UNCONDITIONAL! Y STABLE IF 


1 

6 ~ 4 


CFL 

CFL' 


- 1 


WHERE 

CFL* IS COURANT LIMIT OF THE UNSMOOTHED SCHEME 
CFL IS THE LARGER OPERATING COURANT NUMBER 
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IMPLICIT RESIDUAL SMOOTHING 


TWO-DIMENSIONAL LOCAL SMOOTHING PARAMETER 


c,j = max 





WHERE 

CFL * IS COURANT LIMIT OF THE UNSMOOTIIED SCHEME 
CFLij IS THE LOCAL COURANT NUMBER 
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IMPLICIT RESIDUAL SMOOTHING 


FEATURES 

• IRS DOES NOT EFFECT FORMAL TIME ACCURACY 
OF R-K SCHEME 

• SMOOTHING ONLY OCCURS WHERE LOCAL COURANT NUMBER CFL h 
IS GREATER THAN THE COURANT LIMIT OF THE R-K SCHEME CFL' 
(EG - VISCOUS REGION) 

. CAN USE CFL = 0(10 - 20) IN VISCOUS REGIONS WHILE CFL = O(.l) 

IN INVISCID CORE 

• £ij EASY TO COMPUTE, OR CAN BE STORED 

t IRS REQUIRES SCALAR TRIDIAGONAL INVERSION ALONG EACH 
GRID LINE 

ADDS APPROXIMATELY 2G% TO CPU TIME BUT BOOSTS TIME STEP 
BY 652% FOR A NET DECREASE OF 460% IN CPU TIME 

• EXPLICIT R-K SCHEME WITH IRS COMPETETIVE WITH IMPLICIT 
SCHEMES FOR UNSTEADY PROBLEMS 
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SSME STATOR AND ROTOR GRIDS 
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ORIGINAL PAGE 

BLACK AND WHITE PHOTOGRAPH 


DETAIL OF STATOR/ROTOR GRID OVERLAP 



UNSTEADY LIFT COEFFICIENTS ON 
STATOR AND ROTOR 


STATOR 



ROTOR 


NUMBER OF PITCH ROTATIONS 
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PRESSURE DISTRIBUTION ENVELOPE ON STATOR 


PS/PO 

STATOR 



PRESSURE DISTRIBUTION ENVELOPE ON ROTOR 


PS/PO 

ROTOR 
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ABSOLUTE MACH NUMBER CONTOURS 
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ENTROPY CONTOURS 



CD-89-39453 


245 


ORIGINAL PAGE 

BLACK AND white photograph 



VELOCITY VECTORS AROUND ROTOR 



UNSTEADY ROTOR/STATOR INTERACTION CODE 


CONCLUDING REMARKS 

. TIME-ACCURATE ROTOR/STATOR INTERACTION CODE HAS DEEN DEVE 
OPED ' 


. CODE IS APPLICABLE TO AXIAL OR RADIAL MACHINES 


. VISCOUS RESULTS FOR SSME TURBINE SHOW 
UNSTEADY LOADING it SEPARATION, it MIGRATION OF STATOR WAKES 


. IMPLICIT RESIDUAL SMOOTHING 
DOES NOT EFFECT FORMAL TIME ACCURACY OF R-K SCHEME 

ADDS 26% TO CPU TIME 

INCREASED TIME STEP BY A FACTOR OF 0.0 

RESULTS IN A N ET DECREASE OF ABOUT 150% IN CPU TIME 

MAKES EXPLICIT R-K SCHEME COMPETETIVE WITH IMPLICIT SCHEMES 
FOR UNSTEADY PROBLEMS 
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AVERAGE-PASSAGE FLOW MODEL DEVELOPMENT 




John J. Adamczyk 
NASA Lewis Research Center 

Mark L. Celestina, Tim A. Beach, Kevin Kirtley, 
Sverdrup Technology, Inc. 

Lewis Research Center Group 

and 

Mark Barnett 

United Technologies Research Center 


A research effort has been underway to develop a three-dimensional model 
for simulating multistage turbomachinery flows using today’s supercomputers. 
This model, referred to as the "average passage" flow model, describes the 
time-averaged flow field within a typical passage of a bladed wheel within a 
multistage configuration. To date, a number of inviscid simulations have been 
executed to assess the resolution capabilities of the model. Recently, the 
viscous terms associated with the "average passage" model have been incorpo- 
rated into the inviscid computer code along with an algebraic turbulence 
model. A simulation of a stage-and-one-half , low-speed turbine has been ex- 
ecuted. The results of this simulation, including a comparison with experimen- 
tal data, is the subject of this report. 

A goal of computational fluid dynamics for turbomachinery is the predic- 
tion of performance parameters and the flow processes which set their values. 
Achieving this goal for multistage machinery is made difficult by the wide 
range of length and time scales in the associated flow fields. Currently, the 
procedure used in the design and off-design analysis is based on a quasi-three- 
dimensional flow model whose origins can be traced to the late forties and 
early fifties. 

Although proven useful, this flow model has its limitations. Among these 
are the inability to analyze off-design performance and unconventional machin- 
ery where extrapolation of the underlying empirical data base is required. 

Other problems arise whenever there are large local variations in the radial 
velocity component within a blade passage. It is generally agreed that a way 
of overcoming these shortcomings is the development of true three-dimensional 
flow models. The "average passage" flow model under development at Lewis 
Research Center is such a model. The objective of this paper is to present 
the status of the development of this model for multistage turbines. This 
will include several comparisons with measurements obtained from a recently 
completed experimental program at the United Technologies Research Center. 
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SIMULATION RESULTS 


The simulation executed was the low-speed rotating rig at United Technolo- 
gies Research Center. The low-speed rotating rig (LSRR) is a stage-and-one- 
half turbine consisting of an inlet guide vane, a rotor, and a stator. The in- 
let guide vane contains 22 blades, and the rotor and stator both contain 
28 blades. The flow coefficient, <J), is 0.78, and the spacing between blades, 
B x , is 0.5. The LSRR grid contains 228 axial, 25 radial, and 41 circumferen- 
tial points. Each blade row contains 40 axial points distributed along the 
chord with 26 axial points between each blade row, the inlet and exit. 

The results presented required 11 hours of Cray 2 CPU time. They repre- 
sent but a small fraction of the information obtained from the simulation. 

They are intended to illustrate the degree to which one can quantitatively pre- 
dict performance parameters that are of interest to designers and to reveal 
qualitative information identifying flow phenomena that may have an effect on 
performance. These results also reflect the current state of model develop- 
ment. The first series of results shows the predicted pressure distribution 
on the surface of each blade row of the turbine as a function of axial chord 
length and percent of span height. The span locations measured from the hub 
are 1.3, 12.5, 50, 87.5, and 98.7 percent, respectively. The experimental 
measurements taken at these locations are also shown. Experimental data were 
also available for 25 and 75 percent of span, but were not used, since they 
provided little additional information relative to the current discussion. 

The results for the first vane are shown in figure 1. The predicted loading 
level is in good agreement with the measurements of Dring (1988). The predict- 
ed pressure-surface pressure distribution is in excellent agreement with the 
experimental results. For the suction surface the agreement between measure- 
ment and simulation is good for the region forward of the minimum pressure 
peak. Aft of the peak the agreement between experiment and simulation deterio- 
rates, This deterioration is believed to be related to viscous effects (i.e., 
turbulence and transition modeling) whose modeling could be improved. Some 
exploratory calculations suggest that the boundary layer aft of the suction- 
surface minimum pressure is growing too rapidly and, as a result of the radial 
pressure gradient, is being transported toward the hub to an extent greater 
than that suggested by flow visualization studies. Improvements in the agree- 
ment between simulation and experiment have been obtained by incorporating a 
simple transition model in which the flow remains laminar forward of the mini- 
mum pressure peak and Baldwin-Lomax turbulence model as implemented by Dawes 
(1986). 

Figure 2 shows the results for the first blade row, which incorporated 
these changes in the turbulence model. The improvement is obvious. The re- 
maining results for the rotor and second stator incorporated the modified tur- 
bulence model. Figure 3 shows the predicted and measured pressure distribu- 
tion for the rotor. The predicted loading levels appear to be in good agree- 
ment with measurements, with the exception of the hub and tip region. The 
present simulation does not include a clearance region, which should account 
for some of the discrepancy in the tip region. The pressure distribution 
along the pressure surface is once more in excellent agreement with the mea- 
surements. At the midspan and at 25 and 75 percent (not shown) of span, the 
predicted pressure distribution along the suction surface is in good agreement 
with the data. At 1.3 and 12.5 percent of span, the suction-surface pressure 
coefficient is lower than that measured. As a result, the loading is lower 
over the forward portion of the rotor than what has been measured. Although 
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the cause of this discrepancy is, at present, unknown, one could speculate 
that it may be due to an improper estimate of the magnitude and extent of the 
low momentum fluid exiting the first vane. 

The pressure distribution for the last vane is shown in figure 4. Once 
again, the loading level is well predicted, with the exception of the 1.3 per- 
cent of span location. The underpredicted suction-surface pressure coeffi- 
cient at 1.3 and 12.5 percent of span suggests that the flow incidence to 
these sections is underestimated. There also appears to be a shift of the pre- 
dicted pressure distribution relative to the measured distribution. This 
shift is believed to be caused by an overestimate of the loss generated by the 
first two blade rows. With the exception of this discrepancy, the pressure 
distribution on the pressure surface is in good agreement with measurements. 
Similarly, the predicted suction-surface distribution at midspan agrees well 
with the experimental distribution. 

Currently, the results of this simulation are being reduced to a format 
which will allow for the comparison with measurements of total pressure level 
and flow angle. This comparison should provide additional information for 
judging the accuracy of the simulation. 


SUMMARY AND CONCLUSION 

Given the early state of the average passage model development, the re- 
sults presented in this report are very encouraging. The amount of empirical 
information used in the stage-and-one-half turbine simulation is considerably 
less than that required to achieve comparable results using today's quasi- 
three-dimensional flow models. 
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FIRST BLADE PRESSURE DISTRIBUTIONS 

FIRST VANE 



X/BX 


BLADE PRESSURE DISTRIBUTIONS 

FIRST BLADE ROW 


PRESSURE 

COEFFICIENT 




O SUCTION SURFACE 
□ PRESSURE SURFACE 



Figure 2 
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BLADE PRESSURE DISTRIBUTIONS 

ROTOR 



BLADE PRESSURE DISTRIBUTIONS 

LAST VANE 


PRESSURE 

COEFFICIENT 






O SUCTION SURFACE 
□ PRESSURE SURFACE 



Figure 4 
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